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ABSTRACT 


This  thesis  presents  an  investigation  of  a  variety  of 
equation  of  state  forms  for  pure  components.  The  investiga¬ 
tion  was  undertaken  in  order  to  develop  an  equation  of  state 
which  would  adequately  describe,  not  only  the  P-V-T  behavior 
over  wide  ranges  of  pressure  and  temperature,  but  also  derived 
properties  such  as  enthalpy  departure  and  fugacity. 

The  best  equation  of  state  form  was  found  to  be  a 
fifteen  constant  expansion,  based  upon  a  modified  form  of  the 
Martin-Hou  equation.  The  equation  of  state  was  fitted,  over 
wide  ranges  of  pressure  and  temperature,  for  eight  gases,  with 
good  results.  The  equation  of  state  was  restricted  to  the 
vapor  phase  because  of  a  drastic  loss  of  accuracy  encountered 
when  trying  to  represent  two  phases  with  one  set  of  equation 
constants . 

A  method  of  obtaining  improved  thermodynamic  property 
^representation  is  presented.  The  method  showed  merit,  but 
was  hampered  by  a  lack  of  adequate  experimental  thermodynamic 
pisoperties.  The  usefulness  of  the  method,  however,  is  apparent. 

All  data  fitting  required  in  this  study  was  accomp¬ 
lished  using  the  Chebyshev  criterion  of  best  fit.  Examples 
of  the  data  fitting  procedure  are  presented  in  the  form  of 
computer  programs.  Explanation  of  the  use  of  the  computer 
programs  is  given. 
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INTRODUCTION 

The  motivation  for  this  investigation  of  equations 
of  state  was  supplied  by  the  development  of  an  alternate  con¬ 
cept  in  data  fitting.  This  fitting  technique  is  described  by 
P.K.  Leung (13)  and  frames  data  fitting  as  a  linear  programming 
problem.  The  criterion  of  best  fit  used  in  this  procedure  is 
minimization  of  the  maximum  absolute  or  absolute  percentage 
deviation.  This  is  generally  known  as  the  Chebyshev  criterion. 
Linear  programming  data  fitting  has  several  features  which 
make  it  especially  suitable  for  investigating  equations  of 
state. 

The  purpose  of  this  investigation  was  to  develop  some 
general  framework  for  an  equation  of  state  applicable  to  pure 
components.  To  be  completely  general,  the  equation  form  had 
to  be  applicable  to  all  types  of  gases  and  able  to  cover  as 
wide  a  range  of  pressures  and  temperatures  as  are  normally 
encountered  industrially. 

An  equation  of  state  must  satisfy  two  conditions  in 
order  to  be  applicable  to  a  particular  component.  These  con¬ 
ditions  are: 

1.  It  must  provide  an  accurate  representation  of  the  volu¬ 
metric  behavior  over  the  range  of  pressure  and  tempera¬ 
ture  to  be  fitted. 

2.  It  should  provide  a  reasonably  accurate  representation  of 
the  thermodynamic  properties  which  can  be  derived  from  the 
volumetric  behavior. 
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The  first  condition  can  usually  be  met  with  standard 
".least  squares"  fitting  techniques,  although  the  accuracy  of 
the  fit  depends  upon  the  range  of  variables  required  and  the 
complexity  of  the  approximating  function  used.  Meeting  the 
second  condition  is  not  so  easy.  Accurate  representation  of 
the  volumetric  properties  does  not  guarantee  equally  accurate 
representation  of  all  of  the  derived  thermodynamic  properties. 

An  analysis  of  the  relative  accuracy  of  calculated 
thermodynamic  quantities  was  made  in  order  to  obtain  some  in¬ 
sight  into  what  kind  of  accuracy  to  expect  for  the  various 
properties.  A  method  for  improving  the  accuracy  of  calculated 

thermodynamic  properties  was  developed  utilizing  the  linear 

♦ 

programming  data  fitting  procedure.  This  procedure,  however, 
requires  accurate  experimentally  determined  thermodynamic 
properties  which  are  consistent  with  a  set  of  P-V-T  data. 
Experimental  data  of  this  kind  are  generally  hard  to  find. 
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LINEAR  PROGRAMMING  DATA  FITTING 

A  detailed  explanation  of  the  linear  programming  data 
fitting  procedure  is  given  in  the  thesis  by  Leung (13).  It  is 
appropriate,  however,  to  give  a  brief  description  of  the  con¬ 
cept,  its  solution  and  advantages,  since  linear  programming 
data  fitting  was  used  almost  exclusively  in  this  thesis. 

Consider  the  approximating  function, 


a  .  g  .  (x  . ) 
3  3  i 


where 


g^ (x^)  are  specified  functions  of  x. 

Using  the  Chebyshev  criterion  of  best  fit,  the  follow¬ 
ing  problem  can  be  formulated. 

n 

| f (x . )  -  y  a.g.(x.)|  <  X  l<i<m 

1  j=l  3  1 

where  we  wi$h  to  minimize  X  maximum. 

The  absolute  value  restriction  may  be  removed  by  the 
formation  of  positive  and  negative  equations. 


n 

f  (x.  )  -  l  a  .g  .  (x.  )  <  X  1  <i  <m 

'  j=l  J  J 


n 

-  f(x.)  +  l  a.g.(x.)  >  -X 

j  =  l  3  J 


Rearranging  the  two  equations  gives, 
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n 

x  +  l  aug-iU..)  >  f(x.) 
j=l  3  3  1  1 


n 

X  -  T  a  .  g  .  ( x .  )  >■  -  f  ( x .  ) 
j£x  3yD  i  i 


The  data  fitting  problem  can  now  be  formulated  in  the 


following  linear  programming  form: 


Minimize  (X)  maximum,  subject  to; 


n 

1)  X  +  y  a.g.(x.)>f(x.)  l<i<m 

j±l  i  i 


n 

2)  X  -  l  a  .g .  (x  . )  >  -  f  (x  . ) 

j _-l  J  J  ^ 


With  the  introduction  of  slack  and  artificial  variables 


the  problem  becomes; 


m 

minimize,  Z  =  X  +  M  £  u. 

i=l  1 

where  M  is  a  penalty  coefficient,  subject  to; 


n 


+  J,  (yl.  -  y2.> 


j  =  l 


V,  +  u.  =  f  (x.  ) 

1 .  i  l 

l 


X  - 


n 

I 

j=i 


l  (y 


i . 
3 


y2j)  g j  Cx±)  + 


V, 


-f (X.) 


l<i<m 


(  :  «  '  ■  -V  " 


rnr  *  '  j  / '  XBe  I*  I  PfJ  .V,"  Xi. 
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where  V,  ,  V0  are  the  slack  variables  and  u.  are  the  arti- 

1.2.  i 

l  l 

ficial  variables.  The  variables  are  not  restricted  in 
sign  so  the  substitution 


ai  =  yl  “  y2 
J  D  D 


was  made.  Thus,  all  the  variables  X ,  y^  ,  y2  ,  ,  V2  and 

ij  ^  j  xi  i 


u^  are  non-negative. 

The  problem  formulated  can  be  solved  directly  by  the 
simplex  method  to  obtain  the  optimum  equation  coefficients  a^. 
However,  if  the  number  of  data  points  m  is  very  large,  the 
problem  quickly  becomes  too  large  for  practical  computation. 

A  means  of  improving  the  computational  efficiency  of  the  prob¬ 
lem  is  available  through  the  application  of  dual  linear  pro¬ 
gramming. 

The  rules  for  obtaining  the  dual  of  a  linear  programm¬ 
ing  problem  may  be  summarized  in  the  following  table  where; 

A  -  coefficient  matrix  for  the  constraint  equations. 

b  -  constraint  equation  values  for  the  primal  problem. 

x  -  variables  of  the  primal  problem. 

C  -  objective  function  coefficients  for  the  primal 
problem. 

Superscript  T  refers  to  the  transpose  of  the  matrix. 

y  -  variables  of  the  dual  problem. 
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Primal 


non-negative 


x  x 


unrestricted 

. x 
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C 
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k+r 

b,  * 
k+r 
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bk+r+l 

• 
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E 
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• 
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k+r+s 

i 

bk+r+s 

C, . c 

c  i . C 

min.  Z 

_i _ E 

_ E±i _ " 

Dual 


unrestricted  non-negative  non-negative 
^1 . ^k+1*  *  ^k+1  ^k+r+1  *  *  *^k+r+s 


at 

T 

-C 

T 

E 

< 

cx 

• 

• 

• 

C 
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bt 
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Cp+1 
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Use  of  these  tables  requires  that  all  b  values  in  the 
primal  be  adjusted  so  that  they  are  non-negative. 

The  duality  theorem  states  that  when  the  optimum  solu¬ 
tion  to  the  dual  problem  has  been  found, the  objective  function 
value  is  equal  to  that  of  the  optimum  solution  of  the  primal. 
The  values  of  the  dual  problem  variables  bear  no  direct  re¬ 
lationship  to  the  primal  variables.  The  primal  variables, 
however,  can  be  solved  for  by  solving  the  set  of  constraint 
equations  in  the  primal  which  correspond  to  non-zero  dual 
variable  values.  The  constraint  equations  become  equalities 
when  the  dual  variables  corresponding  to  those  equations  are 
non-zero.  The  only  non-zero  variables  in  the  linear  program 
solution  are  the  basic  variables.  Therefore,  the  constraint 
equations  corresponding  to  the  optimum  solution  basic  variables 
can  be  made  equality  equations  and  solved  for  the  optimum 
primal  variable  values. 

The  dual  of  the  data  fitting  problem  has  the  follow¬ 
ing  form. 

yl  •  •  •  •  ym  ym+l . y2m 


1 

1  .  .  .  .  1 

1 . 1 

< 

1 

• 

gl<xl>  •  '  gl(xm> 

~gl*xm+l)-  •  *  “gl(x2m) 

= 

0 

• 
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•  •  " 

— 
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• 

• 

•  • 
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•  • 

= 

0 

n 

g  (x,  )  .  .  g  (x  ) 

1  ^n  m 

_gn(xm+l)  -  "gn(x2m) 

= 

0 
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The  objective  function  Z  is  to  be  maximized.  All  of 
the  dual  variables  to  y ^ are  non-negative.  To  obtain  a 
starting  solution  a  slack  variable  is  added  to  equation  1  and 
artificial  variables  are  added  to  equations  2  to  n.  For  the 
non-trivial  case  where  A>0  the  slack  variable  introduced  in 
equation  1  must  take  a  value  of  zero  in  the  final  solution. 

To  insure  this,  equation  1  is  made  an  equality  equation  and 
an  artificial  rather  than  a  slack  variable  is  added  to  the 
equation. 

The  optimum  solution  of  the  dual  problem  gives  the 
optimum  value  of  A  and  the  variables  in  the  optimum  basic 
solution,  tell  at  which  data  points  the  maximum  A  value  occurs. 
The  constraint  equations  corresponding  to  the  basic  data 
points  can  then  be  made  equality  equations  and  solved  for  A 
and  the  equation  coefficients  a^.  This  requires  the  solution 
of  a  set  of  n+1  linear  equations. 

In  the  solution  of  a  linear  programming  problem  addi¬ 
tion  of  variables  costs  very  little  in  computation  time  while 
addition  of  constraints  adds  a  great  deal  to  computation 
time.  For  the  data  fitting  problem  the  number  of  primal  con¬ 
straints  is  2m,  while  for  the  dual  the  number  of  constraints 
is  n+1,  where  m  is  the  number  of  data  points  and  n  is  the 
number  of  equation  coefficients.  Thus,  solution  of  the  dual 
of  the  data  fitting  problem  represents  a  very  significant 
gain  in  computational  efficiency. 
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There  are  several  characteristics  of  the  linear  pro¬ 
gramming  method  of  data  fitting  which  are  worthy  of  note: 

1.  The  approximating  function  can  be  any  linear  combination 
of  specified  functions. 

2.  No  specific  data  point  spacing  is  required. 

3.  The  data  to  be  fitted  must  be  well  smoothed  and  error 
free  for  the  results  to  be  meaningful.  However,  large 
random  errors  can  be  detected  by  careful  examination  of 
the  data  points  which  control  the  fit.  Several  controlling 
data  points  occurring  within  a  small  range  of  the  variables 
indicates  a  questionable  data  point.  Data  points  in  such 

a  region  should  be  graphically  plotted  in  order  to  pin¬ 
point  a  possible  error. 

4.  There  are  no  computational  difficulties  encountered  in 
solving  for  the  equation  coefficients,  even  for  high 
order  equations.  The  method  does  not  suffer  from  the  ill- 
conditioning  problem  encountered  when  fitting  high  order 
equations  using  "least  squares". 

Various  constraints  on  the  function  fit  such  as  deriva¬ 
tives  of  the  function  can  readily  be  incorporated  into 
the  fitting  procedure. 


5. 
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B-W-R  EQUATION 


The  Benedict-Webb-Rubin  equation  of  state  has  gained 
the  widest  acceptance;  therefore,  it  was  chosen  as  a  conveni¬ 
ent  reference  in  evaluating  the  equations  of  state  studied  in 
this  thesis.  An  extensive  study  of  the  B-W-R  equation  for 
propane  was  undertaken  in  order  to  establish  the  fitting 
accuracy  of  the  equation.  Various  methods  of  fitting  were 
used,  including  methods  for  improving  enthalpy  departure  pre¬ 
dictions. 

The  following  results  establish  the  accuracy  of  the 
pressure  and  enthalpy  departure  predictions  of  the  B-W-R 
equation  for  propane.  They  also  indicate  the  utility  of  the 
various  fitting  criteria. 

Equations : 


A  C 

Z  =  1.0  +  (B - - - ^)d  +  (b 

°  RT  RT 

aad~*  cd2  0  J2 

+  -  +  — j  [(1  +  ya  )  e  yQ  ] 


— )  dJ 
RT 


RT 


RT' 


(H  -  Vt 


4C 

J  {  (B  RT  -  2A - 4 

o  O 


-)  d  +  (2bRT  - 


d 

3a)  — 
2 


6aad 

+  -  + 


cd2  3  -  3e  yd 


-Yd‘ 


[• 


T‘ 


Yd‘ 


+ 


y  d 


2 


]  > 


J 


conversion  factor  from  units  of  PV  to  BTU/lb. 


' 

ne*  or  n  at- 

j 

■'  V.  n  :•  ■■  I  ■  ■  9>  HC  JOT: 

b  yq.U  '  3r.  °  ’■  T  "  tct  ?r  *  i:-  ...  %i‘ 

.  mo-  J’o  tt> 


*.i...  *  ■  c/s  32  ,  . »: .  : 

o  nr  J  •  .  .9%.  J  :  •*  c.  a  b  y .1*  *-**•■  '  >nR  !  Jl 

.  4'  ::  i"  a  '  L3  15  *uo i  i*.  / 

:  ..  aoiif  O'  •?’ 


(i  £  -  TJ  •;  )  +  b(-  -  AS  -  ’I ’ :  ) 


4 


J-  \U '  3  o  V  c  v3:tw  me  j  io;Joj  j  fltoj  sx^vr.  jo 
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Propane  -  vapor  region 

P-V-T  data (6)  p  -  14.696  - >2057. 44  psia 

(237  points)  T  -  545.67  - >1096.76  °R 

(H  -  Ho)t  data (21)  P  -  200.00  - >1000.0  psia 

(11  points)  T  -  659.63  - >  859.63  °R 


Fitting 

Criteria 

1 

2 

3 

4 

A 

o 

=r 

.  25980775xl05 

.  884  714  87xl04 

. 77471258x10^ 

. 81509683xl04 

Bo 

= 

.14833089X101 

.12128731x10° 

.18420332x10° 

.17099770x10° 

Co 

— 

. 56938962xl010 

. 91613260xl010 

.  96361525xl010 

. 95383007xl010 

a 

= 

.40377788x10s 

. 11381319xl06 

.78955402xl05 

.68293160xl05 

b 

= 

.47479564X101 

. 10426864xl02 

.73152495X101 

. 66016528X101 

c 

.20117470X1011 

. 34855595X1011 

.28184636X1011 

.26289272X1011 

aa 

= 

. 11229799xl06 

. 16273998x10° 

.13971655x10° 

.13186063x10° 

Y 

= 

.5654  x  101 

.4510  x  101 

.4600  x  101 

.4800  x  101 

1 

2 

3 

4 

Max. 

%  dev.  in  Z 

ss 

4.4413 

3.8930 

4.1442 

3.4140 

Max.  abs 

.  dev.  in  Z 

= 

0.01876 

0.01387 

0.01409 

0.0200 

Ave. 

%  dev.  in  Z 

= 

0.7105 

1.0755 

0.9545 

0.9595 

Ave.  abs 

.  dev.  in  Z 

= 

0.0035 

0.0058 

0.0051 

0.0053 

Max.  %  dev. 

in 

<H-ho>t 

= 

12.1024 

14.8361 

17.1674 

17.1674 

Max.  abs.  dev. 

in 

(h-Vt 

= 

7.0811 

3.8905 

2.0000 

2.0000 

Ave.  %  dev. 

in 

(H-Ho>  T 

— 

5.1321 

6.2422 

6.6384 

6.6097 

Ave.  abs.  dev. 

in 

(H-Ho>T 

= 

1.6915 

1.5107 

1.4681 

1.4746 

:  «:  %a 

tvr:oj  j  >c.  >  -  on;  qoos  I 

siaq  M>.V20S-< - -  de&.*I  -  <3  (d)ei»b  T-V-<3 

a°  £d.C'c3 - r.Ccc  -  T  (eJtnioq  XI) 
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xa 
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oa 
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.  9VA 

ni  .vr*fc  . 

*soi.£i 

T? 

(  H-l)  fli 

0  )(.  . 

T 

,Vfc>  .ad*  .x*M 

S3 

T 

■  )  ni 

.  r<  3  v  .f./A 

{  H-H)  ni 
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1  -  Least  Squares 

2  -  Min.  max.  abs.  dev.  in  Z 

3  -  Min.  max.  abs.  dev.  in  Z,  subject  to; 

Max.  abs.  dev.  in  (H  -  Hq)t  4  2.0  BTU/lb. 

4  -  Min.  max.  %  dev.  in  Z  subject  to; 

Max.  abs.  dev.  in  (H  -  Hq)t  ^  2.0  BTU/lb. 

The  B-W-R  equation  is  difficult  to  fit  because  it  con¬ 
tains  a  non-linear  exponential  term.  This  necessitates  the 
use  of  a  trial . and .error  procedure  to  solve  for  the  optimum 
Y  value. 

To  obtain  optimum  y  values  for  the  various  fitting 
criteria,  a  trial  and  error  search  procedure  was  used.  A 
plot  of  fitting  error  versus  y  was  made  to  obtain  the  optimum 
Y.  A  final  solution  was  then  made  using  the  optimum  y  value. 

Fitting  criteria  2 ,  3  and  4  required  use  of  the  linear 
programming  method  of  data  fitting.  The  computer  programs 
required  for  this  method  are  given  in  Appendix  B  and  C,  along 
with  an  explanation  of  their  usage. 


;0^ 

.dI\UT£  O.S  t(oH 

Si  £  aii  oed  •  '  ci  i-'u  tiif>  ei  n-  i  -  :p  ■  - "  ed. 

,eso&a  ; ;  i ff T 

ftji 1  ?> .  frfd  i<  .  >r  o-  .'O'  i  icy  a<  l  .  ■>  '  xiP  l  io  ;*rj 

n  •  it  '.jo  •  ri:;  r-.c* .  o  ox  <  i  •'  r  uo  3  "ii'  ~  3 

I/--  :i  ..  -o'-  j  x  qt  tr  p'-r  i  *oi.  k  ,£  si  :i 
■  I.  TO;  .  •  O  X  ,  i ;;  i  -i'*  6-  io 

■  ;  £  )  .  ,  !b,:  >  .0  >  '"V  f-  I  •-  : 

.  opsax  i.  sniJ  .'O  x.  xn-  o  ns  tirrlw 
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Z  =  4>(P,T)  FORMS 

Equations  of  state  of  the  form  Z  =  <f>(P,T)  or  V  =  <j>(P,T) 
have  one  restriction.  They  may  only  be  applied  to  a  single 
phase  region,  since  the  function  would  otherwise  have  to  be 
dual  valued  at  pressures  and  temperatures  along  the  vapor  pres¬ 
sure  curve.  The  Z  =  <J>(P,T)  form,  however,  has  some  features 
which  make  it  attractive  for  engineering  application. 

In  almost  all  engineering  applications,  the  measured 
or  given  variables  are  pressure  and  temperature;  thus,  an 
equation  of  the  form  V  =  <f>(P,T)  can  be  solved  explicitly. 
Equations  of  the  form  P  =  <f>(V,T)  require  iterative  solutions 
for  volume.  In  terms  of  calculational  efficiency,  it  is  pre¬ 
ferable  to  have  a  fairly  complex  form  V  =  <j>(P,T)  ,  rather  than 
some  simpler  form  P  =  4>  ( V , T)  .  This,  of  course,  assumes  that 
the  same  degree  of  accuracy  is  obtainable  with  both  equation 
forms.  Furthermore,  the  derived  thermodynamic  properties  most 
useful  to  engineers  are  defined  as  functions  of  P  and  T.  As 
a  result,  equations  of  the  form  V  =  <J>(P,T)  give  simpler  ex¬ 
pressions  for  such  quantities  as  (H  -  Hq)t,  (S  -  Sq)t, 

(C  -  C  )  and  In  (f/P). 

P  Pq 

The  thermodynamic  quantities  can  be  expressed  as  func¬ 
tions  of  V  and  T,  but  the  functions  become  complicated  and  re¬ 
quire  an  intermediate  step  in  calculating.  The  given  variables 
are  generally  P  and  T,  requiring  an  iterative  solution  of  the 
equation  of  state  for  V  before  any  of  the  thermodynamic  pro- 


>.  n<  a*  tu$  jjeqiu&J  *>x  waj  -z  :  .  •  bt>  /.  sv  AJit:b 
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perties  can  be  calculated. 

Eight  equations  of  the  form  Z  »  <f>(P,T)  were  fitted  to 
propane  P-V-T  data  using  minimization  of  the  maximum  absolute 
deviation  in  Z  as  the  fitting  criterion.  The  equation  forms 
fitted  contained  from  six  to  twenty  terms.  However,  the  ac¬ 
curacy  of  fit  obtained  with  the  equations  tested  did  not  ap¬ 
proach  anywhere  near  the  accuracy  of  fit  obtained  using  the 
B-W-R  equation. 

The  investigation  indicated  that  equations  of  the  form 
Z  =  <f>(P,T)  may  require  forty  or  more  terms  in  order  to  approach 
the  accuracy  obtainable  with  the  B-W-R  equation.  The  investi¬ 
gation  was  carried  out  using  propane  P-V-T  data  limited  to 
pressures  of  two  thousand  psia.  To  make  an  equation  of  the 
form  Z  =  <J>(P,T)  applicable  to  pressures  of  ten  thousand  psia 
would  require  many  more  terms.  Thus,  it  was  concluded,  that 
equations  of  the  form  Z  =  4>  ( P , T )  would  be  far  too  complex. 

Any  advantage  gained  by  an  explicit  form  of  equation  would  be 
lost  due  to  the  highly  complex  nature  of  the  required  function. 

The  following  results  are  a  summary  of  the  various 
Z  =  <f>(P,T)  equations  fitted.  The  data  fitted  allows  direct 
comparison  of  fitting  accuracy  with  the  results  presented  for 
the  B-W-R  equation. 


otte  1/ ' 


b:  1  *3  ; 
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Analysis  of  Various  Forms  of  Z  =  $(p,T). 
Data : 

Propane  vapor(6),  237  points. 

P:  14.696 - >2057. 44  psia 

T:  545.67 - >-1096.76  °R 

1)  Equation:  (18  terms) 


Z  —  1.0  +  A,  —  +  A~  — -f  4-  A,.  — sr  +  A  „  — -r  +  Ac  — p-  +  Ar  — 7- 

1  T  2  3  T3  4  T4  5  T5  6  T6 


+  B,  +  B-  +  B  *  — -t  +  Bc  — =-  +  B>-  — 7* 

1  T  2  3  t3  4  t4  5  t5  6 


+  C,  —  +  C0  +  C-  — =■  +  C.  — t  +  Cc  — 7"  +  Q.  c  — 7" 
1  T  2  t2  3  t3  4  t4  5  t5  6  t6 


Maximum  absolute  error  in  Z  =  0D9624 

/ 

2)  Equation:  (20  terms) 


Z  =  1. 0  +  A,  —  +  A~  —  +  A-3  —  +  A.  — 

T 


+  Ar 


+  A, 


T 


+  A. 


p®  p  p2  p3  p4 

+  An  +  B,  — rr  +  B0  — rr  +  B — jr  +  B  -  — =* 

8  T  1  t2  2  t2  3  t2  4  t2 


P5  P^  p7  P8  P  P2 

+  Bc  - o'  +  - 7T  +  B_  - 7T  +  BQ  - 7T  +  C,  - ~  - ZT 

«p2  0  ip2  7  T2  “  t8  2  T8 


P  P 

+  -rr  +  C 


3  t3 


4  T3 


Maximum  absolute  error  in  Z 


0.06822 


' 

;  x  .  i  b  K  ,  v  c 

°  rat.)i - s 

XI  i  b i  • 


q 


‘q  'q 


q 


a  h 


S"r  +  +  + 


- 


:i  + 


.0  + 


T 


tsae.;-. 


"  '■  '  ;r 

I  0!  )  :  HO  '7T>  (S 


T 


8. 


t;  e  1 


A  + 


.A 


££850.0 
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3)  Equation:  (7  terms) 


Z  = 


1.0  +  A 

2 


1  -  +  A~  —  +  At  - T 

RT  R  RT 


P  P‘ 

+  A ,  — 


RT 


+  Ar  -  +  A,-  - t"  +  A-  - t" 

R  RT  RT 


Maximum  absolute  error  in  Z  =  0.1426 

4)  Equation:  (12  terms) 


P  P  P  P 

Z  —  1.0  +  A-.  —  +  A~  - 7T  +  A,.  - t-  "t  A.  - -r 

1  RT  A  RT  J  RTJ  RT 

2  2  2 
P  P  P  P 

+  Ar  g-  +  B.  +  B2  J  +  B-  g 

D  RT  RT  RT°  J  RT 


+  C,  —  +  C 

RT  “  RT 


^  - r  D,  - +  E,  — 

A  —  RT J  1  RT 


Maximum  absolute  error  in  Z  =  0.09965 

5)  Equation:  (11  terms) 


Z  —  1.0  +  A,  —  +  A~  -  +  At  - j  +  B,  — 

1  R  RT  RT  1  R 


+  B~  -  +  Bt  - t  +  C,  -  +  C~  - t" 

^  RT  J  RT  RT  RT 


+  D,  -  +  D_  - t"  +  E,  - 

1  R  ^  RTJ  RT 


Maximum  absolute  error  in  Z 


0.09434 


■ 

3  nol  3fpa 


0  -  nt  on*'  r  J:  ot  b  mamlxaM 
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6)  Equation:  (15  terms) 


P  P  P  P 

Z  =  1.0  +  A,  —  +  Ap  —  +  A~  sr  +  B,  — 

R  RT  J  RTJ  R 


+  —  +  B _  - 5-  +  C,  —  +  Cp  — 

RT  J  RT  R  RT 


+  C-  — — p-  +  D,  —  +  Dp  - 

RT  R  RT 


P  P 

-  +  D 


3  3 

J  RT J 


+  E. 


+  E. 


+  E 


R 


RT 


3  3 

RT 


Maximum  absolute  error  in  Z  =  0.09294 


7)  Equation:  (6  terms) 


Z  —  1.0+  A,  —  +  Ap  — 

1  R  RT 


+  B. 


+  B. 


RT 


+  C,  -  +  D,  - 

RT  R 


Maximum  absolute  error  in  Z  =  0.1125 


8)  Equation:  (8  terms) 


Z  =  1.0  +  A,  —  +  Ap  —  +  A 

1  RT  R 


3  3  +  ®i  +  B2 

RT  RT  R 


-y/p' 


+  C,  —  +  D,  P  - 

1  RT  1  RT 


3  +  D2 


P4  e-^P' 


RT' 


Optimum  y  =  10.0 


Maximum  absolute  error  in  Z 


0.1358 


T* 


'  3 


ji 
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Z  =  <f>(V,T)  FORMS 

In  the  search  for  a  better  equation  of  state  than  the 
B-W-R  equation,  the  equation  of  state  proposed  by  Martin  and 
Hou(15,16)  was  chosen  as  a  starting  point.  The  reasons  for 
this  choice  were  : 

1.  The  equation  is  easily  expanded  to  include  as  many  terms 
as  desired. 

2.  Extensive  physical  argument  is  given  for  the  form  of  the 
equation  by  Martin  and  Hou. 

3.  The  equation  is  calculationally  more  efficient  than  the 
B-W-R  equation.  Exponential  terms  in  the  equation  require 
evaluation  of  only  one  exponential  quantity  at  a  given 
temperature.  Thus,  the  exponential  quantity  is  excluded 
from  iterative  solutions  of  the  equation  for  volume. 

A  preliminary  investigation  of  the  Martin-Hou  equa¬ 
tion  in  its  original  form  was  undertaken  in  order  to  determine 
the  ability  of  the  equation  to  fit  the  P-V-T  data  for  propane. 
Results  of  this  investigation  are  as  follows: 

Preliminary  Investigation  of  the  Martin-Hou  Equation 

Data: 

Propane  vapor(6),  237  points. 

P:  14.696  - ►2057. 44  psia 


T:  545.67 


>1096.76  °R 


,  i>»j  i  .j  ; 
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1)  Original  Equation:  (8  terms) 


b  =  0.714598  (15) 


V 


Z  = 


A1V 


A2v 


A3V 


(V-b) 


RT(V-b) 


RT  (V-b) 


RT  (V-b) 


B1V 

+  - -  j  + 

R(V-b) 


B2V 


B4V 


R  (V-b) 


R  ( V-b) 


-5.475T 


CXV  e 


R 


-5.475T 


C2V  e 


R 


RT  (V-b) 


RT  (V-b) 


Maximum  absolute  error  in  Z  =  0.01524 

2)  Original  Equation:  suggested  exponential  term  added; 


b  =  0.714598  (15) 


V 


A±V 


a2v 


A3V 


(V-b) 


RT ( V-b) 


RT (V-b) 


RT (V-b) 


B1V 


B2V 


B  „V 
4 


R  ( V-b) 


R (V-b) 


RT  (V-b) 


-5.475T 


CXV  e 


R 


-5. 475T 


C2V  e 


R 


RT ( V-b) 


(V-b) 


-5.475T 


C4V  e 


R 


(V-b) 


Maximum  absolute  error  in  Z  =  0.01300 
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These  results  indicate  that  the  original  Martin-Hou 
equation  is  comparable  in  accuracy  of  fit  to  the  B-W-R  equation. 

The  original  Martin-Hou  equation  can  be  written  in 
the  form: 

-5.475T 

V  V  n  (A-  +  B.T  +  C.e  ) 

Z  =  - -  +  —  l  - i - i - ~ - 

(V-b)  RT  i=l  (V-b) 

This,  of  course,  is  an  expansion  in  l/(V-b)m.  In  order  to 
determine  the  best  expansion  form,  twelve  term  expansions  of 
the  equation  in  l/(V)ra,  l/(V-b)m  and  l/(V+b)m  were  made  with 
the  following  results. 


Data: 


Propane (6)  in  the  vapor  region 


-5.475T 


1) 


V  V  4 

z  =  -  +  —  l 

(V-b)  RT  i=l 


(A.  +  B.T  +  C.e 

i  l  l 


R, 


(V-b) 


T+T 


b  =  0.714598  (15) 

Maximum  absolute  error  in  Z  =  0.01212 

-5.475T 

1  4  (A.  +  B.T  +  C.e  ) 

2)  Z  =  1.0  +  —  l  - - - - - - 

RT  i=l  (V) 


Maximum  absolute  error  in  Z  =  0.01122 

-5.475T 

V  V  4  (A.  +  B.T  +  C.e  K) 

3)  Z  =  -  +  -  l  - 

(V+b)  RT  i=l  (V+b) x  x 


b  =  0.714598  (15) 


Maximum  absolute  error  in  Z 


0.01951 


V. 

7  1 


■  • : c  >.  V  0 
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The  results  indicate  that  for  the  vapor  region,  the 

i 

best  expansion  of  the  Martin-Hou  equation  is 

-5.475T 

1  n  (A.  +  B.T  +  C.e  ) 

Z  =  1.0  +  —  l  — i - ± - - 

RT  i=l  (V) 1 


Expansion  of  the  equation  in  l/V™  gives  the  best  re¬ 
sults  for  the  vapor  region.  However,  the  expansion  in  1/(V- 
ran  a  close  second  and  it  was  postulated  that  this  expansion 
might  have  merit  for  fitting  both  liquid  and  vapor  regions. 
In  order  to  check  out  this  possibility,  fifteen  constant  ex¬ 
pansions  of  the  equation  in  l/V^  and  l/(V-b)m  were  made  with 
the  following  results. 

Data : 


Propane(6),  in  the  vapor  and  liquid  regions. 

-5.475T 

1  5  (A.  +  B.T  +  C.e  K) 

1)  Z  =  1.0  +  —  l  — - - - - t— - - 

RT  i=l  (V) 


Maximum  absolute  error  in  Z  =  0.09227 


-5.475T 


2) 


V  V  5 

Z  =  -  +  —  1 

(V-b)  RT  i=l 


(A.  +  B.T  +  C.e 

l  1  1 


R, 


(V-b) 


ITT 


b  =  0.714598  (15) 


Maximum  absolute  error  in  Z  =  0.2949 

The  results  indicate  that  the  best  expansion  of  the 
Martin-Hou  equation  is,  regardless  of  the  region  to  be  fitted 


m 


-  S  3  n  •  srtF-q  • 


■* 
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of  the  form: 


-5.475T 

1  n  (A.  +  B.T  +  C.e  ) 

Z  =  1.0  +  —  j  — i - — - ri - 

RT  i=l  (V)1 


In  order  to  obtain  a  direct  comparison  of  fitting 
accuracy  for  fitting  the  gas  and  gas-liquid  regions,  a  fifteen 
constant  expansion  of  the  above  form  was  fitted  to  the  vapor 
region  only. 

Equation : 


-5.475T 

V  5  (A.  +  B.T  +  C.e  ) 

Z  =  1.0  +  —  l  — i - i - r± - 

RT  i=l  (V)1 


1)  Data: 

Propane (6),  in  the  vapor  region  including  data  along 
the  vapor  pressure  curve. 

Maximum  absolute  error  in  Z  =  0.009177 

2)  Data: 

Propane(6),  in  the  vapor  region,  no  vapor  pressure 
data  included. 

Maximum  absolute  error  in  Z  =  0.009170 

A  direct  comparison  of  fitting  accuracy  can  be  made 
between  these  results  and  those  shown  previously  for  the  fif¬ 
teen  constant  expansion  in  l/V™,  fitted  to  the  vapor  and 
liquid  regions.  Comparison  of  the  results  reveals  an  order 
of  magnitude  loss  of  accuracy  upon  inclusion  of  the  liquid 
region  P-V-T  data.  There  is,  however,  no  appreciable  loss  of 
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accuracy  when  vapor  pressure  data  are  included  in  the  equation 
fit. 

It  was  concluded,  that,  because  of  the  serious  loss  of 
accuracy  created  by  fitting  both  phases  with  the  same  equation, 
the  equation  should  be  limited  to  the  vapor  phase.  This  does 
not  seriously  limit  the  usefulness  of  the  equation,  as  it  can 
always  be  tied  to  an  equation  for  the  liquid  phase  by  making  the 
vapor  and  liquid  fugacities  match  along  the  vapor  pressure 
curve . 

The  fact  that  fifteen  terms  were  used  in  the  equation 
has  no  real  significance.  The  choice  was  made  for  computation¬ 
al  reasons.  In  formulating  an  equation  of  state,  one  must 
reach  a  compromise  between  computational  efficiency  and  fit¬ 
ting  accuracy. 

The  best  nine  and  twelve  constant  expansions  of  the 
equation  are  given  for  comparative  purposes.  These  expansions 
were  obtained  by  trying  various  nine  and  twelve  term  combina¬ 
tions  of  the  terms  in  the  fifteen  constant  expansion. 

It  is  important  to  remember  that  the  equation  coef¬ 
ficients  are  not  linearly  independent,  therefore,  they  will 
have  very  different  values,  depending  upon  the  number  of  terms 
used.  Thus,  to  obtain  constants  for  the  nine  and  twelve  con¬ 
stant  equations,  the  data  must  be  fitted  using  these  equations. 
Data : 

Propane (6)  in  the  vapor  region,  including  data  along 

the  vapor  pressure  curve. 
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1)  15  Constant  equation 

-5.475T 

1  5  (A.  +  B.T  +  C.e  ) 

Z  =  1 #  o  +  _  l  -  1 _ i _ ,_i _ 

RT  i=l  (V) 1 


Maximum  absolute  error  in  Z  =  0.009177 

2)  Best  9  constants 


Maximum  absolute  error  in  Z  = 

3)  Best  12  constants 


C 


5 


0.01126 


A^ ,  A2 ,  A^ ,  At-  ,  B^ ,  B2,  B^/  Clt  C2  /  C^/ 

Maximum  absolute  error  in  Z  =  0.01072 

An  equation  explicit  in  Z  is  generally  not  a  very  use¬ 
ful  form?  therefore,  the  equation  was  rearranged  so  as  to  be 
explicit  in  pressure.  Fifteen  and  twenty-one  term  expansions 
of  the  equation,  explicit  in  pressure,  were  obtained  and 
fitted  to  the  propane  data.  These  results  give  an  indication 
of  the  pressure  errors  of  the  fit  and  at  the  same  time  show 
that  higher  order  expansions  do  give  a  significant  improvement 
in  fitting  accuracy. 

Data : 

Propane (6),  in  the  vapor  region  including  vapor 
pressure  data. 
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1)  Equation:  (15  terms) 


-5.475T 


R 


P 


RT 


+  l 


5 


(A.  +  B.T  +  C.e 
l  l  l 


V 


i=l 


Maximum  absolute  error  in  P  =  29.10  psia 

2)  Equation:  (21  terms) 


-5.475T 


R 


P 


RT 


+  l 


7 


(A.  +  B.T  +  C.e 

3.  1  X 


v 


i=l 


Maximum  absolute  error  in  P  =  23.68  psia 

It  was  discovered  from  these  results,  that 
minimizing  the  maximum  absolute  deviation  in  pressure  is  not 
a  very  good  fitting  criterion.  Some  of  the  maximum  deviations 
in  pressure  occurred  at  points  on  the  vapor  pressure  curve, 
where  pressures  are  relatively  low,  causing  large  percentage 
errors  in  pressure  (10%)  in  this  region. 


Percentage  deviations  have  more  significance  than  ab¬ 


solute  deviations  when  the  function  being  fitted  has  a  very 
large  range  of  values.  The  pressure  range  of  greatest  inter¬ 
est  runs  from  about  ten  to  ten  thousand  psia;  hence,  minimiza¬ 
tion  of  the  maximum  absolute  percentage  deviation  in  pressure 
was  chosen  as  the  better  fitting  criterion.  This  fitting 
criterion  was  applied  to  the  fifteen  constant  expansion  of  the 
equation  with  the  following  results. 
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Propane  data(6),  for  the  vapor  region,  including 
the  vapor  pressure  data. 


P 


RT  A. 

—  + 

V  V 


+ 


+ 


B.  T  B-T  B.T  B„T  B-T 

1  2  3  4  5 

— +  — o—  +  +  —e—  +  ~4— 

V2  V3  V  V5  V6 


+ 


-5 . 475Tr 
e 


Maximum  absolute  percentage  deviation  in  P  =  2.0260% 

Up  to  this  point  in  the  investigation  of  equation  of 
state  forms,  the  entire  study  was  restricted  to  the  fitting 
of  propane  P-V-T  data.  Restricting  the  investigation  to  one 
vapor  was  very  useful  for  establishing  the  best  form  of 
equation,  the  best  fitting  criterion,  and  the  range  of  appli¬ 
cation.  Having  established  this  much  for  propane,  it  was 
necessary  to  establish  the  usefulness  of  the  equation  for  a 
wide  variety  of  vapors. 

To  establish  the  general  utility  of  the  equation  a  de¬ 
tailed  fit  of  eight  different  vapors  was  carried  out.  The 
eight  vapors  fitted  were  methane,  propane,  ethylene,  benzene, 
carbon  dioxide,  hydrogen  sulfide,  nitrogen  and  argon.  The 
results  of  these  equation  fits  are  as  follows: 
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Equation: 


P  = 


RT 


V 


V 


V 


a3 

7 


V 


A 


V 


B,T 

+  —  + 


b9t  b^t  b.t 

-4—  +  -  +  — p — 

3  „4  „5 


V 


V 


V' 


+ 


b5t 

v^" 


c,  c9 

[4  +  4  +  -r  + 

V  V  V 


+ 


V' 
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R 
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Fitting  criterion: 

Minimize  the  maximum  absolute  percentage  deviation  in 

pressure 

Range  of  application: 

T  >  Tc  ;  All  pressures  up  to  the  highest  pressure  of 
interest. 

T  <  T  ;  All  pressures  up  to  and  including  the  vapor 

V— » 

pressure . 

(1)  Material:  Methane 
Equation  Constants : 


T 

— 

34  3 . 19  °R 

c 

Ai 

— 

-0.67068028 

X 

104 

a2 

— 

-0.13433947 

X 

105 

A3 

— 

+0.26809945 

X 

105 

A4 

= 

-0.19713974 

X 

105 

A5 

= 

+0.48385398 

X 

104 

Bn 

= 

+0.64013752 

X 

101 

:  1  iro  .  (I) 
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B2  =  +0.23022082  x  102 

B3  =  -0.30632364  x  102 

B4  =  +0.18616297  x  102 

B5  =  -0.11969532  x  101 

Cx  =  -0.79486746  x  106 

C2  =  +0.36788516  x  107 

C3  =  -0.57315354  x  107 

C4  =  +0.35269793  x  107 

C5  =  -0.72916887  x  106 

Molecular  weight  =  16.043 

Units:  ^ 

R  =  10.73147  (psia) ( ft  /lb-mole) /°R 

P  =  psia 

T  =  °R 

V  =  ft3/lb-mole 

Range  of  Data: 

P  :  14.696  - >14696  psia(12) 

T  :  365.69  - >851.  69°R 

(169  points) 

Max.  %  error  in  P  =  0.7648% 

Max.  abs.  error  in  P  =  112.39  psia 

Ave.  abs.  %  error  in  P  =  0.3923% 


Ave.  abs.  error  in  P 


17.86  psia 
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(2)  Material:  Propane 
Equation  Constants: 


29 


Tc 

ss 

666. 0°R 

A1 

— 

+0.15065950 

X 

105 

a2 

ss 

-0.11126574 

X 

107 

A3 

s: 

+0.79445232 

X 

107 

a4 

ss 

-0.22610278 

X 

00 

o 

1 — 1 

A5 

— 

+0.21827919 

X 

108 

B1 

= 

-0.45698786 

X 

102 

B2 

— 

+0.14768928 

X 

io4 

B3 

ss 

-0.10247113 

X 

io5 

B4 

ss 

+0.28824875 

X 

105 

B5 

ss 

-0.27506201 

X 

105 

C1 

ss 

-0.27058287 

X 

107 

C2 

ss 

+0.31838750 

X 

108 

C3 

ss 

-0.21227985 

X 

io9 

C4 

ss 

+0.67206797 

X 

io9 

C5 

— 

-0.72138224 

X 

io9 

Molecular  Weight  = 

44.097 

Units : 

R  =  10.73147  (psia) ( ft 3/lb-mole) /°R 

P  =  psia 

T  =  °R 

=  ft'Vlb-mole 
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Range  of  Data: 


P  : 

14.696 

2057  psia(6 

T  : 

545.69 

1096. 78°R 

(245 

points) 

Max. 

%  error 

in 

P 

=  2.0263% 

Max.  abs 

.  error 

in 

P 

=  41.69  psi 

Ave.  abs. 

%  error 

in 

P 

=  1.0292% 

Ave.  abs 

.  error 

in 

P 

=  10.33  psi 

(3)  Material:  Ethylene 


Equation  Constants: 


Tc 

508. 79°R 

A1 

= 

-0.22753538 

X 

10 

a2 

— 

+0.79991381 

X 

10 

A3 

— 

-0.26390817 

X 

10 

A4 

= 

+0.28374683 

X 

10 

A5 

= 

-0.94597719 

X 

10 

B1 

= 

+0.26812945 

X 

10 

b2 

= 

-0.14424227 

X 

10 

B3 

= 

+0.47236436 

X 

10 

B4 

= 

-0.49954036 

X 

10 

B5 

— 

+0.17726410 

X 

10 

C1 

rs 

-0.59373207 

X 

10 

C2 

= 

-0.17246599 

X 

10 

C3 

= 

+0.13444364 

X 

10 

C4 

s 

-0.18112673 

X 

10 

Cc 

= 

+0.68660545 

X 

10 

5 

5 

6 
6 

5 
2 
3 
3 
3 
3 

6 

7 

8 
8 
7 


4~ 

T 

*  -.s’ 

V  i  ..  ■  •  // 

as  -  .  — 

x  :as:*o- 
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31 


Molecular  weight  =  28.054 

Units : 

R  =  10.73147  (psia) (ft3/lb-mole)/°R 

P  =  psia 
T  =  °R 

V  =  ft3/ lb -mole 
Range  of  Data: 


P  ;  14.878  - — >-37196  psia(7) 

T  :  491.69  - *761.69°R 

(231)  points) 


Max. 

% 

error  in  P 

ss 

2.3992% 

Max.  abs. 

error  in  P 

— 

892.44  psia 

Ave .  Abs . 

% 

error  in  P 

ss 

1.3188% 

Ave .  abs . 

error  in  P 

ss 

90.68  psia 

(4)  Material: 

Benzene 

Equation  Constants : 

T 

— 

1012. 79°R 

c 

Ai 

ss 

+0.81343905 

X 

105 

A2 

= 

-0.38590771 

X 

107 

A3 

ss 

+0.42759643 

X 

OO 

o 

i — i 

A4 

= 

-0.26575865 

X 

io9 

A5 

— 

+0.62507683 

X 

cr> 

o 

i — 1 

B1 

ss 

-0.11161890 

X 

ro 

O 

i- 1 

B2 

ss 

+0.35639773 

X 

o 

i — i 

B3 

ss 

-0.38597862 

X 

105 

ir  r  • ;  -J  e-lCil 


T 


eq  0  fji  IO^V.9  .  J  v€>vA 
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B4  =  +0.23327206  x  106 

B5  =  -0.53642336  x  106 

Cx  =  -0.47924205  x  107 

C2  =  +0.56489854  x  108 

C3  =  -0.64244265  x  109 

C4  =  +0.62202228  x  1010 

C5  =  -0.18566092  x  1011 

Molecular  weight  =  78.114 

Units: 

R  =  10.73147  (psia) ( ft 3/lb-mole) /°R 

P  =  psia 

T  =  °R 

V  =  ft8/lb-mole 

Range  of  Data: 


P  :  375.6  - ►930.4  psia(9) 

T  :  952.02  - ►1134. 32°R 


(187  points) 

Max.  %  error  in  P 
Max.  abs.  error  in  P 
Ave.  abs.  %  error  in  P 
Ave.  abs.  error  in  P 


0.3865% 
3.60  psia 
0.1865% 
1.27  psia 


(5)  Material:  Carbon  Dioxide 

Equation  Constants : 

T  =  54  7. 56  °R 

c 

A±  =  -0.26280181  x  105 


. 


< 

Cds> 

f  ~Sf.  •  :.{  * ■;  ■ 
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A2 

+0.96207443 
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X 

•  ^ 

105 

A3 

= 

-0.21501424 

X 

106 

A4 

s - 

+0.17530775 

X 

I06 

A5 

— 

-0.47329448 

»• 

X 

105 

BX 

ss 

+0.27269119 

X 

102 

B2 

ss 

-0.12725418 

X 

ro 

O 

i — 1 

B3 

-5 

+0.29627500 

X 

o 

U) 

B4 

ss 

-0.24136579 

X 

103 

B5 

= 

+0.66861749 

X 

102 

C1 

ss 

+0.92167559 

X 

105 

C2 

ss 

-0.49328657 

X 

107 

C3 

ss 

+0.14323348 

X 

108 

C4 

ss 

-0.13172478 

X 

108 

C5 

ss 

+0.38213556 

X 

107 

Molecular  weight  = 

44.011 

Units  : 

R  =  10.73147  (psia) (ft3/lb-mole)/°R 

P  =  psia 

T  =  °  R 

V  =  ft3/lb-mole 

Range  of  Data: 

P  :  7.35  - *  44088  psia(7) 

T  :  401.69  - *-761.69°R 


(242  points) 
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Max.  %  error  in  P  = 
Max.  abs.  error  in  P  = 
Ave.  abs.  %  error  in  P  = 
Ave.  abs.  error  in  P  = 


1.5957% 
703.52  psia 
0.5670% 
50.26  psia 


(6)  Material:  Hydrogen  Sulfide 


Equation  Constants: 


Tc 

672. 4°R 

A1 

= 

-0.17219137 

X 

105 

a2 

= 

-0.12356671 

X 

106 

A3 

= 

+0.53008304 

X 

106 

a4 

= 

-0.72407331 

X 

106 

a5 

ss 

+0.31145722 

X 

106 

B1 

= 

-0.24436326 

X 

101 

B2 

ss 

+0.24040259 

X 

103 

B3 

ss 

-0.83492352 

X 

103 

B4 

— 

+0.10516778 

X 

104 

B5 

ss 

-0.43005513 

X 

103 

C1 

ss 

+0.29752710 

X 

106 

C2 

ss 

-0.18326696 

X 

107 

C3 

ss 

-0.11910975 

X 

107 

C4 

= 

+0.91815119 

X 

107 

C5 

— 

-0.61695012 

X 

107 

Molecular  weight  : 

34.082 

Units : 


10.73147  (psia)  ( ft3/lb-mole) /°R 


R 
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P  =  psia 
T  =  °R 

V  =  ft'Vlb-mole 

Range  of  Data: 


P  :  14.696 

T  :  499.69 

(202  points) 


-*10,000  psia (20) 


-*799 . 69°R 


Max 

.  % 

error  in  P 

ss 

4.6444% 

Max.  abs. 

error  in  P 

= 

464.45  psia 

Ave.  abs 

.  % 

error  in  P 

= 

2.0817% 

Ave.  abs. 

error  in  P 

= 

44.62  psia 

(7)  Material: 

Nitrogen 

Equation  Constants : 

Tc 

= 

226. 8°R 

A1 

= 

-0.67355560 

X 

104 

a2 

= 

+0.75599856 

X 

104 

A3 

= 

-0.14046270 

X 

105 

A4 

55 

+0.12607564 

X 

105 

A5 

ss 

-0.40924410 

X 

104 

B1 

ss 

+0.11494687 

X 

102 

B2 

— 

-0.99447537 

X 

101 

B3 

ss 

+0.27365071 

X 

102 

B4 

ss 

-0.24403024 

X 

102 

Be 

ss 

+0.91193732 

X 

101 

...  .  ' 
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C±  =  -0.99546378  x  105 

C0  =  -0.21271695  x  106 

C3  =  +0.21250628  x  107 

C.  =  -0.41404779  x  107 

4 

Cr  =  +0.19459221  x  107 

5 

Molecular  weight  =  28.016 

Units : 

R  =  10.73147  (psia) ( ft3/lb-mole) /°R 

P  =  psia 

T  =  °  R 

V  =  ft3/lb-mole 

Range  of  Data: 


P  :  6.6  - *  14741  psia(17) 

T  :  221.69  ►  761.69°R 


(109  points) 

Max.  %  error  in  P 
Max.  abs.  error  in  P 
Ave.  abs.  %  error  in  P 
Ave.  abs.  error  in  P 


0.5223% 
76.53  psia 
0.2591% 
11.78  psia 


(8)  Material:  Argon 

Equation  Constants : 

T  =  271. 3°R 

c 

A1  =  -0.43563009  x  104 

A2  =  +0.50761453  x  103 

A3  =  -0.90829412  x  104 


-  df 

0,  x  8£c5(  5S1S.  .0+ 

01  X  .'.ssce^^x.o-: 

310 . bS  lii  iifc  1c  . 
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BiA  -tv  0 ,  ■  •' 


i Dfiiecj  ITU 

: 

8  dl>  0VA 

(8) 

♦«noD  nc  i;*fiwpS 

T 

?oi  £* 

si^es8oe.o- 

* 

-  37  - 

a4  = 

+0.11013453  x  105 

A5  = 

-0.35006883  x  104 

B1  = 

+0.79613487  x  101 

b2  = 

-0.85922136  x  101 

B3  = 

+0.34658412  x  102 

B4  = 

-0.32599490  x  102 

II 

in 

+0.98217947  x  101 

ci  = 

-0.15747792  x  106 

c2  = 

-0.44802253  x  106 

c3  = 

+0.21020384  x  107 

n 

II 

-0.20411990  x  107 

C5  = 

+0.58566869  x  106 

Molecular  weight  =  39.944 


Units : 

R  = 

10.73147  (psia) (ft3/lb-mole)/°R 

P  = 

psia 

T  = 

°R 

V  = 

ft3/lb-mole 

Range  of  Data: 


P  : 

14.696  - >14696  psia(7) 

T  : 

180.0  - >1080 °R 

Max.  % 

error  in  P  =  3.7363% 

Max.  abs. 

error  in  P  =  549.1  psia 

Ave .  abs .  % 

error  in  P  =  1.5214% 

Ave.  abs. 

error  in  P  =  47.23  psia 
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The  fitting  accuracy  of  the  fifteen  constant  equation, 
indicated  by  the  previous  results,  was  very  good  considering 
the  wide  range  of  pressures  and  temperature  covered.  It  is 
important  to  remember,  however,  that  the  proposed  equation  is 
empirical  and,  therefore,  the  accuracy  of  the  equation  is  only 
guaranteed  within  the  range  of  the  variables  fitted. 

There  is  one  further  modification  of  the  equation  which 

could  be  made  to  improve  the  P-V-T  fit  of  a  particular  vapor. 

The  exponential  term  in  the  equation  could  be  expressed  in 
~BTR 

the  form  e  ,  with  B  being  another  constant  to  be  determined 
for  a  particular  vapor.  This  seemingly  minor  modification, 
however,  makes  the  equation  non-linear,  requiring  an  iterative 
search  procedure  for  obtaining  the  optimum  B  value.  Since 
this  study  was  aimed  primarily  at  establishing  the  best  equa¬ 
tion  form,  the  approach  taken  by  Martin  and  Hou  was  adhered  to. 

Martin  and  Hou  suggested  that  one  value  of  B  (5.475) 
is  applicable  to  all  vapors  and  this  was  their  primary  reason 
for  using  the  exponential  temperature  term  in  their  equation. 
They  did  considerable  work  to  establish  the  best  value  of 
B(5.475);  hence,  this  value  was  used  in  all  work  done  on  the 
equation  in  this  thesis. 

The  proposed  equation  of  state  gives  an  improved  re¬ 
presentation  of  the  volumetric  behavior  of  pure  vapors,  but, 
it  suffers  from  a  deficiency  common  to  all  equations  of  the 
form  P  =  <J>(V,T).  Equations  of  this  form  can  predict  pressures 
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to  an  accuracy  predicted  by  the  fitting  procedure.  Calcula¬ 
tion  of  volume  values,  however,  is  an  entirely  different  story. 
The  accuracy  of  calculated  volumes  is  dependent  upon 
In  regions  where  l(•^^)vl  is  less  than  0.5,  calculated 

volumes  will  be  less  accurate  than  calculated  pressures,  where 

3P, 


(*5^r)v I  is  greater  than  0.5,  calculated  volumes  will  be  more 
accurate  than  calculated  pressures.  This  means  that  in  most 
of  the  vapor  region,  calculated  volumes  will  be  quite  accur¬ 
ate.  In  the  vicinity  of  the  critical  point,  however,  such  is 

3P 


not  the  case.  In  this  region  | (jrj^v'  is  approaching  zero  or 
very  small  values.  This  means  that  volumes  calculated  in  the 
region  of  the  critical  point  will  be  in  error  by  a  factor  of 
several  times  the  predicted  error  in  pressure.  For  this  reason, 
equations  of  state  which  are  explicit  in  pressure  are  not  re¬ 
commended  for  calculating  volumes  in  the  region  of  the  critical 
point. 


Having  established  the  proposed  equation's  ability  to 
predict  accurate  volumetric  relationships  in  all  but  the  critical 
region,  work  was  undertaken  to  establish  the  equation's  ability 
to  predict  accurate  derived  thermodynamic  properties. 
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GENERAL  CHARACTERISTICS  OF  APPROXIMATING  FUNCTIONS 


Much  information  about  the  accuracy  of  various  derived 
thermodynamic  properties  can  be  gained  by  looking  at  some 
general  characteristics  of  approximating  functions. 

Approximating  functions  are  exactly  what  the  name  im¬ 
plies.  They  are  an  approximation  to  some  exact  function. 

This  means  that  there  will  be  some  error  associated  with  each 
point  on  the  functional  curve  obtained  by  fitting  the  approxi¬ 
mating  function  to  the  data.  As  a  result  of  the  associated 
errors,  approximating  functions  are  cyclic  in  nature  in  that 
they  follow  the  true  function  but  cycle  above  and  below  it. 

As  the  number  of  terms  in  the  approximating  function  is  in¬ 
creased,  the  amplitude  of  the  deviations  decreases,  but  at  the 
same  time  the  frequency  of  cycling  increases. 

The  cyclic  nature  of  the  approximating  function  causes 
large  deviations  in  the  predicted  derivative  values.  Maximum 
deviations  in  the  derivatives  occur  where  the  approximating 
function  has  zero  deviation  (crosses  the  actual  function) . 
Integrals  of  the  approximating  function  will  usually  be  very 
accurate  due  to  the  partial  cancellation  of  positive  and  nega¬ 
tive  deviations. 

Derived  thermodynamic  properties  depend  upon  integration 
and  differentiation  of  an  approximating  function  (equation  of 
state) .  If  we  look  at  the  mathematical  processes  that  each 
derived  thermodynamic  property  is  dependent  upon,  we  can  obtain 
some  idea  of  the  relative  accuracy  of  the  calculated  values. 


>J  J  uli-.  :  ws  -  t  WC A  Ioi‘ 

' 


41 


RELATIVE  ACCURACY  OF  CALCULATED  THERMODYNAMIC  PROPERTIES 


i)  (aa)t 


f  P  dv 


(AF)t  =  (AA)t  +  A(PV)t 

f  (AF) 

In  (— )  =  - 

f  RT 

o 

The  accuracy  of  these  quantities  will  be  greater  than 
the  accuracy  of  calculated  pressures  since  all  of  the  quantities 
are  dependent  upon  the  process  of  integration. 

9P 

2)  (AE)  =  /  [T  ( — )  -  P]  dV 

3T 

(ah)t  =  (ae)t  +  a(pv)t 

3P 

(AS)  =  /  (— )  dV 

9T 

These  quantities  will  be  less  accurate  than  calculated 
pressures.  They  are  partially  dependent  upon  derivatives  of 
the  pressure  function.  Some  improvement  in  accuracy  is  attained 
by  integration  of  the  derivatives  but  not  enough  to  overcome 
the  inaccuracy  of  the  derivatives  themselves 


y 


3T 

(— ) 
9P 


H 
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9H 
(— ) 


l  av 
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Joule-Thompson  coefficients  will  be  quite  inaccurate 
because  of  their  direct  dependency  upon  derivatives. 


4) 


dzV 


( C  -  C  )  „  =  -T  f  ( - )  dP 


'p  '  T 
*0 


3T 


2  P 


Errors  in  the  functional  relationship  are  magnified  to 

a  very  great  extent  when  second  derivatives  of  the  function 

are  taken.  As  a  result,  calculated  values  of  C  at  elevated 

P 

pressure  will  be  unreliable. 
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IMPROVING  THE  PREDICTION  OF  THERMODYNAMIC  PROPERTIES 

One  of  the  primary  reasons  for  developing  equations 
of  state  is  their  usefulness  in  calculating  derived  thermo¬ 
dynamic  properties.  Having  an  equation  of  state  allows  the 
development  of  explicit  relationships  for  the  important  derived 
thermodynamic  properties.  The  most  important  derived  properties 
from  an  engineering  standpoint  are  (H  -  H  )T,  (S  -  Sq)t  and 
f/P.  Expressions  for  all  of  these  are  easily  derived  from  an 
equation  of  state. 

Most  of  the  derived  thermodynamic  property  values  re¬ 
ported  in  the  literature  have  been  obtained  by  graphical  inte¬ 
gration  and  differentiation  of  experimental  P-V-T  data.  Values 
obtained  in  this  manner  are  highly  dependent  upon  the  accuracy 
of  the  P-V-T  data  and  the  graphical  techniques  employed  in 
their  calculation.  In  order  to  exclude  any  dependency  upon 
calculational  techniques,  it  was  decided  that  directly  measured 
experimental  values  for  the  thermodynamic  properties  would  be 
used  for  determining  the  accuracy  of  values  calculated  from  an 
equation  of  state. 

Of  the  three  important  derived  properties,  (H  -  Hq)t 
is  the  only  one  for  which  there  are  any  experimentally  deter¬ 
mined  data  available.  As  a  result,  the  study  of  derived  thermo¬ 
dynamic  properties  was  restricted  to  enthalpy  departure. 

A  method  of  improving  the  ability  of  equations  of  state 
to  predict  thermodynamic  properties  was  developed  and  applied 
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specifically  to  the  improvement  of  (H  -  Hq)t  predictions.  The 
method  was  first  applied  to  the  B-W-R  equation  in  order  to 
evaluate  its  usefulness.  Results  of  this  application  were 
presented  in  the  previous  discussion  of  the  B-W-R  equation  and 
they  justify  the  further  application  of  the  method. 

The  method  used  for  improving  enthalpy  departure  pre¬ 
dictions  includes  enthalpy  departure  data  into  the  equation 
of  state  fitting  procedure  and  can  be  formulated  as  a  linear 
programming  problem  in  the  following  manner: 

Minimize  the  maximum  absolute  deviations 

X  and  B,  subject  to 

|  P  -  f (V,T) |  i  «  * 

|  (H  -  Hq)t  -  g(V,T)  |  i  ^  6 

i  refers  to  a  particular  data  point. 

The  set  of  constraint  equations  formulated  cannot  be 
solved  directly,  as  linear  programming  allows  the  optimization 
of  only  one  objective  function.  The  formulated  problem  con¬ 
tains  two  objective  functions,  X  and  6.  There  are  two  means 
of  overcoming  this  difficulty. 

1)  Let  8  =  wX,  where  is  some  arbitrary  weighting  factor. 

2)  Let  8  =  K  ,  where  K  is  some  arbitrary  constant. 

Caution  must  be  exercised  when  using  method  (2) .  If 
the  constant  K  specified  is  too  small,  there  may  not  be  a 
solution  to  the  linear  programming  problem.  In  actual  applica— 
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tion,  the  computer  program  given  in  Appendix  C  for  solution 
of  this  problem,  will  print  "unbounded  solution"  and  cease 
calculating. 

The  most  useful  formulation  of  the  linear  programming 
problem  for  fitting  equations,  of  state  was  found  to  be 

P  -  f ( V,T) 

V“ 

-  „<  x 

P  i 


| (H  -  ho)t  -  g (V,T) \L  <  K 

i  refers  to  a  particular  data  point. 

Objective:  minimize  X  (maximum  %  error  in  P). 

Specify  K,  the  maximum  allowable  error  in  (H  -  Hq)t. 

Computer  programs  required  for  the  solution  of  this 
linear  programming  problem  are  given  in  Appendix  C,  along 
with  an  explanation  of  their  usage. 

Improvement  of  the  enthalpy  departure  predictions  was 
carried  out  on  the  proposed  fifteen  constant  equation  of  state 
Experimental  enthalpy  departures  could  only  be  found  for  meth¬ 
ane,  propane  and  benzene;  thus,  the  fitting  procedure  was  only 
applied  to  these  vapors.  Results  of  the  fitting  procedure  are 
as  follows: 

Equation  of  State: 


P  = 


RT 
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B1T 

Jr 
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b9t  b~t 

-4-  +  -4- 


V' 


V 


b4t 

v5 


b5t 

v6 


V 


C2  c3  C4  C5 ,  “5 . 4  75Tr 


V* 


V 


V* 
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Enthalpy  Departure  Equation: 


(H  -  Hq)t  =  J  [2 


A, 


V 


B^T 

V 


3  A2 

4 

+  — 

3 

A3 

5 

A4 

6 

A5 

2  V2 

V3 

i  — 

4 

v4 

T 

5 

V6 

b2t 

b3t 

4. 

b4t 

4. 

b5t 

V2 

v3 

v4 

T 

v5 

C-,  (2  +  5 . 475T_.) 


+  ( 


V 

C3(4 

+  5.47  5Tr) 

v3 

C5(6 

+  5.475TJ 
- , - 2-) 

C2(3  +  5 . 475Tr) 

7 


C4  (5  +  5 . 475Tr) 


V 


T 


R-, 


V' 


J  -  units  conversion  from  units  of  PV  to  BTU/lb. 
Fitting  Criterion: 

Minimize  the  maximum  percentage  error  in  P  subject 

to  the  maximum  absolute  error  in  (H  -  H  )m  <  2.0  BTU/lb. 
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(1)  Material:  Methane 
Equation  Constants: 
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To 

: s 

34  3 . 19  °R 

A1 

= 

-0.10485025 

X 

105 

a2 

a 

+0.50240222 

X 

■'T 

o 

1 — 1 

a3 

a 

-0.48803061 

X 

104 

a4 

a 

+0.32307928 

X 

104 

A5 

a 

-0.10977803 

X 

104 

B1 

a 

+0.12119481 

X 

102 

B2 

a 

-0.45553829 

X 

101 

b3 

a 

+0.16332147 

X 

102 

B4 

a 

-0.15146344 

X 

102 

B5 

a 

+0.74635761 

X 

101 

C1 

a 

-0.21888554 

X 

106 

C2 

a 

+0.82097618 

X 

106 

C3 

a 

-0.77795798 

X 

106 

C4 

a 

-0.89290670 

X 

105 

C5 

a 

+0.21571796 

X 

106 

Molecular  weight  = 

16.043 

Units : 


R  = 

J  = 

P  = 


10.73147  (psia 
144.0/ (mol.  wt 

psia 
°  R 


)  ( f t3/lb-mole) /°R 

.  x  778.173)  (in/  (lb-mole)  (BTU) 

(ft  ) (lb) ( ft-lbf ) 

V  =  ft2/lb-mole 

(H_Ho) T  = 


T 


BTU/lb 
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:  I  . 

itfrc  Oi  21.0- 

I3 

ss 

r3oi 

a3 

5/J.  t.vi.s1  v  '  .^w  .  £c .oso.^i-  *  x. 


T 


-  48  - 

Range  of  Data: 

(a)  PVT  data  :  P  :  14.696  - > 


(169  points)  T  :  365.69  - > 


(b) 

(H 

[-Hq)t  data 

• 

P  :  600 

(20  points) 

T  :  399. 

Max 

.  % 

error  in  P 

— 

0.7672% 

Max.  abs. 

error  in  P 

— 

112 . 74  psia 

Ave.  abs.% 

error  in  P 

— 

0.3692% 

Ave.  abs. 

error  in  P 

= 

16.92  psia 

Max.  %  error 

in  (H-Ho)t 

= 

5.0085% 

Max 

.  abs.  error 

in  (H-Hq)t 

= 

2.00  BTU/lb 

Ave. 

,  abs.%  error 

in  (H-Hq)t 

= 

2.6900% 

Ave 

.  abs.  error 

in  (H-Hq)t 

= 

1.17  BTU/lb 

(2) 

Material : 

Propane 

Equation  Constants: 

T 

c 

=r 

666. 0°R 

Ai 

ss 

-0.25705848 

X 

105 

A2 

= 

-0.46766287 

X 

106 

A3 

= 

+0.41155780 

X 

107 

A4 

= 

-0.13127352 

X 

00 

o 

1 — 1 

A5 

= 

+0.13617211 

X 

108 

B1 

ss 

+0.63466624 

X 

101 

B2 

= 

+0.62921939 

X 

103 

B3 

= 

-0.52240027 

X 

104 

B4 

= 

+0.16491636 

X 

105 

Bc 

-0.16918478 

X 

105 

5 

14696  psia(12) 

851 . 6  9  °R 

- >-2000  psia  (11) 

- >489 . 69°R 


(S  i  .  i.aq  iQd’I  * - 


{:  ...  . 

I  .  ftN 
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C1 

ss 

-0.16189794 

X 

107 

C2 

ss 

+0.16956579 

X 

108 

C3 

as 

-0.11384937 

X 

109 

C4 

— 

+0.39697617 

X 

109 

C5 

= 

-0.46090494 

X 

109 

Molecular  weight 

44.097 

Units : 

R 

— 

10.73147 

(psia) ( f t2/lb-mole) /°R 

(in2) (lb-mole) (BTU) 

J 

ss 

144.0/  (mol .' 

wt . 

x  778.173)  - 

(ftz) (lb) (ft-lbf) 

P 

ss 

psia 

V  = 

f t2/lb-mole 

T 

ss 

OR 

(H-Ho>  T  = 

BTU/lb. 

Range  of  Data: 

(a) 

PVT  data 

P  :  14.696 

- ►  2057  psia(6) 

(245  points) 

T  :  545.69 

- *  1096. 78°R 

(b) 

(H- 

-Hq)t  data 

P  :  200  — 

- >*1000  psia(21) 

(11 

points) 

T  :  659.63 

— — >-85  9 . 6  3  °R 

Max, 

.  % 

error  in 

P 

= 

2.1379% 

Max.  abs. 

error  in 

P 

ss 

43.99  psia 

Ave . abs 

i.  % 

error  in 

P 

— 

1.0220% 

Ave .  abs . 

error  in 

P 

ss 

10.72  psia 

Max.  %  error 

in  (H-HQ) 

T 

3= 

16.2192% 

Max. abs.  error 

in  (H-HQ) 

T 

ss 

2.00  BTU/lb 

Ave. abs. %  error 

in  (H-H  ) 

T 

ss 

6.2454% 

Ave.  abs.  error 

in  (H-H  ) 

T 

ss 

1.36  BTU/lb 

“ 
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(3)  Material:  Benzene 
Equation  Constants: 


Tc 

ss 

1012. 79°R 

A1 

— 

-0.14282083 

X 

106 

a2 

3S 

+0.37980274 

X 

107 

a3 

= 

-0.73156804 

X 

108 

A4 

St 

+0.47849122 

X 

109 

A5 

S3 

-0.96045465 

X 

109 

B1 

S3 

+0.92470515 

X 

102 

B2 

S3 

-0.33672543 

X 

104 

B3 

ts 

+0.63698605 

X 

105 

B4 

— 

-0.41041176 

X 

106 

B5 

ts 

+0.81809824 

X 

106 

C1 

= 

-0.16911234 

X 

106 

C2 

= 

-0.11270206 

X 

109 

C3 

St 

+0.24748603 

X 

1010 

C4 

St 

-0.16527034 

X 

1011 

C5 

S3 

+0.33384717 

X 

1011 

Molecular  Weight  = 

78.114 

Units : 

R  =  10.73147  (psia) (ft3/lb-mole)/°R 

2 

i-tox  (in  )  (lb-mole)  (BTU) 
J  =  144 . 0/ (mol.wt.  x  778.173)  - - ^ - 

(ftJ) (lb) (ft-lbf) 

P  =  psia  V  =  ft3/lb-mole 

T  =  °R  (h"h0}T  =  BTU/lb 


&;i  j^a  '  [ 


'  i  ,  i 


8 


a 


a 

a 


, 


Range  of  Data: 


P  :  375.6 


(a)  PVT  data 
(187  points) 

(b)  (H-Hq)t  data 
(8  points) 

Max.  %  error  in  P 

Max.  abs.  error  in  P 

Ave.  abs.  %  error  in  P 

Ave.  abs.  error  in  P 

Max.  %  error  in  (H-H  )m 

Max. abs.  error  in  (H-Hq) T 

Ave.  abs.%  error  in  (H-H  )_ 

o  T 

Ave. abs.  error  in  (H-H  )m 

O  T 


930.4  psia(9) 

T  :  952.02  - ►1134. 32°R 

P  :  429.42  - ^929. 69  psia(8) 

T  :  1012.79- — ►  1215. 35°R 
0.5717% 

5.32  psia 
0.2814% 

1.90  psia 
17.0068% 

2.00  BTU/lb 
5.4083% 

1.50  BTU/lb. 


The  results  of  the  improved  enthalpy  departure  fitting 
procedure  have  little  meaning  without  some  means  of  comparing 
the  results  with  the  enthalpy  departure  predicted  by  the 
equation  of  state  when  fitted  with  just  P-V-T  data.  To  make 
a  comparison  of  the  results,  programs  were  set  up  to  calculate 
enthalpy  departures  using  equation  constants  from  the  two 
fitting  procedures. 

Two  methods  of  calculating  enthalpy  departures  were 
employed.  For  comparison  of  the  data  fitting  procedures  re¬ 
sults  obtained  by  specifying  V  and  T  are  best.  Results  were 
also  obtained  by  specifying  P  and  T.  These  results  have 
greater  practical  interest  as  P  and  T  are  the  variables  which 
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are  usually  specified.  Examples  of  the  computer  programs 
used  for  both  calculational  procedures  are  given  in  Appendix 

D. 

(1)  Material:  Methane 

Calculation  of  Enthalpy  Departures  (20  points) : 

(a)  Given  P  and  T, calculate  V  and  (H-Hq)t. 


Constraint : 

A  B 


Max.  % 

error  in 

V 

3.5571% 

3.9477 

Max. 

abs . 

error  in 

V 

= 

0.0713  f t3/lb-mole 

0.0533 

Ave.abs.  % 

error  in 

V 

— 

0.9038% 

0.7657 

Ave . 

abs . 

error  in 

V 

= 

0.0310  ft3/lb-mole 

0.0210 

Max.  % 

error 

in 

<H-V 

T 

= 

16.8890% 

5.7782 

Max.  abs. 

error 

in 

(H-H0) 

T 

= 

5.8841  BTU/lb 

1.9529 

Ave .  abs .  % 

error 

in 

(H-H0) 

T 

— 

4.5577% 

2.9461 

Ave.  abs. 

error 

in 

(H-Ho) 

T 

— 

2.2317  BTU/lb 

1.2688 

(b)  Given  V  and  T,  calculate  P  and  (H-Hq) 

Constraint : 


A  B 

Max.  %  error  in  P  =  3.3190%  3.6189 

Max.  abs.  error  in  P  =  54.5964  psia  54.5127 

Ave.abs.  %  error  in  P  =  0.8359%  0.7153 

Ave.  abs.  error  in  P  =  11.0164  psia  10.2640 
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Max  %  error 

in 

(h-ho)t  = 

17.2982% 

5.0085 

Max.  abs. error 

in 

<h-Vt  = 

7.6532  BTU/lb 

2.0000 

Ave.abs.%  error 

in 

(H-Ho>t  = 

4.7163% 

2.6900 

Ave.  abs.  error 

in 

= 

2.4741  BTU/lb 

1.1662 

Constraints : 

A  -  PVT  fit;  minimize  the  maximum  percentage  error 
in  P. 

B  -  PVT  fit;  minimize  the  maximum  percentage  error 
in  P, subject  to, the  maximum  absolute  error  in 
(H-Hq)t  £  2.0  BTU/lb. 

(2)  Material:  Propane 

Calculation  of  Enthalpy  Departures ( 11  points): 

(a)  Given  P  and  T, calculate  V  and  (H-Hq)t. 

Constraint : 

A  B 


Max.  % 

error  in 

V 

8.4277% 

8.0664 

Max.  abs. 

error  in 

V 

= 

0.5523  ft3/lb-mole 

0.5286 

Ave. abs.  % 

error  in 

V 

— 

1.7192% 

1.4979 

Ave .  abs . 

error  in 

V 

= 

0.1823  ft3/lb-mole 

0.1469 

Max.  %  error 

in 

(H-H0) 

T 

= 

48.5140% 

15.3019 

Max.  abs.  error 

in 

(H“Hq) 

T 

— 

8.4032  BTU/lb 

2.5586 

Ave. abs.  %  error 

in 

(H-Ho) 

T 

= 

19.7735% 

6.3195 

Ave.  abs.  error 

in 

(H-H0) 

T 

= 

4.4893  BTU/lb 

1.3565 

V 
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(b)  Given  V  and  T, calculate  P  and  (H-Hq)t. 

Constraint : 

A  B 


Max.  % 

error  in 

P 

ss 

2.2995% 

2.3991 

Max.  abs. 

error  in 

P 

as 

19.5796  psia 

13.4348 

Ave . abs .  % 

error  in 

P 

= 

0.9391% 

0.7440 

Ave.  abs. 

error  in 

P 

ss 

5.9783  psia 

4.5741 

Max.  %  error 

in 

(H-Ho) 

T 

= 

48.8229% 

16.2192 

Max.  abs.  error 

in 

(H-H0) 

T 

= 

8.0130  BTU/lb 

2.0000 

Ave. abs.  %  error 

in 

T 

= 

19.1293% 

6.2454 

Ave.  abs.  error 

in 

T 

= 

4.1812  BTU/lb 

1.3614 

Constraints : 

A  -  PVT  fit;  minimize  the  maximum  percentage  error 
in  P. 

B  -  PVT  fit;  minimize  the  maximum  percentage  error 
in  P, subject  to, the  maximum  absolute  error  in 
(H-Hq)t  $2.0  BTU/lb. 

(3)  Material:  Benzene 

Calculation  of  Enthalpy  Departures (8  points): 

(a)  Given  P  and  T, calculate  V  and  (H-Hq)t. 

Constraint : 


A 

B 

Max.  % 

error 

in 

V  = 

0.0024% 

4.4733 

Max.  abs. 

error 

in 

V  = 

0.0003  f t3/lb-mole 

0.7229 

Ave. abs.  % 

error 

in 

V  = 

0.000305% 

1.4762 

Ave.  abs. 

error 

in 

V  = 

0.000031  ft3/lb-mole 

0.1833 

r.i 

i  »x 

' 

ad*>  .  ^ 

i  .6  5 .  VA 

“SOI  IB 

io.3»  .adii  *evA 

ft  ''  '  '  *■-  ' 

•  •  •  •  -•  ’■  ■  ■  ‘  -  >  •  •• 
rij  w  '  n  ^ 

.drVTTft  .  T(rH«L) 

.-rrfttjn ‘  -Z  :  i'.'fslh?  I 

I  -*  :  '  ,  ...i.  ,j  . 


* 

**  ::h  .  - 

V 

.X 

55 


Max.  %  error 

in 

(h-Vt 

= 

51.9400% 

14.672 

Max.  abs.  error 

in 

(B-Vt 

— 

14.1720  BTU/lb 

7.0248 

Ave.abs.  %  error 

in 

(B-Vt 

SB 

28.2720% 

5.5060 

Ave.  abs.  error 

in 

(B-Vt 

SB 

8.6639  BTU/lb 

1.4762 

(b)  Given 

V  and  T, 

calculate  P  and  (H-Hq)t. 

Constraint 

• 

• 

A 

B 

Max.  % 

error  in  P 

— 

0.0019% 

22.4251 

Max.  abs. 

error  in  P 

= 

0.0181  psia 

208.4839 

Ave . abs .  % 

error  in  P 

SB 

0.0010% 

3.4630 

Ave .  abs . 

error  in  P 

SB 

0.0072  psia 

29.7753 

Max.  %  error 

in 

(B-Vt 

= 

51.9400% 

17.0068 

Max.  abs.  error 

in 

<B-Ho>t 

= 

14.1714  BTU/lb 

2.0000 

Ave.abs.  %  error 

in 

<B-Ho>t 

— 

28.2717% 

1.5043 

Ave.  abs.  error 

in 

(B-Vt 

= 

8.6638  BTU/lb 

5.4083 

Constraints : 

A  -  PVT  fit,  minimize  the  maximum  percentage  error 
in  P 

B  -  PVT  fit,  minimize  the  maximum  percentage  error 
in  P  subject  to  the  maximum  absolute  error  in 
(H-Hq)t  .<2.0  BTU/lb . 

The  enthalpy  departures,  calculated  from  the  proposed 
equation  of  state,  using  equation  constants  obtained  by  both 
fitting  procedures,  indicate  the  usefulness  of  including 


-  '.’2 


ni 

1  .  >  •  • 

i  j 

as 

■ti 

T .  T-  .  •  B,  <7  A 

1 011-0  .3!  >'  •  JVA 

ft  ©vxD 

&■  .  xr* 

•  q  X810.0 

«f  it  "  a^!> 

*1  /  .'  tcne 

.  t  r 

)(  f  - 

t‘.'V  i 

- 

-  •  ,B  ■  ' 

as 

10119  .sd©  .xsM 

;  (o‘:“ 

10*1  t£»  £  ,  ntdc.svA 

KX 

,t  If  £-  A 

:  eS/ii.E  :.  isnoD 

' 


56 


enthalpy  departure  data  in  the  fitting  procedure.  Very  sub¬ 
stantial  improvements  in  calculated  enthalpy  departures  can  be 
obtained  with  only  a  small  loss  in  accuracy  in  predicted 
volumetric  behavior. 

There  are  obvious  discrepancies  in  the  pressures  predicted 
by  the  equation  of  state  using  either  set  of  equation  constants. 
The  pressure  discrepancies  become  evident  when  the  calculated 
maximum  percentage  errors  in  pressure  are  compared  to  the  maxi¬ 
mum  errors  predicted  by  the  linear  programming  data  fitting 
procedures.  The  pressure  discrepancies,  however,  occur  only 
at  points  in  the  enthalpy  departure  data  sets. 

The  pressure  discrepancies  in  the  enthalpy  departure 
data  sets  are  a  result  of  errors  in  the  volume  values  used. 

The  available  experimental  enthalpy  departure  data  was  re¬ 
ported  at  fixed  pressures  and  temperatures  with  no  correspond¬ 
ing  volume  values  given.  Unfortunately,  the  expression  for 
enthalpy  departure,  obtained  from  an  equation  of  state  expli¬ 
cit  in  pressure,  is  a  function  of  volume  and  temperature. 

Volume  values,  used  with  the  enthalpy  departure  data, 
for  propane  and  methane  were  obtained  from  the  P-V-T  data  by 
Lagrangian  interpolation.  Benzene  volumes  used  were  calcu¬ 
lated  from  the  proposed  equation  of  state  fitted  to  the 
P-V-T  data.  Use  of  either  of  these  methods  for  obtaining  the 
volume  values  results  in  some  error  in  the  volumes. 

Small  inconsistencies  between  the  P-V-T  and  enthalpy 
departure  data  do  not  disrupt  the  fitting  procedure  to  any 
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great  extent.  This  fact  is  borne  out  in  the  results  obtained 
by  calculating  enthalpy  departures  from  the  proposed  equation 
of  state  at  given  pressures  and  temperatures.  The  enthalpy 
departures  calculated  in  this  manner,  with  equation  constants 
from  the  enthalpy  constrained  fitting  procedure,  retain  their 
improvement,  even  though  volume  errors  are  incorporated  into 
the  calculated  enthalpy  departures. 

The  enthalpy  data  used  for  benzene  was  found  to  contain 
serious  errors.  Evidence  of  the  errors  in  these  data  is  pre¬ 
sented  by  Yarborough  and  Edmister ( 21) .  Inclusion  of  benzene 
enthalpy  departure  data  into  the  fitting  procedure  completely 
disrupted  the  equation  fit.  The  enthalpy  departure  constrained 
fit  for  benzene  is  quite  useless  but  does  illustrate  the  need 
for  careful  examination  of  the  thermodynamic  data  to  be  in¬ 
cluded  into  the  fitting  procedure. 

To  obtain  completely  consistent  results,  when  includ¬ 
ing  thermodynamic  data  in  the  fitting  procedure,  the  P-V-T 
and  thermodynamic  data  used  must  be  a  complete  set.  This 
means  that  there  must  be  a  P-V-T  data  point  for  every  thermo¬ 
dynamic  data  point  used  in  the  fitting  procedure.  If  data 
of  this  type  are  used,  the  maximum  errors  in  the  calculated 
properties  will  not  exceed  the  maximum  errors  predicted  by 
the  fitting  procedure. 

The  accuracy  of  a  particular  property,  however,  is 
only  guaranteed  over  the  range  of  variables  for  which  data 
for  that  property  has  been  included  in  the  fitting  procedure. 
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CONCLUSIONS 

The  proposed  fifteen  constant  modified  version  of  the 
Martin-Hou  equation  of  state  is  a  very  effective  means  of 
representing  the  p-v-T  behavior  of  a  wide  variety  of  vapors. 

The  form  of  the  equation  is  such  that  it  allows  the  use  of 
more  or  less  than  fifteen  constants,  depending  on  the  range 
of  data  to  be  fitted  and  the  degree  of  accuracy  desired.  The 
fifteen  constant  form  was  chosen  as  a  compromise  between  ac¬ 
curacy  and  calculational  efficiency. 

The  proposed  equation  of  state  contains  high  powers  in 
volume  and,  therefore,  the  method  of  linear  programming  data 
fitting  should  be  used  for  fitting  the  equation.  "Least  squares" 
data  fitting  will  result  in  an  ill-conditioned  set  of  equations 
when  applied  to  this  equation. 

Pressure  predictions  of  the  proposed  equation  of  state 
are  quite  accurate  over  the  entire  vapor  region.  The  equation, 
however,  will  not  predict  accurate  volumes  in  the  critical 
region.  This  is  a  common  failing  of  all  equations  of  state 
which  are  explicit  in  pressure. 

Volumes  can  be  calculated,  by  iterative  means,  from 
the  proposed  equation  of  state  more  efficiently  than  from  the 
B-W-R  equation.  The  calculational  efficiency  of  the  proposed 
equation  is  due  to  the  fact  that  the  exponential  quantities  in 
the  equation  are  not  involved  in  iterative  procedures  required 
for  calculating  volumes. 
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Thermodynamic  properties  calculated  from  the  proposed 
equation  of  state  are  no  more  accurate  than  those  calculated 
from  lower  order  equations  of  state.  Improved  thermodynamic 
properties,  however,  may  be  obtained  by  including  accurate 
thermodynamic  properties  in  the  equation  fitting  procedure. 
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RECOMMENDATIONS 

Further  work  is  required  to  make  the  proposed  equation 
of  state  more  useful,  since  the  primary  objective  of  this  in¬ 
vestigation  was  to  establish  the  basic  form  for  an  equation 
of  state. 

The  equation  of  state  presented  is  only  applicable  to 
the  vapor  phase;  therefore,  a  complimentary  equation  of  state 
for  the  liquid  phase  must  be  developed  for  a  complete  descrip¬ 
tion  of  the  volumetric  behavior  of  pure  fluids.  An  equation 
of  state  for  the  liquid  phase  could  be  made  compatible  with 
the  vapor  phase  equation  by  use  of  the  two  constraint  linear 
programming  data  fitting  procedure.  Pressures  or  fugacities 
at  points  on  the  vapor  pressure  curve  could  be  used  as  the 
binding  constraint  for  tying  in  the  liquid  phase  equation. 

The  usefulness  of,  and  method  for,  including  thermo¬ 
dynamic  properties  in  the  data  fitting  procedure  has  been  in¬ 
dicated.  Application  of  this  data  fitting  procedure  is  recom¬ 
mended  whenever  the  equation  of  state  may  be  used  for  the 
calculation  of  thermodynamic  properties.  Equations  of  state 
do  not  predict  accurate  thermodynamic  properties  when  fitted 
exclusively  to  P-V-T  data.  This,  of  course,  depends  upon  the 
particular  thermodynamic  property. 

Development  of  some  form  of  mixing  rule  for  the  equa¬ 
tion  is  required  to  make  it  useful  for  predicting  the  volumetric 
behavior  of  mixtures. 
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APPENDIX  A 


PROPOSED  EQUATION  OF  STATE 


A- 2 


Use  of  this  equation  is  only  recommended  for  pure 
vapors  within  the  range  of  P-V-T  data  used  to  obtain  the  equa¬ 
tion  constants.  The  equation,  however,  is  not  recommended  for 
calculating  volumes  in  the  critical  region. 

The  linear  programming  method  of  data  fitting  should 
be  used  for  obtaining  the  equation  coefficients.  Minimizing 
of  the  maximum  percentage  deviation  in  pressure  is  the  best 
fitting  criterion  in  most  instances. 
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APPENDIX  B 


P-V-T  DATA  FITTING  BY  LINEAR  PROGRAMMING 


B-2 


Part  I  Program 


Approximating 

Function : 
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]  e 

Fitting  Criterion 

• 

Minimize  the  maximum  percentage  deviation  in  pressure. 
Input  Information: 

M  -  Number  of  data  points 
N  -  Number  of  equation  coefficients 
TC  -  Critical  temperature 
P  -  Pressure 
T  -  Temperature 
V  -  Volume 
R  -  Gas  Constant 
Output  Information: 

Objective  value  -  maximum  percentage  deviation  in 

pressure  . 

Basic  Variables  -  indexes  of  the  controlling  data 

points.  Controlling  data  points 
being  those  where  the  maximum 
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deviations  occur, of  which  there 
will  be  N+l . 

For  indexes  greater  than  M,  attach 
a  minus  sign  and  add  M.  These  in¬ 
dexes  and  their  signs  are  required 
for  Part  II  program.  Data  points 
with  indexes  greater  than  M  are 
points  at  which  P-equation  is 
greater  than  P-experimental . 


Part  II  Program 


This  program  is  basically  a  Gaussian-Elimination  solu¬ 
tion  of  a  set  of  linear  equations.  The  linear  equations  to  be 
solved  are  the  approximating  function  equations  at  the  con¬ 
trolling  data  points. 

Input  Information: 

N  -  number  of  coefficients  plus  one 


R 

TC 

*  *p 

* 

**V 


gas  constant 

critical  temperature 

pressure 

temperature 

volume 


**SG  -  sign  of  the  index  (+1.0  or  -1.0) 

**  These  values  are  for  the  controlling  data  points  given  by 


the  Part  I  program. 
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Input  Information  for  Checking: 

MM  -  total  number  of  data  points. 

(same  as  M  for  Part  I) 

P / T , V  -  set  of  data  used  in  the  Part  I  program. 

Output  Information: 

1.  The  first  value  in  the  output  solution  is  the 
maximum  pressure  deviation.  This  value  should  check  with  that 
obtained  from  the  Part  I  program,  if  the  solution  for  the 
equation  coefficients  is  accurate. 

2.  The  remaining  values  in  the  solution  list  are  the 
equation  coefficients  in  the  order  which  they  appear  in  the 
equation . 

3.  A  check  of  the  equation  fit  is  printed  out.  This 
listing  provides  a  further  check  on  the  solution  for  the 
equation  coefficients.  The  maximum  deviations  should  occur 

at  the  controlling  data  points  if  the  equation  fit  is  correct. 
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P-V-T  FITTING  PROGRAM  -  PART  I 


$IBFTC  MAIN  NODECK, NOLIST 


*  * 

COMMON  A(22 , 550) ,B (22 

•k  k 

1  M,MM, N ,NA,BB 

*  * 

READ (5,1)  M,N 

*  * 

FORMAT (2X, 214) 

kk 

READ (5,300)  TC 

**300 

FORMAT (1H  , F15 . 8) 

*  * 

R=10. 73147 

** 

NA=N+1 

*  * 

MM=M+M 

** 

DO  102  1=1, M 

*  * 

READ (5,101)  P ,T , V, 

**101 

FORMAT ( 1H  , 3F15.8) 

*  * 

WRITE (6,400)  I , P ,T , V 

**400 

FORMAT ( 1H  , 13 , 3F15 . 8) 

** 

AT=1. 0/(V*V) 

*  * 

AU=AT/V 

•k  k 

AV=AU/V 

k  k 

AW=AV/V 

** 

AX=AW/V 

** 

EA=EXP (-5 . 475  *T/TC) 

** 

F ( I ) = (P-R*T/V) /P 

*  * 

A ( 2,1) =AT/P 

** 

A  (  3 , 1 ) =AU/P 

*  * 

A  (  4 , 1 ) =AV/P 

kk 

A ( 5 , I ) = AW/P 

** 

A ( 6 , I ) =AX/P 

** 

A ( 7 , I ) =AT*T/P 

*  k 

A ( 8 , I ) =AU*T/P 

k  k 

A ( 9 , I ) =AV*T/P 

k  k 

A (10 ,1) =AW*T/P 

k  k 

A ( 11 , I ) =AX*T/P 

kk 

A  ( 1 2 , 1 )  =  AT  *E  A/P 

k  k 

A  ( 13 , 1 ) =AU*EA/P 

k  k 

A  ( 14 , 1 ) =AV*EA/P 

k  k 

A  ( 1 5 , 1 )  =  AW  *E A/P 

**102 

A (16 ,1) =AX*EA/P 

*  ★ 

DO  5  J=1 ,M 

*  * 

JM= J+M 

*  k 

F ( JM) =-F ( J) 

k  k 

A ( 1 , J) =0 . 01 

kk 

A ( 1 , JM) =0.01 

DO  5  1=2, NA 
5  A (I , JM) =-A(I , J) 

BB=0. 00 
K(l)=l 
B  ( 1)  =1 
DO  7  1=2, NA 
XX= A (1,1) 

B (I) =-XX 
A (I , 1) =0 . 00 


, C ( 5  50 )  , F ( 55  0 )  ,K (22 )  ,RX(22)  , 
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DO  7  J=2,M M 

7  A(I, J)=A(I, J)-XX*A(1,J) 

DO  8  1=2, NA 

IF  (B  (I) . GE .0.0)  GO  TO  8 
B  ( I ) =-B ( I ) 

DO  9  J=1  ,MM 

9  A(I , J) =-A(I , J) 

8  CONTINUE 

DO  10  L=2 ,NA 
XX=A (L , 1) 

JA=1 

DO  11  J=2,MM 

IF  (A(L, J) .LE.XX)  GO  TO  11 
XX=A(L, J) 

JA=  J 

11  CONTINUE 
K(L)=JA 
B(L)=B(L)/XX 
DO  12  J=1 ,MM 

12  A(L , J) =A(L, J) /XX 
DO  13  1=1 ,NA 

IF  (I.EQ.L)  GO  TO  13 
YY=A ( I , JA) 

B ( I ) =B ( I ) -B (L) *YY 
DO  14  J=1 ,MM 

14  A ( I , J ) =A ( I , J ) -A ( L , J ) *YY 

13  CONTINUE 

10  CONTINUE 

40  DO  30  1=1, NA 

I F  ( B  ( I )  .  GE .  0 . 0 )  GO  TO  30 
I  A=I 

GO  TO  50 
30  CONTINUE 
GO  TO  31 

50  B  ( I  A)  =-B  ( IA) 

DO  51  J=1 , MM 

51  A(IA, J) =-A(IA, J) 

XX=0. 00 

JA=0 

DO  52  J=1 , MM 
DO  53  1=1, NA 
IF  ( J.EQoK (I) )  GO  TO  52 
53  CONTINUE 

IF  (A(IA,J) .LE.XX)  GO  TO  52 
XX=A(IA, J) 

JA=  J 

52  CONTINUE 
K ( I A) = JA 

B ( I A) =B ( I A) /XX 
DO  54  J=1 ,MM 
A(IA, J) =A(IA, J) /XX 
DO  55  1=1, NA 
IF  (I.EQ.IA)  GO  TO  55 
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56 

YY=A(I,JA) 

B ( I ) =B ( I ) -B ( I A) *YY 

DO  56  J=1,MM 

A  ( I ,  J) =A ( I  ,  J) -A ( IA, J) *YY 

55 

CONTINUE 

31 

GO  TO  40 

DO  20  1=1, NA 

20 

IA=K ( I ) 

BB=BB+F(IA) *B(I) 

22 

DO  21  J=1,MM 

C(J)=0.00 

DO  22  1=1, NA 

IA=K ( I ) 

C ( J) =C ( J) +F ( I A) *A ( I , J) 

21 

C  ( J)  =C  ( J)  -F  ( J) 

25 

WRITE  (6,25) 

FORMAT  ( 1HO ,  15X,  22H  FOR  INITIAL  TABLEAU 

-  ) 

26 

WRITE  (6,26)  BB 

FORMAT  ( 1HO ,  16X,  2 1HOB JECTI VE  FUNCTION  = 

,  E15.8) 

27 

WRITE  (6,27)  (K (L) ,  B (L) ,  L=1,NA) 

FORMAT  (1H  ,  18X ,  16HINDEP .  VECTOR  =  ,  14 , 

5X ,  E15 

67 

CALL  SIMPLEX 

CONTINUE 

STOP 

END 
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$IBFTC 


399 


300 


302 

301 


303 
307 

304 


306 

305 

310 

311 


DUAL  NODECK, NOLIST 

SUBROUTINE  SIMPEX 

COMMON  A ( 2 2 , 550 )  ,B(22)  ,C(550)  ,F(550)  ,K(22)  ,RX(22)  , 
1  M,MM,N ,NA,BB 
QX=C  (1) 

JA=1 

DO  300  J=2 , MM 
IF  (C(J).GE.QX)  GO  TO  300 
JA=  J 
QX=C ( J) 

CONTINUE 

IF  (QX.GE.0.0)  GO  TO  330 
DO  301  1=1, NA 

IF  ( A ( I , JA) . LE .0.0)  GO  TO  302 
RX ( I ) =B ( I ) /A ( I , JA) 

GO  TO  301 
RX  ( I )  =  - 1 . 0  0 
CONTINUE 
DO  303  1=1, NA 

IF  (RX(I) oLT.0.0)  GO  TO  303 
QA=RX ( I ) 

II=I 

GO  TO  304 
CONTINUE 
WRITE  (6,307) 

FORMAT  ( 1HO ,  15X,  20H  UNBOUNDED  SOLUTION.) 

GO  TO  398 

IF  (II.GE.NA)  GO  TO  310 
IC=II+1 

DO  305  I=IC ,NA 

IF  (RX(I) .LT.O.O.OR.RX(I) .GT.QA)  GO  TO  305 
IF  (RX(I) .EQ.QA)  GO  TO  306 
QA=RX (I) 

II=I 

GO  TO  305 

IF  (K  ( II )  .  GE.K(I))  GO  TO  305 
QA=RX  ( I ) 

II=I 

CONTINUE 
DX=A(II, JA) 

B ( I I ) =B (II) /DX 
DO  311  J=1 , MM 
A ( II , J)=A(II , J) /DX 
DO  312  1=1, NA 
IF  (I.EQ.II)  GO  TO  312 
DXX=A(I, JA) 

B(I)=B(I) -DXX*B(II) 

DO  313  J=1 , MM 
IF  (J.EQ.JA)  GO  TO  314 
A(I , J) =A(I , J) -DXX*A (II , J) 

GO  TO  313 
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314  A ( I , J) =0 . 00 

313  CONTINUE 

312  CONTINUE 

BB=BB-C ( JA) *B ( II ) 

DO  315  J=1,MM 
IF  (J.EQ.JA)  GO  TO  315 
C(J)=C(J) -C ( JA) *A (II , J) 

315  CONTINUE 

C ( JA) =0.00 
K ( II ) =JA 
GO  TO  399 

330  WRITE  (6,331) 

331  FORMAT  (1HO,  30H  THIS  IS  THE  OPTIMUM  SOLUTION.  ) 
WRITE  (6,321)  BB 

321  FORMAT  (1H  ,  20X,  22H  OBJECTIVE  FUNCTION  =  ,  E15.8) 
WRITE  (6,322) 

322  FORMAT  (1HO,  25X,  4  8HBASIC  VARIABLES  AND  VALUES 

1ARE  AS  FOLLOWS) 

WRITE  (6,323)  (K(I),  B(I),  1=1, NA) 

323  FORMAT  (1H  ,  29X,  14,  12X,  E15.8) 

398  CONTINUE 

RETURN 

END 
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P-V-T  FITTING  PROGRAM  -  PART  I 
SAMPLE  OUTPUT 


P  (psia) 


T  ( °  R) 


PROGRAM  IS  BEING  ENTERED  INTO  STORAGE 


99.15390848 

120.85990144 

145.90188544 

174.58847744 

207.25768448 

244.21812224 

285.82249728 

332.37942784 

384.25630208 

441.83523840 

505.49830144 

575.67170560 

652.86979584 

737.62162688 

830.69138944 

933.00494336 

1070.60359168 

7.34799992 

14.69599984 

29.39199968 

44.08799872 

58.78399936 

73.47999936 

88.17599744 

7.34799992 

14.69599984 

29.39199968 

44.08799872 

58.78399936 

73.47999936 

88.17599744 

117.56799872 

146.95999872 

176.35199488 

7.34799992 

14.69599984 

29.39199968 

44.08799872 

58.78399936 

73.47999936 

88.17599744 

117.56799872 


401.68798720 

410.68799488 

419.68798720 

428.68799488 

437.68798720 

446.68799488 

455.68798720 

464.68799488 

473.68798720 

482.68799488 

491.68798720 

500.68799488 

509.68798720 

518.68799488 
527.68799232 

536.68799488 
547.55998720 

401. 68798720 

401.68798720 

401.68798720 

401.68798720 

401.68798720 

401.68798720 

401.68798720 

437.68798720 

437.68798720 

437.68798720 

437.68798720 

437.68798720 

437.68798720 

437.68798720 

437.68798720 

437.68798720 

437.68798720 

473.68798720 

473.68798720 

473.68798720 

473.68798720 

473.68798720 

473.68798720 

473.68798720 

473.68798720 


V ( f t^/mole) 


39.06036032- 

32.29404288 

26.90453504 

22.56472992 

19.03493760 

16.13327472 

13.72297680 

11.70887280 

10.00636672 

8.55905992 

7.31971984+ 

6.14622424 

5.19888340 

4.38689720 

3.63574996 

2.91308352 

1.50666524- 

582.54621696 

289.17934336 

142.50295552 

93.56382720 

69.05901504 

54.34625856 

44.50485504 

635.96191232 

316.35599872 

156.55303936 

103.22193792 

76.58811264 

60.59230656 

49.89788864 

36.48227872 

28.48085184 

23.08781888 

689.08856320 

343.20836608 

170.25768832 

112.57691264 

83.75062464 

66.44357184 

54.88908672 

40.40897120 
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P  (psia) 


T  (  °R) 


V(ft'Vmole) 


146. 95999872 

176.35199488 

220.43999488 

293.91999744 

367.39999744 

7.34799992 

14.69599984 

29.39199968 
44.08799872 

58.78399936 

73. 47999936 

88.17599744 

117.56799872 

146.95999872 

176.35199488 

220.43999488 

293.91999744 

367. 39999744 
440.87998976 
587.83999488 

7.34799992 

14.69599984 

29.39199968 
44.08799872 

58.78399936 

73.47999936 

88.17599744 

117.56799872 

146.95999872 

176.35199488 

220.43999488 

293.91999744 

367.39999744 
440.87998976 
587.83999488 
734.79999488 
881.75997952 

7.34799992 

14.69599984 

29.39199968 
44.08799872 

58.78399936 

73.47999936 

88.17599744 

117.56799872 

146.95999872 

176.35199488 

220.43999488 


473.68798720 

473.68798720 

473.68798720 

473.68798720 

473.68798720 

509.68798720 

509.68798720 

509.68798720 

509.68798720 

509.68798720 

509.68798720 

509.68798720 

509.68798720 

509.68798720 

509.68798720 

509.68798720 

509.68798720 

509.68798720 

509.68798720 

509.68798720 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

545.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 


31.76601952 
25.93590400 
20.12693872 
14.26157520 
10.61687192 
742,09538048 
369.93382912 
183.82838784 
121.78384000 
90.76509184 
72.14679232 
59.72519296 
44.18761856 
34.88199360 
28.63594528 
22.39694688 
16.12974992 
12.32995280 
9.76385640 
6.41383152 
795.01758464 
396.56061440 
197.31449344 
130.89912320 
97. 69496064 
77.76541568 
64.46965120 
47.85347072 
37.88517280 
31.21614144 
24.55416032 
17.87807936 
13.84564192 
11.13855072 
7.70533936 
5.57632512 
3.99154912 
821.44697344 
409.84227328 
204.03287296 
135.42504064 
101.12817280 
80.54300288 
66.81015744 
49.65114816 
39.35151328 
32.47099104 
25.59751776 
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P  (psia)  T  ( °R)  V ( f t 3/mole) 


293.91999744 
367. 39999744 

440. 87998976 

587.83999488 

734.79999488 

881.75997952 

1175.67998976 

1469.59998976 
2204.39998464 
2939.19997952 
4408.79996928 
5878.39995904 
7347.99986688 
8817.59993856 
11756.79991808 
14695.99973376 
17635.19987712 
22043.99927296 
29391.99946752 
14.69599984 
73.47999936 
146.95999872 

293.91999744 

440.87998976 

587.83999488 

734.79999488 

881.75997952 

1175.67998976 

1469.59998976 
2204.39998464 
2939.19997952 
4408.79996928 
5878.39995904 
7347.99986688 
8817.59993856 
11756.79991808 
14695.99973376 
17635.19987712 
22043.99927296 
29391.99946752 
14.69599984 
73.47999936 
146.95999872 

293.91999744 

440.87998976 

587. 83999488 

734.79999488 

881.75997952 

1175.67998976 

1469.59998976 


563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

563.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 
581.69899232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 


18.70994560 

14.55766320 

11.77302528 

8.26226696 

6.11280864 

4.60910420 

2.42510226- 

1.09623080 

0.90024870 

0.83891618 

0.76982896 

0.73528535+ 

0.70849642 

0.68699481 

0.65964191 

0.63792878 

0.62122096 

0.60112927 

0.57617328- 

423.10984704 

83.29944192 

40.79670560 

19.51361344 

12.37930080 

8.78394576 

6.58654952 

5.10681412 

3.12513900 

1.74550962 

1.00669941+ 

0.89954373+ 

0.80507755 

0.76207427 

0.72964558 

0.70454858- 

0.67395283 

0.65026579 

0.63186603 

0.61001191 

0.58315250 

449.60971776 

88.77002112 

43.65184064 

21.06455104 

13.52840480 

9.74975704 

7.45154968 

5.92528616 

3.97533476 

2.78040788 
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P  (psia) 


2204.39998464 

2939.19997952 

4408. 79996928 

5878.39995904 

7347.99986688 

8817.59993856 
11756.79991808 

14695.99973376 
17635.19987712 
22043.99927296 

29391.99946752 
36739.99933440 

14.69599984 
73. 47999936 
146.95999872 
293.91999744 
440.87998976 
587. 83999488 
734.79999488 
881.75997952 
1175.67998976 
1469.59998976 

2204.39998464 

2939.19997952 

4408.79996928 

5878.39995904 

7347.99986688 

8817.59993856 
11756.79991808 

14695.99973376 
17635.19987712 
22043.99927296 

29391.99946752 
36739.99933440 

14.69599984 
73.47999936 
146.95999872 
293.91999744 
440.87998976 
587.83999488 
734.79999488 
881.75997952 
1175.67998976 
1469.59998976 

2204.39998464 

2939.19997952 

4408.79996928 

5878. 39995904 

7347.99986688 

8817.59993856 


T  (°R) 


617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

617.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

653.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 


V ( f t ^/mole ) 


1.38033436 

1.06309713 

0.89037910 

0.82058691 

0.77617371 

0.74444998 

0.70342063 

0.67550376 

0,65350866 

0.62812968+ 

0.59746341+ 

0.57469284 

476,06731264 

94.20535168 

46.47172672 

22.58024000 

14.62111040 

10.63802120 

8.23406800 

6.62955272 

4.61333404 

3.41065248 

1.85901004 

1.29785269 

0.99682983 

0.88755922 

0.82834161 

0.78815823 

0.73528535 

0,70158772 

0.67550376 

0.64659993 

0.60860197- 

0.58745282 

502.50374144 

99.61248384 

49.25636480 

24.06772992 

15.67151808 

11.47693728 

8.95313904 

7.27530680 

5.16814660 

3.92105188 

2.27987810 

1.58477610 

1.12231475 

0.96510609 

0.88614927 

0.83398137 
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192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 

211 

212 
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218 
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229 

230 
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233 

234 

235 
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P  (psia) 


T  ( °  R) 


V ( f t3/mole) 


11756.79991808 

14695.99973376 

17635.19987712 
22043.99927296 

29391.99946752 

36739.99933440 

44087.99854592 
14.69599984 
73.47999936 
146.95999872 
293.91999744 
440.87998976 
587.83999488 
734.79999488 
881.75997952 
1175.67998976 
1469.59998976 
2204.39998464 
2939.19997952 
4408.79996928 
5878.39995904 
7347.99986688 
8817.59993856 

11756.79991808 

14695.99973376 

17635.19987712 
22043.99927296 

29391.99946752 

36739.99933440 

44087.99854592 
14.69599984 
73.47999936 
146.95999872 
293.91999744 
440.87998976 
587. 83999488 
734.79999488 
881.75997952 
1175.67998976 
1469.59998976 
2204.39998464 
2939.19997952 
4408.79996928 
5878.39995904 
7347.99986688 
8817.59993856 

11756.79991808 

14695.99973376 

17635.19987712 

29391.99946752 

36739.99933440 

44087.99854592 


689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

689.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

725.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 


0.76982896 

0.72823564 

0.69813335 

0.66542268 

0.62714272 

0.60042430 

0.58026211+ 

528.93312512 

105.01256448 

52.02690304 

25.52702112 

16.70077680 

12.28765472 

9.63696152 

7.88863200 

5.67502112 

4.36306904 

2.64505338 

1.86958462 

1.26401405 

1.05393250 

0.94959671 

0.88473933 

0.80578251 

0.75643450 

0.72189089 

0.68452741 

0.64244061 

0.61339577 

0.59133016 

555.34841856 

110.40559872 

54.79039104 

26.97926208 

17.70888608 

13.07722272 

10.29963480 

8.45965896 

6.14876216 

4.76772284 

2.97004524 

2.12266942+ 

1.41628788 

1.15474343 

1.01868394 

0.93620225- 

0.84385097 

0.78463336 

0.74585992 

0.65802047+ 

0.62643775 

0.60246871 


(sc  'd  1 


xsbnl 


.-•V  v 

iseevc.  4 
s  r s  «ev  3i  sv 
:n  w-  . .. f  V 

.:  -Vi 

;r  •  : 

£!•  ::.£  :  ..'8  l£V 
•  c; 

•.:  r .  .<  a 

rsa.2  :r 

-a  r  '  •  r 

<*>  *  ?  <■  *  c*. 

/  26  V  >  dT 


r  if.  ee\ .  .  v  i  i 

181 

eex 

;?££€$  J  -OSS 

eex 

e  es.e.  vac 

88X 

8P  V80,l> 

vex 

8t  £  <?V  i  -  V 

"628.3*1 

re  ;eie.ees 

£0£ 

r  se  ?vs.o*- 

a;  .V82 

fc  0  5 

: *  .188 

20S 

;.  •  -  Ur  .  3  rVi  : 

80S 

8CS 

i  sees: 

80S 

IIS 

;8  >3:<  r:  V>£V 

£15 

3r  >  88  ;  .VI88 

EIS 

80  jeev.jcvii 

US 

2IS 

..  •  v  aeei  vi 

ae&vseee.etoss 

VIS 

BIS 

o*#t  f.oe«.e£vac 

exs 

58  :•  ■’  9r  .  -.■*■■ 

OSS 

xss 

sss 

ESS 

frSS 

2SS 

•  <v:e  c  .fccv 

vss 

18  < 

avp^?  t  *.  2V..X 

ess 

XES 

*E£ 

2ES 

de8eece2.vi88 

3ES 

:?08i  m  r .  (i2vxx 

VES 

8ES 

sxvvaeex  .£a>vx 

ees 

r.ev  (  3  .xexes 

ots 

i  .  a  •  a^eevae 

x*s 

:r  .  3'  :  .V8  0/-  ■: 

s*s 

B-15 


FOR  INITIAL  TABLEAU  - 


OBJECTIVE  FUNCTION 

— 

-0.40152110E- 

02 

INDEP.  VECTOR 

= 

1 

0 . 987242E  00 

INDEP.  VECTOR 

— 

98 

0 . 2  89 7  34E-01 

INDEP.  VECTOR 

— 

343 

0 . 122462E-01 

INDEP  c  VECTOR 

= 

14 

0.106308E  00 

INDEP.  VECTOR 

S3 

109 

0 . 615342E-02 

INDEP.  VECTOR 

SS 

259 

0. 316666E-01 

INDEP.  VECTOR 

= 

477 

0 . 664858E-02 

INDEP.  VECTOR 

S3 

232 

0 . 22  0156E-01 

INDEP.  VECTOR 

= 

484 

0 . 34  3212E-02 

INDEP.  VECTOR 

— 

471 

0 . 404424E-01 

INDEP.  VECTOR 

s 

239 

0 . 771331E-02 

INDEP  o  VECTOR 

= 

384 

0 . 1306 18E-01 

INDEP »  VECTOR 

= 

145 

0.274734E-02 

INDEP.  VECTOR 

s 

245 

0.980893E  00 

INDEP.  VECTOR 

— 

394 

0 . 332443E-04 

INDEP.  VECTOR 

= 

382 

0 . 1 34  6  06E-01 

THIS  IS  THE  OPTIMUM  SOLUTION 

OBJECTIVE  FUNCTION  =  0.15957026E  01 

BASIC  VARIABLES  AND  VALUES  ARE  AS  FOLLOWS 


259  i 

(-17) 

0 . 11389451E 

00 

415  i 

(-173) 

0. 36296422E 

02 

232 

0 . 13209053E 

00 

367  i 

(-125) 

0. 9570035 3E 

01 

339 

(-97) 

0 . 23864708E 

00 

351  i 

(-109) 

0 . 22129839E 

01 

102 

0 . 21339985E 

01 

240 

0 . 2  816  7650E 

01 

243 

(-1) 

0. 43130310E 

00 

150 

0. 33762138E 

02 

11 

0 . 51209398E 

00 

151 

0 . 65132731E 

00 

120 

0. 10159683E 

00 

197 

0. 99317753E 

01 

478  - 

(-236) 

0. 38013893E 

00 

121 

0. 71483271E 

00 

Cf  x  : 

U  x  i 
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P-V-T  FITTING  PROGRAM  -  PART  II 


$IBFTC  CLARE  NODECK 


** 

DOUBLE  PRECISION  G ( 16 , 1 7 ) , S , D , B , SIG, DIV ( 16 ) , Y(16 ) 

** 

DIMENSION  X (16) 

4 

FORMAT  (1H  , 8HSINGULAR) 

5 

FORMAT  (1H  , 5HERROR) 

** 

N=16 

** 

R=10. 73147 

** 

READ (5,550)  TC 

★  * 

550 

FORMAT (1H  , F15 . 8) 

** 

52 

1 — 1 
It 

H 

CN 

in 

O 

Q 

** 

READ (5,53)  P,T,V,SG 

*  * 

53 

FORMAT ( 1H  ,3F15.8,F4.1) 

** 

AT=1 . 0/ ( V*V) 

** 

AU=AT/V 

*  * 

AV=AU/V 

'** 

AW=AV/V 

** 

AX=AW/V 

*  * 

EA=EXP (-5. 475*T/TC) 

•k  k 

G (I , 1) =0 . 01 

kk 

G  ( 1 , 2 ) =SG*AT/P 

kk 

G(I , 3) =SG*AU/P 

kk 

G (1 , 4)  =SG*AV/P 

kk 

G  ( 1 , 5 ) =SG*AW/P 

** 

G  ( 1 , 6 ) =SG*AX/P 

★  * 

G  (1 , 7) =SG*AT*T/P 

** 

G  (1 , 8 ) =SG*AU*T/P 

*  * 

G  ( 1 , 9 ) =SG*AV*T/P 

** 

G ( I , 10 ) =SG*AW*T/P 

kk 

G (I , 11) =SG*AX*T/P 

kk 

G ( I ,  12 ) =SG*AT*EA/P 

** 

G ( 1 , 1 3) =SG*AU*EA/P 

** 

G ( I  ,  14 ) =SG*AV*EA/P 

kk 

G ( I , 15 ) =SG*AW*EA/P 

** 

G ( I , 16 ) =SG*AX*EA/P 

** 

52 

G (I , 17) =SG* (P-R*T/V) /P 

NN=N+1 

DO  70  J=1 / N 

70 

DI V ( J) =ABS ( G ( 1 / J) ) 

DO  71  J=1,N 

DO  71  1=1 , N 

71 

G  (I ,  J)  =G  (I ,  J)  /DIV(  J) 

K=N 

60 

MAX=  1 

Q=ABS (G ( 1 , K) ) 

DO  65  1=2, K 

IF(ABS(G(I,K) ) . LT.Q)  GO  TO  65 

MAX=I 

Q=ABS ( G ( MAX , K ) ) 

65 

CONTINUE 

V\VA*WA 
r  WA«X A 

v  •<> 

q'  :ja*OC~;  .  , '. 


** 

** 

** 

w,.  «l  r; 

A  ,  1  '-f  C'Z 

:  '  .Tii . ,  iX,i)c  '■•A 
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DO  62  J=1,NN 
AST=G (MAX, J) 

G ( MAX , J ) =G ( K  ,  J ) 

62  G(K,J)=AST 

K=K- 1 

I F ( K . GE . 2 )  GO  TO  60 
9  K=N 

11  IF(G(K,K) .NE.0.0)  GO  TO  14 

16  M=K-1 

22  IF(G(M,K) .EQ.0.0)  GO  TO  20 

18  DO  19  J=1,K 
S=G(M,J) 

G  (M,  J)  =G  (K  ,  J) 

19  G(K,J)=S 
S=G (M,NN) 

G(M,NN) =G(K,NN) 

G (K ,NN) =S 

GO  TO  14 

20  M=M-1 
IF(M.LT.O)  GO  TO  99 
IF(M.GT.O)  GO  TO  22 

23  WRITE (6, 4) 

GO  TO  401 

14  D=G (K ,K) 

DO  26  J=1 ,K 

26  G(K,J)=G(K,J)/D 
G(K,NN)=G(K,NN)/D 
I=K-1 

30  IF(G(I,K) .EQ.0.0)  GO  TO  28 

27  B=G (I , K) 

DO  29  J=1 , K 

29  G(I, J)=G(I, J) -B*G(K, J) 

G  ( I , NN ) =G(I ,NN) -B*G(K,NN) 

28  1=1-1 
IF(I.LT.O)  GO  TO  99 
IF(I.GT.O)  GO  TO  30 

32  K=K-1 

IF(K.LT.l)  GO  TO  99 
IF(I.GT.l)  GO  TO  11 

34  IF(G(1,1) .NE.0.0)  GO  TO  36 

35  WRITE (6 ,4) 

GO  TO  401 

36  G ( 1 , NN ) =G ( 1 , NN ) /G ( 1 , 1 ) 

DO  40  1=2, N 

SIG=0 . 0 

NI=I-1 

DO  38  J=1 ,NI 

38  SIG=SIG+G(I, J) *G ( J ,NN) 

40  G ( I , NN) =G ( I , NN) -SIG 

DO  400  1=1, N 
Y (I) =G (I ,NN) 


400 


"  -  ,y 

Q? 

■  .  - 

CfiT>  s3 

•  v^.  .  n 
M  CY  OD 
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(Ktf ,  ••:  Dj  :  1  .-•  ■  ,  .)• 
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63 

55 

56 

****C 

57 


213 

1 


59 


1 

2 


61 

58 


212 

99 

401 


DO  63  1=1, N 
X  ( I ) =Y ( I ) /DI V ( I ) 

WRITE  (6,55) 

FORMAT  ( 1HO , 10X ,  1 6HFINAL  SOLUTIONS.  ) 

WRITE (6,56)  (X(I) ,1=1, N) 

FORMAT (1H  ,E15 . 8) 

CHECK  OF  SOLUTIONS. 

READ  (5,57)  MM 
FORMAT (1H  ,13) 

AM=MM 
S  R=  0  o  0 
SP=0. 0 
WRITE (6,213) 

FORMAT ( 1H  ,15HV(CU. FT./LB.M. )  , 6X, 9HT (DEG. R.  )  , 

9X ,  6HP  .  EXP  .  ,  8X ,  7HP  .  CALC .  ,  7X ,  8IIABS  .  DEV.  ,  7X ,  8HPER.  DEV.  ) 
DO  58  J=1 ,MM 
READ (5,59)  P,T,V 
FORMAT (1H  , 3F15.8) 

AT=1 o  0/ ( V*V) 

AU=AT/V 

AV=AU/V 

AW=AV/V 

AX=AW/V 

EA=EXP (-5. 475*T/TC) 

PK=R*T/V+X ( 2 ) *AT+X ( 3 ) *AU+X ( 4 ) *AV+X ( 5 ) *AW+X ( 6 ) *AX+X ( 7 ) * 
AT*T+X ( 8) *AU*T+X ( 9 ) *AV*T+X(10) *AW*T+X(11) *AX*T+X(12) * 
AT*EA+X (13) *AU*EA+X ( 14 ) *AV*EA+X(15) *AW*EA+X(16) *AX*EA 
DEV=ABS (P-PK) 

PERC=DEV*100 . 0/P 

WRITE (6, 61)  V, T , P , PK , DEV, PERC 

FORMAT ( 1H  , F15 . 4 , F15 . 2 , 4F15 . 4 ) 

SR=SR+DEV 
SP=SP+PERC 
SRA=SR/AM 
SPE=SP/ AM 

WRITE (6 , 212 )  SRA, SPE 

FORMAT (1H  , 9HAVE.DE V.  =  , FI 0. 6, 13HAVE . PER. DE' V.  =  , FI 0. 6) 

GO  TO  401 
WRITE (6, 5) 

CONTINUE 

STOP 

END 
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P-V-T  FITTING  PROGRAM  -  PART  I 


SAMPLE  OUTPUT 


OBJECT  PROGRAM  IS  BEING  ENTERED  INTO  STORAGE 

FINAL  SOLUTIONS 
0 . 1595 7249E  01 
-0.26280181E  05 
0. 9620744 3E  05 
-0.21501424E  06 
0 . 17530775E  06 
-0.47329448E  05 
0 . 2  72 6 9119E  02 
-0 . 12725418E  03 
0 . 2962 7500E  03 
-0 . 24136579E  03 
0.66861749E  02 
0. 92167559E  05 
-0. 49328657E  07 
0. 14323348E  08 
-0 . 13172478E  08 
0 . 38213556E  07 


V(ft  /lbm) 

T(°R) 

P.EXP. 

P.CALC. 

ABS.DEV. 

PER.  DEV. 

39.0604 

401.69 

99.1539 

100.7361 

1.5822 

1.5957- 

32.2940 

410.69 

120.8599 

122.4884 

1.6285 

1.3474 

26.9045 

419.69 

145.9019 

147.4381 

1.5363 

1.0529 

22.5647 

428. 69 

174.5885 

175.8349 

1.2464 

0.7139 

19.0349 

437. 69 

207.2577 

207.9322 

0.6745 

0.3254 

16.1333 

446.69 

244.2181 

244.0367 

0.1815 

0.0743 

13. 7230 

455.69 

285.8225 

284.5093 

1.3132 

0.4594 

11.7089 

464.69 

332.3794 

329.6126 

2.7668 

0.8324 

10.0064 

473.69 

384.2563 

379.8604 

4.3959 

1.1440 

8.5591 

482.69 

441.8352 

435.6334 

6.2018 

1.4037 

7.3197 

491.69 

505.4983 

497.4319 

8.0664 

1.5957+ 

6.1462 

500. 59 

575.6717 

570.8441 

4.8276 

0.8386 

5.1989 

509.69 

652.8698 

648.5577 

4.3121 

0.6605 

4.3869 

518.69 

737.6216 

733.3795 

4.2421 

0.5751 

3.6357 

527.69 

830.6914 

829.5576 

1.1338 

0.1365 

2.9131 

536.69 

933.0049 

939.2325 

6.2276 

0.6675 

1.5067 

547.56 

1070.6036 

1087.6871 

17.0835 

1.5957- 

582.5462 

401.69 

7.3480 

7.3593 

0.0113 

0.1534 

289.1793 

401.69 

14.6960 

14.7415 

0.0455 

0.3094 

142.5030 

401.69 

29.3920 

29.5622 

0.1702 

0.5791 

93.5638 

401.69 

44.0880 

44.4599 

0.3719 

0.8435 

69.0590 

401.69 

58.7840 

58.4292 

0.64  52 

1.0975 

54.3463 

401.69 

73.4800 

74.4377 

0.9577 

1.3033 

44.5049 

401.69 

88.1760 

89.5032 

1.3272 

1.5052 

635.9619 

437.69 

7. 3480 

7.3530 

0.0050 

0.0685 

316.3560 

437.69 

14.6960 

14.7149 

0.0189 

0.1285 

156.5530 

437.69 

29.3920 

29.4594 

0.0674 

0.2292 
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V ( ft  /lbm) 

T(°R) 

P.EXP. 

P . CALC . 

ABS.DEV. 

PER. DEV. 

103.2219 

437.69 

44.0880 

44.2479 

0.1599 

0.3627 

76.5881 

437.69 

58.7840 

59.0354 

0.2514 

0.4276 

60.5923 

437.69 

73.4800 

73.8374 

0.3574 

0.4863 

49.8979 

437.69 

88.1760 

88.6788 

0.5028 

0.5702 

36.4823 

437.69 

117.5680 

118.4503 

0.8823 

0.7505 

28.4809 

437.69 

146.9600 

147.8713 

0.9113 

0.6201 

23.0878 

437.69 

176.3520 

177.2767 

0.9247 

0.5244 

689.0886 

473.69 

7.3480 

7.3505 

0.0025 

0.0338 

343.2084 

473.69 

14.6960 

14.7046 

0.0086 

0.0582 

170.2577 

473.69 

29.3920 

29.4224 

0.0304 

0.1033 

112.5769 

473.69 

44.0880 

44.1592 

0.0712 

0.1614 

83.7506 

473.69 

58.7840 

58.8951 

0.1111 

0.1890 

66.4436 

473. 69 

73.4800 

73.6403 

0.1603 

0.2181 

54.8891 

473.69 

88.1760 

88.4057 

0.2297 

0.2605 

40.4090 

473.69 

117.5680 

118.0167 

0.4487 

0.3817 

31. 7660 

473.69 

146.9600 

147.4031 

0.4431 

0.3015 

25.9359 

473.69 

176.3520 

177.0209 

0.6689 

0.3793 

20.1269 

473.69 

220.4400 

220.9677 

0.5277 

0.2394 

14.2616 

473.69 

293.9200 

293.3238 

0.5962 

0.2028 

10.6169 

473.69 

367.4000 

364.9076 

2.4923 

0.6784 

742.0954 

509.69 

7.3480 

7.3492 

0.0012 

0.0158 

369.9338 

509.69 

14.6960 

14.6993 

0.0033 

0.0224 

183.8284 

509.69 

29.3920 

29.4049 

0.0129 

0.0439 

121.7838 

509.69 

44.0880 

44.1172 

0.0292 

0.0661 

90.7651 

509.69 

58.7840 

58.8296 

0.0456 

0.0776 

72.1468 

509.69 

73.4800 

73.5471 

0.0671 

0.0913 

59.7252 

509.69 

88.1760 

88.2755 

0.0995 

0.1129 

44.1876 

509.69 

117.5680 

117.7507 

0.1827 

0.1554 

34.8820 

509.69 

146.9600 

147.1364 

0.1764 

0.1200 

28.6359 

509.69 

176.3520 

176.6802 

0.3282 

0.1861 

22.3969 

509.69 

220.4400 

220.8504 

0.4104 

0.1862 

16.1298 

509.69 

293.9200 

294.2850 

0.3650 

0.1242 

12.3300 

509.69 

367.4000 

367.4251 

0.0251 

0.0068 

9.7639 

509.69 

440.8800 

439.7815 

1.0985 

0.2492 

6.4138 

509.69 

587. 8400 

581.4405 

6.3995 

1.0886 

795.0176 

545.69 

7.3480 

7.3485 

0.0005 

0.0070 

396.5606 

5  45.  69 

14.6960 

14.6972 

0.0012 

0.0079 

197.3145 

545.69 

29.3920 

29.3967 

0.0047 

0.0161 

130.8991 

545.69 

44.0880 

44.0972 

0.0092 

0.0209 

97.6950 

545.69 

58.7840 

58.7951 

0.0111 

0.0188 

77.7654 

545.69 

73.4800 

73.4962 

0.0162 

0.0220 

64.4697 

545.69 

88.1760 

88.2076 

0.0316 

0.0359 

47. 8535 

545.69 

117.5680 

117.6215 

0.0535 

0.0455 

37. 8852 

545.69 

146.9600 

147.0135 

0.0535 

0.0364 

31.2161 

545.69 

176.3520 

176.4969 

0.1449 

0.0821 

24.5542 

545.69 

220.4400 

220.6448 

0.2048 

0.0929 

17. 8781 

545.69 

293.9200 

294.1873 

0.2673 

0.0909 

13.8456 

545.69 

367.4000 

367.8763 

0.4763 

0.1296 

11.1386 

545.69 

440.8800 

441.5984 

0.7184 

0.1629 

7.7053 

545.69 

587.8400 

588.8692 

1.0292 

0.1751 

5.5763 

545.69 

734.8000 

735.0322 

0.2322 

0.0316 
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V ( f t  /lbm) 

T(°R) 

P.EXP. 

P.CALC. 

ABS.DEV. 

PER. DEV. 

3.9915 

545.69 

881.7600 

885.6750 

3.9150 

0.4440 

821.4470 

563.69 

7.3480 

7.3484 

0.0004 

0.0056 

409.8423 

563.69 

14.6960 

14.6969 

0.0009 

0.0064 

204.0329 

563.69 

29.3920 

29.3948 

0.0028 

0.0096 

135.4250 

563.69 

44.0880 

44.0941 

0.0061 

0.0139 

101.1282 

563.69 

58.7840 

58.7893 

0.0053 

0.0090 

80.5430 

563.69 

73.4800 

73.4875 

0.0075 

0.0102 

66.8102 

563.69 

88.1760 

88.1960 

0.0200 

0.0227 

49.6511 

563.69 

117.5680 

117.5990 

0.0310 

0.0264 

39.3515 

563.69 

146.9600 

147.0040 

0.0440 

0.0299 

32.4710 

563. 69 

176.3520 

176.4645 

0.1125 

0.0638 

25.5975 

563.69 

220.4400 

220.5895 

0.1495 

0.0678 

18.7099 

563.69 

293.9200 

294.1637 

0.2437 

0.0829 

14.5577 

563.69 

367.4000 

367.8980 

0.4980 

0.1355 

11.7730 

563.69 

440.8800 

441.8307 

0.9507 

0.2156 

8.2623 

563.69 

587. 8400 

589.8373 

1.9973 

0.3398 

6.1128 

563.69 

734.8000 

737.9580 

3.1581 

0.4298 

4.6091 

563.69 

881.7600 

888.2284 

6.4684 

0.7336 

2.4251 

563.69 

1175.6800 

1194.4406 

18.7606 

1.5957- 

1.0962 

563.69 

1469.6000 

1486.4647 

16.8647 

1.1476 

0.9002 

563.69 

2204.4000 

2207.6287 

3.2287 

0.1465 

0.8389 

563.69 

2939.2000 

2892.9321 

46.2679 

1.5742 

0.7698 

563.69 

4408.8000 

4435.6738 

26.8738 

0.6095 

0.7353 

563.69 

5878.4000 

5784.5967 

93.8033 

1.5957+ 

0.7085 

563.69 

7347.9998 

7276.1425 

71.8573 

0.9779 

0.6870 

563.69 

8817.5999 

8874.1005 

56.5006 

0.6408 

0.6596 

563.69 

11756.7999 

11630.5234 

126.2765 

1.0741 

0.6379 

563.69 

14695.9996 

14617.1074 

78.8923 

0.5368 

0.6212 

563.69 

17635.1998 

17569.6738 

65.5261 

0.3716 

0.6011 

563.69 

22043. 9994 

22122.2538 

78.2544 

0.3550 

0.5762 

563.69 

29391.9992 

29861.0312 

469.0317 

1.5958- 

423.1098 

581.69 

14.6960 

14.6970 

0.0010 

0.0066 

83.2994 

581.69 

73.4800 

73.4900 

0.0100 

0.0136 

40.7967 

581.69 

146.9600 

147.0273 

0.0673 

0.0458 

19.5136 

581.69 

293.9200 

294.2660 

0.3460 

0.1177 

12.3793 

581.69 

440.8800 

442.0022 

1.1223 

0.2545 

8.7839 

581.69 

587.8400 

590.1670 

2.3270 

0.3959 

6.5865 

581.69 

734.8000 

739.8288 

5.0288 

0.6844 

5.1068 

581.69 

881.7600 

888.5701 

6.8102 

0.7723 

3.1251 

581.69 

1175.6800 

1192.8554 

17.1754 

1.4609 

1.7455 

581.69 

1469.6000 

1491.6596 

22.0596 

1.5011 

1.0067 

581.69 

2204.4000 

2169.2267 

35.1733 

1.5956+ 

0.8995 

581.69 

2939.2000 

2892.2976 

46.9024 

1.5958+ 

0.8051 

581.69 

4408.8000 

4458.2197 

49.4197 

1.1209 

0.7621 

581. 69 

5878.4000 

5822.3418 

56.0582 

0.9536 

0.7296 

581.69 

7347.9998 

7345.8867 

2.1132 

0.0288 

0.7045 

581.69 

8817.5999 

8958.3134 

140.7135 

1.5958- 

0.6740 

581.69 

11756.7999 

11661.3974 

95.4025 

0.8115 

0.6503 

581.69 

14695.9996 

14537.2460 

158.7537 

1.0803 
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V ( f t  /lbm) 

T(°R) 

P.EXP. 

P.CALC. 

ABS.DEV. 

PER. DEV. 

0.6319 

581.69 

17635.1998 

17420.3476 

214.8523 

1.2183 

0.6100 

581.69 

22043.9994 

21834.1816 

209.8179 

0.9518 

0.5832 

581.69 

29391.9992 

29291.1952 

100.8042 

0.3430 

449.6097 

617.69 

14.6960 

14.6976 

0.0016 

0.0107 

88.7700 

617.69 

73.4800 

73.5100 

0.0300 

0.0409 

43.6518 

617.69 

146.9600 

147.0899 

0.1299 

0.0884 

21.0646 

617.69 

293.9200 

294.6250 

0.7050 

0.2399 

13.5284 

617.69 

440.8800 

442.3709 

1.4909 

0.3382 

9.7498 

617.69 

587. 8400 

590.3028 

2.4628 

0.4190 

7.4515 

617.69 

734.8000 

740.0801 

5.2801 

0.7186 

5.9253 

617.69 

881.7600 

888.8056 

7.0456 

0.7990 

3.9753 

617.69 

1175.6800 

1189.8388 

14.1588 

1.2043 

2.7804 

617.69 

1469.6000 

1489.9418 

20.3418 

1.3842 

1.3803 

617.69 

2204.4000 

2173.9406 

30.4594 

1 • 3818 

1.0631 

617.69 

2939.2000 

2892.6408 

46.5592 

1.5841 

0.8904 

617.69 

4408.8000 

4403.7583 

5.0417 

0.1144 

0.8206 

617.69 

5878.4000 

5865.3301 

13.0699 

0.2223 

0.7762 

617.69 

7347.9998 

7381.3676 

33.3678 

0.4541 

0.7444 

617.69 

8817.5999 

8920.0078 

102.4078 

1.1614 

0.7034 

617.69 

11756.7999 

11776.0625 

19.2626 

0.1638 

0.6755 

617.69 

14695.9996 

14541.8808 

154.1189 

1.0487 

0.6535 

617.69 

17635.1998 

17393.1640 

242.0359 

1.3725 

0.6281 

617.69 

22043.9994 

21692.2508 

351.7485 

1.5957+ 

0.5975 

617.69 

29391.9992 

28922.9684 

469.0308 

1.5958+ 

0.5747 

617.69 

36739.9992 

36352.9824 

387.0171 

1.0534 

476.0673 

653.69 

14.6960 

14.6987 

0.0027 

0.0185 

94.2054 

653.69 

73.4800 

73.5347 

0.0547 

0.0745 

46.4717 

653.69 

146.9600 

147.1580 

0.1980 

0.1348 

22.5802 

653.69 

293.9200 

294.8587 

0.9387 

0.3194 

14.6211 

653.69 

440.8800 

442.7273 

1.8473 

0.4190 

10.6380 

653.69 

587.8400 

590.7226 

2.8827 

0.4904 

8.2341 

653.69 

734.8000 

739.5775 

4.7775 

0.6502 

6.6296 

653.69 

881.7600 

888.5693 

6.8094 

0.7722 

4.6133 

653.69 

1175.6800 

1187.5926 

11.9126 

1.0133 

3.4107 

653.69 

1469.6000 

1483.0101 

13.4101 

0.9125 

1. 8590 

653.69 

2204.4000 

2207.1749 

2.7749 

0.1259 

1.2979 

653.69 

2939.2000 

2922.9897 

16.2103 

0.5515 

0.9968 

653.69 

4408.8000 

4368.1831 

40.6169 

0.9213 

0.8876 

653.69 

5878.4000 

5884.7830 

6.3830 

0.1086 

0.8283 

653.69 

7347.9998 

7362.3947 

14.3949 

0.1959 

0.7882 

653.69 

8817.5999 

8847.2087 

29.6088 

0.3358 

0.7353 

653.69 

11756.7999 

11785.2075 

28.4076 

0.2416 

0.7016 

653.69 

14695.9996 

14570.3588 

125.6409 

0.8549 

0.6755 

653.69 

17635.1998 

17462.9252 

172.2747 

0.9769 

0.6466 

653.69 

22043.9994 

21732.2988 

311.7007 

1.4140 

0.6086 

653.69 

29391.9992 

29860.9968 

468.9976 

1.5957- 

0.5875 

653.69 

36739.9992 

36194.7420 

545.2573 

1.4841 

502.5037 

689.69 

14.6960 

14.6998 

0.0038 

0.0257 
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'  ( f t^/lbm) 

T(°R) 

P.EXP. 

P . CALC . 

ABS.DEV. 

PER.  DEV. 

99.6125 

689.69 

73.4800 

73.5614 

0.0814 

0.1108 

49.2564 

689.69 

146.9600 

147.2499 

0.2899 

0.1973 

24.0677 

689.69 

293.9200 

295.0407 

1. 1207 

0.3813 

15.6715 

689.69 

440.8800 

443.1884 

2.3084 

0.5236 

11.4769 

689.69 

587.8400 

591.3456 

3.5056 

0.5964 

8.9531 

689.69 

734.8000 

739.9436 

5.1436 

0.7000 

7.2753 

689.69 

881.7600 

888.0238 

6.2638 

0.7104 

5.1681 

689.69 

1175.6800 

1185.0877 

9.4077 

0.8002 

3.9211 

689.69 

1469.6000 

1476.8746 

7.2746 

0.4950 

2.2799 

689.69 

2204.4000 

2205.0561 

0.6561 

0.0298 

1.5848 

689.69 

2939.2000 

2922.2157 

16.9843 

0.5779 

1.1223 

689.69 

4408.8000 

4388.0654 

20.7346 

0.4703 

0.9651 

689.69 

5878.4000 

5903.5646 

25.1647 

0.4281 

0.8861 

689.69 

7347.9998 

7363.2120 

15.2122 

0.2070 

0.8340 

689.69 

8817.5999 

8856.3151 

38.7152 

0.4391 

0.7698 

689.69 

11756.7999 

11727.7795 

29.0204 

0.2468 

0.7282 

689.69 

14695.9996 

14612.1960 

83.8037 

0.5702 

0.6981 

689.69 

17635.1998 

17500.5702 

134.6296 

0.7634 

0.6654 

689.69 

22043.9994 

21767.4702 

276.5293 

1.2544 

0.6271 

689.69 

29391.9992 

29007.3600 

384.6392 

1.3087 

0.6004 

689.69 

36739.9992 

36226.5076 

513.4917 

1.3976 

0.5803 

689.69 

44087.9988 

43384.4736 

703.5254 

1.5957+ 

528.9331 

725.69 

14.6960 

14.7004 

0.0044 

0.0302 

105.0126 

725.69 

73.4800 

73.5774 

0.0974 

0.1325 

52.0269 

725.69 

,.146.9600 

147.3159 

0.3559 

0.2422 

25.5270 

725.69 

293.9200 

295.2605  ■ 

1.3405 

0.4561 

16.7008 

725.69 

440.8800 

443.4717 

2.5917 

0.5879 

12.2877 

725.69 

587.8400 

591.8361 

3.9961 

0.6798 

9.6370 

725.69 

734.8000 

740.4007 

5.6007 

0.7622 

7.8886 

725.69 

881.7600 

887.0359 

5.2759 

0.5983 

5.6750 

725.69 

1175.6800 

1183.1112 

7.4312 

0.6321 

4.3631 

725.69 

1469.6000 

1474.5426 

4.9427 

0.3363 

2.6451 

725.69 

2204.4000 

2193.2287 

11.1713 

0.5068 

1. 8696 

725.69 

2939.2000 

2902.1612 

37.0388 

1.2602 

1.2640 

725.69 

4408.8000 

4401.5530 

7.2470 

0.1644 

1.0539 

725.69 

5878.4000 

5916.3138 

37.9138 

0.6450 

0.9496 

725.69 

7347.9998 

7395.0576 

47.0578 

0.6404 

0.8847 

725.69 

8817.5999 

8857.6470 

40.0471 

0.4542 

0.8058 

725.69 

11756.7999 

11728.6789 

28.1210 

0.2392 

0.7564 

725.69 

14695.9996 

14611.7118 

84.2878 

0.5735 

0.7219 

725.69 

17635.1998 

17473.7260 

161.4739 

0.9156 

0.6845 

725.69 

22043.9994 

21781.3154 

262.6841 

1.1916 

0.6424 

725.69 

29391.9992 

28999.7772 

392.2222 

1.3345 

0.6134 

725.69 

36739.9992 

36272.8332 

467.1660 

1.2715 

0.5913 

725.69 

44087.9988 

43671.7504 

416.2485 

0.9441 

555.3484 

761.69 

14.6960 

14.7010 

0.0050 

0.0341 

110.4056 

761.69 

73.4800 

73.5859 

0.1059 

0.1441 

54.7904 

761.69 

146.9600 

147.3499 

0.3899 

0.2653 

26.9793 

761.69 

293.9200 

295.3395 

1.4195 

0.4830 
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V ( ft  /lbm) 

T  ( °R) 

P.EXP. 

P . CALC . 

ABS.DEV. 

PER. DEV. 

17. 7089 

761.69 

440.8800 

443.7584 

2.8784 

0.6529 

13.0772 

761.69 

587.8400 

592.2655 

4.4255 

0.7528 

10.2996 

761.69 

734.8000 

740.6798 

5.8798 

0.8002 

8.4597 

761.69 

881.7600 

887.8070 

6.0470 

0.6858 

6.1488 

761.69 

1175.6800 

1182.0811 

6.4011 

0.5445 

4.7677 

761.69 

1469.6000 

1473.6428 

4.0428 

0.2751 

2.9700 

761.69 

2204.4000 

2183.0381 

21.3619 

0.9691 

2.1227 

761.69 

2939.2000 

2892.2982 

46.9017 

1.5957+ 

1.4163 

761.69 

4408.8000 

4382.2837 

26.5163 

0.6014 

1.1547 

761.69 

5878.4000 

5888.4962 

10.0962 

0.1718 

1.0187 

761.69 

7347.9998 

7431.2811 

83.2813 

1.1334 

0.9362 

761.69 

8817.5999 

8958.3100 

140.7101 

1.5958- 

0.8439 

761.69 

11756.7999 

11758.9757 

2.1758 

0.0185 

0.7846 

761.69 

14695.9996 

14707.7120 

11.7123 

0.0797 

0.7459 

761.69 

17635.1998 

17507.4580 

127.7419 

0.7244 

0.6580 

761.69 

29391.9992 

28922.9736 

469.0259 

1.5958+ 

0.6264 

761.69 

36739.9992 

36243.4256 

496.5737 

1.3516 

0.6025 

761.69 

44087.9988 

43834.0560 

253.9429 

0.5760 

AVE.  DEV.  =  50.261772 
AVE . PER. DEV.  =  0.566961 
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Example  Problem 

The  Part  I  program  is  the  actual  linear  programming 
portion  of  the  procedure.  Use  of  this  program  and  the  meaning 
of  the  output  is  best  illustrated  by  a  simple  example.  The 
statements  in  the  program  which  need  to  be  changed  for  a  par¬ 
ticular  fitting  problem  are  indicated  in  the  program  listing 
by  **. 

Data:  P.  -  T.  i  =  1,  30 

1  l 

Approximating  Function: 

2 

P  =  a  +  a,T  +  a~T 

o  1  z 

Objective : 

Minimize  the  maximum  absolute  deviation  in  P. 

M  =  30  (Number  of  data  points) 

N  =  3  (Number  of  equation  coefficients) 

Statements  in  the  program  listing  indicated  by  **  are 
replaced  by  the  following. 

COMMON  A(4,60),  B(4),  C(60),  F(60),  K(4), 

RX(4) ,  M,  MM,  NN,  NA ,  BB 
READ  (5,1)  M,N 
1  FORMAT  ( 2X ,  214) 

NA  =  N  +  1 

MM  =  M  +  M 

DO  102  I  =  1,  M 


T  -  -  <5 


* -Oii:  ..  or; .Li 


O' 


4  •  u;  jj-.  r  t  r 

;  ;  r  r  c  q  t-r)  V.  '  .*r.a  v  , 

.  !  lo  ■  .  • 

*  foe;  ;oifoni  pnidaiX  ai  h3ci&m&  -  i': 


r<n  fc  e‘i  *rfi  Y.c  f-’ 


O' "I  ,  A  TO  ,WH  vK  ,  (MX* 


(M  :  £)  TAMHC'? 


v  ,i  »  i  so;  j 
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READ  (5,101)  P , T 

101  FORMAT  (1H  ,  2F15.8) 

WRITE  (6,400)  I,  P,  T 

400  FORMAT  ( 1H  ,  13,  2F  15.8) 

F ( I)  =  P 
A ( 2 , 1 )  =  1.0 
A(3 , I)  =  T 

102  A ( 4 , I)  =  T**2 

DO  5  J  =  1,  M 

JM  =  J  +  M 

F  ( JM)  =  -F  ( J) 

A ( 1 , J)  =  1.0 
A ( 1 , JM )  =  1.0 

For  minimization  of  the  maximum  percentage  deviation 
in  P,  the  following  modifications  are  made: 

F  ( I )  =  1.0 
A(2 ,1)  =  1.0/P 
A( 3 , I)  =  T/P 
102  A( 4 , I)  =  T**2/P 

A ( 1 , J)  »  0.01 
A ( 1 , JM)  =  0.01 

Note:  The  COMMON  statement  for  the  subroutine  must  be  made 

identical  to  that  of  the  main  program. 

The  Part  I  program  will  give  the  following  output: 

1)  Input  data  listing. 

I  p  T  1<  I  s<  3  0 


,-j  (ioj  )  -  '  /i 


(a  .SI  It  .  II  ,  Kit  TSB.1OT  OOl 

T  *=  (I,€)A 

M  ,X  *  l  2  OC 
M  +  L  *  f'C* 

(m-  *  (ml)*? 

0,1  »  (L%JL) A 

q\O.X  *  (I,S)A 

X  0 . 0  8=r  (l.,I)A 

XO.O  -  (ML, I) A 

<  •  "  n.  i  •  :o  :■  -7  <: '  -  i  ■  -bi 

:  ituo  oXXo^t  sill  o'"  1  " 
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2)  A  solution  for  the  initial  tableau.  This  is  of  no 
interest  to  the  user. 

3)  The  optimum  solution. 

OBJECTIVE  FUNCTION  =  X.XXXX 
BASIC  VARIABLES  AND  VALUES  ARE  AS  FOLLOWS 
2  X.XXXX 

20  X.XXXX 

7  X.XXXX 

53  X . XXXX 

The  objective  function  value  is  the  maximum  absolute 
or  percentage  deviation  of  the  fit,  depending  upon  which 
criterion  was  used. 

The  basic  variable  indexes  indicate  the  data  points 
at  which  the  maximum  deviation  occurs.  Some  of  the  indexes 
will  be  greater  than  the  number  of  data  points.  For  these 
values  attach  a  minus  sign  and  add  the  number  of  data  points, 
i.e.  53  - *-53  +  30  =  -23 

This  index  corresponds  to  data  point  number  23.  The 
sign  of  the  index  is  required  for  Part  II  program. 

Part  II  Program: 

This  program  uses  the  results  of  the  Part  I  program 
to  solve  for  the  equation  coefficients.  The  program  is  simply 
a  Gaussian  Elimination  solution  method  for  a  set  of  linear 
equations.  The  final  portion  of  the  program  is  a  built  in 
check  of  the  fitting  procedure  and  need  not  be  included. 


-  : 


■  x  • 


if  ■■  y  •  '  '•  "  ;  •  •  ' 
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Statements  which  must  be  changed  for  the  particular  fitting 
problem  are  indicated  by  **.  The  portion  of  the  program  fol¬ 
lowing  ****  may  be  omitted  entirely. 

For  the  particular  example  problem  the  set  of  linear 
equations  to  be  solved  is : 


1.0 

1.0 

T 

2 

<t2)2 

X 

P2 

1.0 

1.0 

o 

CN 

EH 

(t20)2 

a 

o 

= 

P20 

1.0 

1.0 

T7 

(t7)  2 

al 

P7 

1.0 

-1.0 

-T2  3 

-<T23)2_ 

_  a2 _ 

_JP23_ 

N 

=  number  of 

coefficients 

+  1 

=  4 

The  subscripts  refer  to  the  indexes  of  the  data  point 
listing  from  the  Part  I  program.  The  reason  for  obtaining  the 
signs  of  the  controlling  data  points  is  obvious  in  the  above 
example . 

The  X  value  to  be  solved  for  is  the  maximum  deviation 
of  the  fit  and  should  correspond  to  the  value  obtained  from 
the  Part  I  program. 

If  minimizing  of  the  maximum  percentage  deviation 
was  used  in  the  Part  I  program,  the  set  of  equations  to  be 
solved  would  be, 


0.01  1.0/P2  T2/P2  (T2)2/P2 

X 

1.0 

0.01  1.0/P20  t20/p20  <t20)2/p20 

a 

o 

1.0 

0.01  1.0/P7  T7//p  7  (T7)2/P7 

al 

1.0 

0.01  -1.0/P23  _T2  3//P2  3  “(T23)  //P23_ 

_  a2 _ 

-1.0 

-  .  I 


X  «■  8  oi  :X£00  X  y  i&dmua 

v  io  ■- x  *  -oq  1  °';:1 

vr  1  a*  fcaeaa^oo  bic/orla  baa  ±0 

■  -  ..  _  — *•'  3.ri. 

ad  od  BfloiiBt/p©  So  iae  slfcJ  ,r«ip cxq  X  *****  i  l  beau  acw 


10. 0 

<J\S(,  ) 
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This  is  assuming  that  the  same  set  of  controlling  data 
points  were  obtained  from  the  linear  program  using  this  differ¬ 
ent  fitting  criterion. 
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P-V-T  FITTING  FOR  IMPROVED  THERMODYNAMIC  PROPERTIES 
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Part  I  Program 


Approximating  Function 


P  = 


RT 


V 


+  — —  + 

2 


A3 

V3  V4 


V' 


Ar 


V 


B,T  B9T  B^T 

+  —4—  +  -4— 


V 


V' 


V 


b4t 

v5 


BCT 

5 

v6 


Cx  -5.475Tr 

T  e 

V 


C2  -5 . 4  75Tr  C3  -5.475Tr 

“T  e  -4  e 

V  V 


C4  5.475Tr 

“5  e 
V 


-5.475T 


V 


R 


Derived  Function: 


(H-Ho)t  =  J*  [2 


A. 


V 


3  A, 


4  A. 


+ - j  + - 3 

2  V  3  V 


5  A4  6  A5 

4  V5  5  V6 


B,  T 

+  + 


B„T  B  B„T  BCT 

2  3  4  5 

+  — d—  +  — -r-  + 


2  3  4  5 

V  V  V  V  V 


-5.475T 


C1(2  +  5 . 475Tr) e 


R 


V 


-5 . 4  75T, 


C2(3  +  5 . 475Tr) e 


2  V* 


-5.475T 


+ 


C3  (4  +  5 . 4  75Tr) e 
_  _ 

3V° 


R 


b«vt*aC 

ft  £ 

7 

T,a 

©(  '  •»  2)  -3 


)  rar-  .f.  +  t 


4 


C-3 


-5.475T 

C4  (5  +  5.47  5Tr) e 


-5  475T 

Cr ( 6  +  5.475T  )e  "  R 

+  - =-5 -  ] 

5  V5 

J  -  units  conversion  factor. 

The  function  in  brackets  has  the  units  of  PV. 
Fitting  Criterion: 

Minimize  the  maximum  percentage  deviation  in  pressure 

subject  to  the  absolute  deviations  in  (H-H  )m  values  being 

o  T 


less  than 

or 

equal  to  2.0  BTU  per  pound. 

Input  Information: 

M 

— 

Total  number  of  data  points  of  both  types 

Ml 

— 

Number  of  P-V-T  data  points 

M2 

- 

Number  of  (H-H  ) T  data  points 

N 

— 

Number  of  equation  coefficients 

TC 

— 

Critical  temperature 

WT 

- 

Molecular  weight 

CONST 

— 

Specified  maximum  deviation  in  (H-Hq)t 

R 

— 

gas  constant 

CON 

- 

Conversion  factor  for  (H-Hq)t  in  units  of 

PV  to  BTU/lb . 

P 

— 

pressure 

T 

— 

temperature 

V 

— 

volume 

H 

— 

enthalpy  departure  (H-H  ) T 
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Output  Information: 

1)  Objective  value  -  maximum  percentage  deviation  of  the 
P-V-T  fit. 

2)  Indexes  of  the  controlling  data  points. 

For  indexes  greater  than  M,  attach  a  minus  sign  and  add 
M  to  find  the  corresponding  data  point  along  with  its 
sign.  Separate  the  data  points  into  P-V-T  and  (H-Hq)t 
groups . 

Controlling  data  points  are  the  ones  where  the  maximum 
deviations  occur.  Data  points  with  indexes  greater  than  M 
correspond  to 


>  ,  >  P 

calc. 


exp. 


(H-H  ) 


Part  II  Program 

This  program  is  basically  a  Gaussian  Elimination 
solution  of  a  set  of  linear  equations.  The  linear  equations 
to  be  solved  are  the  approximating  or  derived  functions  at 
the  controlling  data  points. 

Input  Information: 

N  -  Number  of  equation  coefficients  plus  one 

R  -  gas  constant 
TC  -  critical  temperature 
WT  -  molecular  weight 
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CONST 

— 

specified  maximum  error  in  (H-H^) T 

N1 

— 

number  of  P-V-T  controlling  data  points 

N2 

— 

N1  +  1 

CON 

conversion  factor  for  (H-H  )m  from  units  of 

o  T 

PV  to  BTU/lb 

*  *p 

— 

pressure 

*  *<ji 

— 

temperature 

*  *y 

— 

volume 

**SG 

— 

sign  of  the  index 

These 

values  are  for  the  controlling  P-V-T  data  points 

given 

by 

the  Part  I  program. 

— 

enthalpy  departure  (H-Hq)t 

*ip 

- 

temperature 

*V 

- 

volume 

*SG 

— 

sign  of  the  index 

*  These  values  are  for  the  controlling  (H-Hq)t  data  points 
given  by  the  Part  I  program. 

For  the  solution  checking  portion  of  the  program: 

M  -  total  number  of  data  points 
Ml  -  number  of  P-V-T  data  points 
Input  the  complete  P-V-T  and  (H-Hq)t  data  deck  from 
the  Part  I  program. 

Output : 

1)  The  maximum  percentage  deviation  in  pressure.  This  value 
should  equal  that  obtained  from  the  Part  I  program. 
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2)  The  equation  coefficients  in  the  order  which  they  appear 


in  the  equation. 

3)  A  check  of  the  calculated  pressures  and  enthalpy  departures. 
A  further  check  of  the  validity  of  the  solution  is  avail¬ 
able  by  noting  that  the  maximum  errors  in  pressure  and 
enthalpy  departure  occur  at  the  controlling  data  points. 
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P-V-T  AND  ENTHALPY  FITTING  PROGRAM  -  PART  I 


$IBFTC  MAIN  NODECK, NOLIST 

COMMON  A ( 16 , 550) ,B(16) ,C(550) ,F(550) ,K(16) ,RX(16) , 
1  M,MM,N,NA,BB 


** 

READ  (5,1)  M ,  Ml ,  M2  ,  N 

** 

1 

FORMAT (2X, 414) 

** 

READ (5,300)  TC,WT 

** 

300 

FORMAT ( 1H  , 2F15 . 8) 

** 

READ (5, 2)  CONST 

** 

2 

FORMAT (1H  , F10 . 4) 

** 

R=10. 73147 

** 

CON=144 . 0/ (WT*778 . 173 

** 

NA=N+1 

** 

M11=M1+1 

** 

MM=M+M 

** 

DO  102  1=1, Ml 

*  * 

READ (5,101)  P , T , V 

** 

101 

FORMAT (1H  , 3F15 . 8) 

★  * 

WRITE (6,400)  I,P,T,V 

*★ 

400 

FORMAT (1H  ,I3,3F15.8) 

** 

AT=1 . 0/ ( V*V) 

** 

AU=AT/V 

** 

AV=AU/V 

** 

AW=AV/V 

** 

AX=AW/V 

** 

EA=EXP (-5 . 475*T/TC) 

** 

F(I)  =  (P-R*T/V)  /P 

** 

A ( 2 , I ) =AT/P 

** 

A  (  3 , 1 ) =AU/P 

** 

A (4, I) =AV/P 

** 

A  ( 5 , 1 ) =AW/P 

** 

A ( 6 , I ) = AX/P 

** 

A ( 7 , I ) = AT*T/P 

*  * 

A ( 8 , I ) =AU*T/P 

*  * 

A(9,I)=A¥*T/P 

** 

a(io,i)=aw*t/p 

** 

A ( 11 , I ) =AX^T/P 

*★ 

A(12,I)=AT*EA/P 

*  * 

A(13,I)=AU*EA/P 

** 

A(14,I)=AV*EA/P 

** 

A ( 1 5 , I ) =AW*EA/P 

** 

102 

A ( 16 , I ) =AX*EA/P 

** 

DO  401  I=M11,M 

** 

READ (5,101)  H , T, V 

** 

WRITE (6,400)  I ,H , T , V 

** 

AT=CON/V 

** 

AU=AT/V 

** 

AV=AU/V 

*  * 

AW=AV/V 

** 

AX=AW/V 

<■  -o 
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*  * 
** 
** 
** 
** 
** 
★  * 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
** 
*  * 
** 
** 
** 
** 


EB=5.475*T/TC 
EC=EXP (-EB) 

ET= ( 2 . 0+EB) *EC 

EU=ET+EC 

EV=EU+EC 

EW=EV+EC 

EX=EW+EC 

F  (I ) =H 

A ( 2 , I ) =2 . 0*AT 
A (3, I) =1. 5*AU 
A  ( 4  f I ) =  4 . 0 * AV/ 3 . 0 
A(5,I)=1.25*AW 
A ( 6 , 1 ) =6 . 0  *AX/5 . 0 
A ( 7 , I ) =T*AT 
A ( 8 , I ) =T*AU 
A ( 9 , I ) =T*AV 
A (10 , I ) =T*AW 
A(ll , I ) =T*AX 
A (12 , I ) =ET*AT 
A ( 13 , I ) =EU*AU/2 . 0 
A ( 1 4 , I ) =E V*AV/ 3.0 
A ( 15 , I ) =EW*AW/4 . 0 
401  A ( 16 , I ) =EX*AX/5 . 0 

WRITE (6 ,3)  Ml, M2, CONST 

3  FORMAT  ( 1H0 ,  5X,  14,  23H  POINTS  LESS  THAN  LAMBDA, /6X, 
1  14,  17H  POINTS  LESS  THAN  ,  E15.6) 

DO  5  J=1 , Ml 
JM= J+M 
F(JM)=-F(J) 

A ( 1 , J) =0 . 0 1 
A(1 , JM) =0 . 01 
DO  5  1=2, NA 

5  A (I , JM) =-A (I , J) 

DO  6 J=M11 ,M 
JM= J+M 

F ( JM) =-F ( J) -CONST 
F ( J) =F ( J) -CONST 
A (1 , J) =0 . 00 
A(1 , JM) =0 . 00 
DO  6  1=2, NA 

6  A  ( I ,  JM)  =  -A  ( I ,  J) 

BB=0. 00 

K  ( 1)  =1 
B  ( 1 ) =1 
DO  7  1=2, NA 
XX=A (1,1) 

B (I) =-XX 
A  (1 , 1) =0 . 00 
DO  7  J=2 , MM 

7  A(I , J) =A(I , J) -XX*A ( 1 , J) 

DO  8  1=2, NA 
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IF  (B  ( I )  .GE. 0. 0)  GO  TO  8 
B ( I) =-B ( I ) 

DO  9  J=1,MM 
9  A  (I ,  J)  =~A  ( I  ,  J) 

8  CONTINUE 

DO  10  L=2 ,  NA 
XX=A (L ,  1 ) 

JA=1 

DO  11  J=2 , MM 

IF  (A(L,J) .LE.XX)  GO  TO  11 
XX=A (L, J) 

JA=  J 

11  CONTINUE 
K (L) = JA 
B(L)=B(L)/XX 
DO  12  J=1,MM 

12  A(L,  J) =A(L, J) /XX 
DO  13  1=1, NA 

IF  (I.EQ.L)  GO  TO  13 
YY=A(I, JA) 

B ( I ) =B ( I ) =B (L) *YY 
DO  14  J=1 , MM 

14  A(I , J) =A(I , J) -A(L, J) *YY 

13  CONTINUE 
10  CONTINUE 

40  DO  30  1=1, NA 

IF  (B (I) .GE.0.0)  GO  TO  30 
I  A=I 

GO  TO  50 
30  CONTINUE 
GO  TO  31 

50  B(IA)=-B(IA) 

DO  51  J=1 , MM 

51  A(IA, J) =-A(IA, J) 

XX=0. 00 

JA=0 

DO  52  J=1 , MM 
DO  53  1=1, NA 
IF  (J.EQ.K(I) )  GO  TO  52 
53  CONTINUE 

IF  (A(IA,J) .LE.XX)  GO  TO  52 
XX=A( IA, J) 

JA=  J 

52  CONTINUE 
K  ( I  A)  =  JA 
B  (IA)  =B  (IA)  /XX 
DO  54  J=1 ,MM 
A ( IA, J) =A (IA, J) /XX 
DO  55  1=1, NA 
IF  (I.EQ.IA)  GO  TO  55 
YY=A(I , JA) 


54 


e  o 


'  ... .  V  (0  .  ■'  .  ;  U  «  '  ■■  "  l 

(  xd)A  X 

1 1  or  co  (  x  .  ( i  ■  ' 7 

auoroio: 

i ..  x 

O  a  03  (J.  . '-U- 

yvMT  x) a*=  ( c  :: ) i 

or  or  oo  ((•  o.  II 

c  r  do 

:  s  a  * " 

,e  c  or  .  .  /■>'  "  ■  r 

ru  m  ■  loo 

.  y.  *.o  f  ;o  oa)  n 

(  C  '•  ■  •. 

xxv  (r  I .  i® (a 

i  c. 

dd  OL  OO  (AI.0a.3> 


C-10 


B  ( I )  =B  ( I )  -B  ( I  A)  *YY 
DO  56  J=1,MM 

56  A(I,  J) =A(I , J) -A(IA, J) *YY 

55  CONTINUE 
GO  TO  40 

31  DO  20  1=1, NA 

IA=K ( I ) 

20  BB=BB+F ( I A) *B  ( I ) 

DO  21  J=1 , MM 
C(J)=0. 00 

DO  22  1=1, NA 
IA=K (I) 

22  C(J)=C(J)+F(IA) *A(I,J) 

21  C(J)=C(J)-F(J) 

WRITE  (6,25) 

25  FORMAT  (1HO,  15X,  22H  FOR  INITIAL  TABLEAU  -  ) 

WRITE  (6,26)  BB 

26  FORMAT  (1HO,  16X,  21HOBJECTIVE  FUNCTION  =  ,  E15.8) 
WRITE  (6,27)  (K  (L)  ,  B(L),  L=1,NA) 

27  FORMAT  (1H,  18X,  16HINDEP .  VECTOR  =  ,  14,  5X,  E15.6) 
CALL  SIMPLEX 

67  CONTINUE 
STOP 
END 
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$IBFTC  DUAL  NODECK , NOLIST 
SUBROUTINE  SIMPEX 

COMMON  A(  16 , 550) ,B(16) ,C(550) ,F(550) ,K(16) ,RX(16) , 
1  M,MM,N,NA,BB 
399  QX=C ( 1) 

JA=1 

DO  300  J=2,MM 
IF  (C(J).GE.QX)  GO  TO  300 
JA=  J 
QX=C ( J) 

300  CONTINUE 

IF  (QX.GE.0.0)  GO  TO  330 
DO  301  1=1, NA 

IF  (A(I,JA) .LE.0.0)  GO  TO  302 
RX ( I ) =B (I)/A(I, JA) 

GO  TO  301 

302  RX ( I ) =-l . 0  0 

301  CONTINUE 

DO  303  1=1, NA 

IF  ( RX ( I ) . LT .0,0)  GO  TO  303 
QA=RX (I ) 

II=I 

GO  TO  304 

303  CONTINUE 
WRITE  (6,307) 

307  FORMAT  (1HO,  15X,  2 OH  UNBOUNDED  SOLUTION.  ) 

GO  TO  398 

304  IF  (II. GE . NA)  GO  TO  310 

IC=II+1 

DO  305  I=IC ,NA 

IF  ( RX ( I ) . LT . 0 . 0 . OR. RX ( I ) . GT . QA)  GO  TO  305 
IF  (RX(I) . EQ.QA)  GO  TO  306 
Q  A=  RX  ( I ) 

II=I 

GO  TO  305 

306  IF  (K(II) .GE.K(I) )  GO  TO  305 

QA=RX ( I ) 

II=I 

305  CONTINUE 

310  DX=A ( II , JA) 

B (II) =B (II) /DX 
DO  311  J=1 , MM 

311  A ( II , J) =A ( II , J) /DX 
DO  312  1=1, NA 

IF  (I.EQ.II)  GO  TO  312 
DXX=A(I, JA) 

B ( I ) =B ( I ) -DXX*B ( II ) 

DO  313  J=1 , MM 
IF  (J.EQ.JA)  GO  TO  314 
A ( I , J) =A (I , J) -DXX*A (II , J) 

GO  TO  313 

314  A (I , J) =0 . 00 

313  CONTINUE 

312  CONTINUE 
BB=BB-C(JA) *B (II) 
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DO  315  J=1,MM 
IF  (J.EQ.JA)  GO  TO  315 
C(J)=C(J)-C(JA) *A(II ,  J) 

315  CONTINUE 

C ( JA) =0.00 
K  ( II ) = JA 
GO  TO  399 

330  WRITE  (6,331) 

331  FORMAT  (1HO,  30H  THIS  IS  THE  OPTIMUM  SOLUTION.  ) 
WRITE  (6,321)  BB 

321  FORMAT  (1H  ,  20X,  22H  OBJECTIVE  FUNCTION  =  ,  E15.8) 
WRITE  (6,322) 

322  FORMAT  (1HO,  25X,  4 8HBASIC  VARIABLES  AND  VALUES 
1  AS  FOLLOWS  ) 

WRITE  (6,323)  (K(I),  B(I),  1=1, NA) 

323  FORMAT  (1H  ,  29X,  14,  12X,  E15.8) 

398  CONTINUE 

RETURN 

END 
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P-V-T  AND  ENTHALPY  FITTING  PROGRAM  -  PART  I 

SAMPLE  OUTPUT 


Index  P  (psia) 


T  ( °  R) 


V  (ft'Vmole) 


OBJECT  PROGRAM  IS  BEING  ENTERED  INTO  STORAGE 


1 

14.69599984 

365.68798720 

265.42469632 

2 

293.91999744 

365.68798720 

11.59307728 

3 

440.87998976 

365.68798720 

7.05650000+ 

4 

587.83999488 

365.68798720 

4.69597540 

5 

734.79999488 

365.68798720 

3.17005208 

6 

881.75997952 

365.68798720 

2.00949060- 

7 

955.23999744 

365.68798720 

1.60472692 

8 

1028.71998464 

365.68798720 

1.34733608+ 

9 

1175.67998976 

365.68798720 

1.14444251 

10 

1469.59998976 

365.68798720 

1.00581990 

11 

1763.51995904 

365.68798720 

0.94773212+ 

12 

2057.43996928 

365.68798720 

0.90649673 

13 

2351. 35997952 

365.68798720 

0.87646567 

14 

2645.27992832 

365.68798720 

0.85330718 

15 

2939.19997952 

365.68798720 

0.83388487 

16 

3673.99993344 

365.68798720 

0.79763524 

17 

4408.79996928 

365.68798720 

0.77108083 

18 

5878.39995904 

365.68798720 

0.73296257 

19 

7347.99986688 

365.68798720 

0.70686787 

20 

8817.59993856 

365.68798720 

0 . 68577Q01 

21 

10287.19992832 

365.68798720 

0.66860211 

22 

11756.79991808 

365.68798720 

0.65509932 

23 

13226.39966208 

365.68798720 

0.64292556 

24 

14695.99973376 

365.68798720 

0.63243434- 

25 

14.69599984 

401.68798720 

291.92793088 

26 

293.91999744 

401.68798720 

13.25675200 

27 

440.87998976 

401.68798720 

8.34710928 

28 

587.83999488 

401.68798720 

5.86273680 

29 

734.79999488 

401.68798720 

4.34775612 

30 

881.75997952 

401.68798720 

3.31388952 

31 

1028.71998464 

401.68798720 

2.58871766 

32 

1175.67998976 

401.68798720 

2.06140524 

33 

2057.43996928 

401.68798720 

1.10119004- 

34 

2351.35997952 

401.68798720 

1.03003101 

35 

2645.27992832 

401.68798720 

0.97985900- 

36 

2939.19997952 

401.68798720 

0.94366154 

37 

3673.99993344 

401.68798720 

0.88001050 

38 

4408.79996928 

401.68798720 

0.83877044 

39 

5878.39995904 

401.68798720 

0.78354889 

40 

7347.99986688 

401.68798720 

0.74740716 

41 

8817.59993856 

401.68798720 

0.72062619 
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42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 
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P  (psia) 


T  (°R) 


V  (ft'Vmole) 


10287.19992832 

11756.79991808 
13226.39966208 

14695.99973376 

14.69599984 
440.87998976 

587.83999488 

881.75997952 

1175.67998976 

1469.59998976 
1763.51995904 
2057. 43996928 
2351. 35997952 
2645.27992832 
2939.19997952 
3673.99993344 
4408. 79996928 
5878.39995904 
7347.99986688 
8817.59993856 
10287.19992832 

11756.79991808 
13226.39966208 

14695.99973376 

14.69599984 
440.87998976 

587.83999488 
734.79999488 

881.75997952 

1175.67998976 

1469.59998976 
1763.51995904 
2057.43996928 
2351.35997952 
2645.27992832 
2939.19997952 
3673.99993344 
4408.79996928 
5878.39995904 
7347.99986688 
8817.59993856 

11756.79991808 

14695.99973376 

14.69599984 
293.91999744 

587.83999488 

881.75997952 

1175.67998976 

1469.59998976 


401.68798720 

401.68798720 

401.68798720 

401.68798720 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

446.68799488 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

491.68798720 

536.68799488 

536.68799488 

536.68799488 

536.68799488 

536.68799488 

536.68799488 


0.69904076 
0.68262781+ 
0.66791201 
0.65502898 
325.07481600 
9.77075192 
7.05045192 
4.32418172 
2.97806144 
2.20907598 
1.75432940 
1.48426528 
1.31888140 
1.20975196 
1 . 1317t>186 
1.01244239 
0.94077669 
0.85620947 
0.80389302 
0.76788105 
0.74026483 
0.71847803 
0.70077653 
0.68453772 
358.18600448 
11.10734784 
8.11739024 
6.32771400 
5.14056616 
3.67095884 
2.81283468 
2.26970526 
1.90785216 
1.66220692 
1.48885980 
1.36665870 
1.17255772 
1.06094695 
0.93737278 
0.86580721 
0.81839531 
0.75639931 
0.71637201 
391.18969344 
18.89150400 
9.13146064 
5.89243936 
4.28412168 
3.33775568 
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92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 
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P  (psia) 


T  ( °  R) 


V  (ft'Vmole) 


1763.51995904 

2057.43996928 

2351.35997952 

2645.27992832 

2939.19997952 

3673.99993344 

4408.79996928 

5878.39995904 

7347.99986688 

8817.59993856 

11756.79991808 

14695.99973376 

14.69599984 

440.87998976 

587.83999488 

881.75997952 

1175.67998976 

1469.59998976 

1763.51995904 

2057.43996928 

2351.35997952 

2645.27992832 

2939.19997952 

3673.99993344 

4408.79996928 

5878.39995904 

7347.99986688 

8817.59993856 

11756.79991808 

14695.99973376 

14.69599984 

440.87998976 

734.79999488 

881.75997952 

1175.67998976 

1469.59998976 

1763.51995904 

2057.43996928 

2645.27992832 

2939.19997952 

3673.99993344 

4408.79996928 

5878.39995904 

7347.99986688 

8817.59993856 

14695.99973376 

14.69599984 

440.87998976 

734.79999488 

881.75997952 


536.68799488 

536.68799488 

536.68799488 

536.68799488 

536.68799488 

536.68799488 

536.68799488 

536.68799488 

536.68799488 

536.68799488 

536.68799488 

536.68799488 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

581.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

671.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 


2.72385796+ 
2.30326792 
2.00327128 
1.78486420 
1.62034626 
1.35659228 
1.20129120 
1.02963979 
0.93581803 
0.87506099 
0.79696574 
0 . 74828191- 
424.26494976 
13.62519120 
10.10427568 
6.57848424 
4.83812480 
3.80960920 
3.13348352 
2.66281720 
2.32078662 
2.06411540 
1.86916558 
1.54390180 
1.34747550 
1.12898534 
1.01049919+ 
0.93443355+ 
0.84087380 
0.78244558 
490.16482816 
16.01052304 
9.50388288 
7.88229656 
5.86381976 
4.66240376 
3.86951816 
3.31007856 
2.58564800 
2.34158822 
1.91452956 
1.64570172 
1.33670738 
1.16592329 
1.05827516 
0.85085824 
556.02884608 
18.34925024 
10.97301888 
9.13582560 
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Index 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 

181 

182 

183 

184 

185 

186 

187 

188 

189 


C-16 


P  (psia)  T  (°R)  V  ( f t3/mole) 


1175.67998976 
1469.59998976 
1763.51995904 
2057.43996928 
2645.27992832 
2939.19997952 
4408.79996928 
5878.39995904 
7347.99986688 
8817.59993856 
11756.79991808 
14695.99973376 
14.69599984 
440.87998976 
734.79999488 
881.75997952 
1175.67998976 
1469.59998976 
1763.51995904 
2057.43996928 
2645.27992832 
2939.19997952 
3673.99993344 
4408.79996928 
5878.39995904 
7347.99986688 
8817.59993856 
11756.79991808 
14695.99973376 
-33. 50000000 
-48.80000000 
-68.30000000 
“108.10000000 
-125.80000000 
-27.90000000 
-39.19999968 
-51.59999936 
-81.69999936 
-104.89999872 
-23.90000000 
-33.00000000 
-42.50000000 
-65.60000000 
-86.69999872 
-20.90000000 
-28.50000000 
-36.19999968 
-55.19999936 
-72.80000000 


761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

761.68799232 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

851.68798720 

399.68999424 

399.68999424 

399.68999424 

399.68999424 

399.68999424 

429.68999424 

429.68999424 

429.68999424 

429.68999424 

429.68999424 

459.68999424 

459.68999424 

459.68999424 

459.68999424 

459.68999424 

489.68999424 

489.68999424 

489.68999424 

489.68999424 

489.68999424 


6.82635104 

5.45714104 

4.54881148 

3.90691176 

3.06377968 

2.77638364 

1.93807244 

1.54613626 

1.32816581 

1.18995544 

1.02198280 

0.92392119 

621.85705984 

20.66648912 

12.39846096 

10.33026000 

7.75485800 

6.21642176 

5.19736348 

4.47304640 

3.51544812 

3.18037816 

2.59401392 

2.21265500 

1.75368486- 

1.49012178 

1.32103884 

1.11693771 

0.99784366- 

5.64929032 

3.78499320 

2.65582174 

1.34928770+ 

1.09298807 

6.45917832 

4.49586588 

3.35531672 

1.85024428+ 

1.35572630 

7.19107024 

5.10695504 

3.89865168 

2.33520766 

1.65392418 

7.88701552 

5.70335096 

4.39865784 

2.72028216 

1.94810168 
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169  POINTS  LESS  THAN  LAMBDA 
20  POINTS  LESS  THAN  0.200000E  01 

FOR  INITIAL  TABLEAU  - 


OBJECTIVE  FUNCTION 

= 

-0. 

12188398E  00 

INDEP.  VECTOR 

= 

1 

0 . 999904E  00 

INDEP.  VECTOR 

= 

174 

0 . 370885E-03 

INDEP.  VECTOR 

= 

366 

0 . 980051E-02 

INDEP.  VECTOR 

— 

179 

0 . 646026E-03 

INDEP.  VECTOR 

~ 

175 

0 . 182  36  8E-01 

INDEP.  VECTOR 

= 

213 

0 . 977551E-03 

INDEP.  VECTOR 

= 

172 

0 . 487458E-02 

INDEP.  VECTOR 

— 

362 

0 . 126763E-02 

INDEP.  VECTOR 

aa 

188 

0. 406641E-02 

INDEP.  VECTOR 

= 

291 

0 . 158537E-02 

INDEP.  VECTOR 

ss 

374 

0. 800862E-02 

INDEP.  VECTOR 

as 

378 

0 . 128549E-02 

INDEP.  VECTOR 

ss 

152 

0. 278642E-02 

INDEP.  VECTOR 

— 

359 

0. 813789E-02 

INDEP.  VECTOR 

sr 

42 

0.1624  9  8E-0  2 

INDEP.  VECTOR 

as 

225 

0 . 261378E-02 

THIS  IS  THE  OPTIMUM  SOLUTION. 


OBJECTIVE  FUNCTION  = 

BASIC  VARIABLES  AND 

0 . 76712749E  00 

VALUES  ARE  AS 

FOLLOWS 

3 

0 . 95600888E 

00 

358  (-169) 

0. 26743441E 

02 

291  (-102) 

0. 17338647E 

00 

222  (-33) 

0.37807111E 

01 

8 

0 . 27551992E 

01 

213  (-24) 

0.22545391E 

00 

43 

0.13031078E 

00 

11 

0. 52058789E 

01 

91 

0. 46948989E 

00 

354  (-165) 

0. 30890546E 

00 

173 

0. 46342028E- 

03 

178 

0. 37755528E- 

03 

118 

0 . 60112164E 

01 

195  (-6) 

0 . 18879286E 

01 

117 

0. 35123013E 

02 

224  (-35) 

0 . 16229045E 

02 
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P-V-T  AND  ENTHALPY  FITTING  PROGRAM  -  PART  II 


$IBFTC  CLARE  NODECK 

DOUBLE  PRECISION  G ( 16 , 17) , S , D,B , SIG, DIV ( 16 ) , Y ( 16 ) ,Q , AST 
DIMENSION  X ( 16 ) 

4  FORMAT (1H  , 8HSINGULAR) 

5  FORMAT  (1H  ,  5HERR0R) 

N=16 

R=10a  73147 

READ (5,550)  TC,WT, CONST 

550  FORMAT ( 1H  ,3F15.8) 

READ (5,551)  N1,N2 

551  FORMAT ( 1H  ,  213) 

CON=144 . 0/ (WT*778 .173) 

DO  52  1=1, N1 

READ (5,53)P,T,V,SG 
53  FORMAT ( 1H  , 3F15 . 8 , F4 . 1) 

AT=1 . 0/ ( V*V) 

AU=AT/'V 

AV=AU/V 

AW=AV/V 

AX=AW/V 

EA=EXP (-5. 475*T/TC) 

G (I , 1) =0 . 01 
G  ( 1 , 2 )  =  AT/P  *S  G 
G  ( 1 , 3 ) =AU/P  *SG 
G (I , 4 ) =AV/P*SG 
G(I,5)=AW/P*SG 
G  ( 1 , 6 ) = AX/P  *SG 
G  ( I , 7) -AT*T/P  *SG 
G(I , 8) =AU*T/P*SG 
G (I , 9) =AV*T/P*SG 
G ( 1 , 10 ) =AW*T/P  *SG 
G ( 1 , 11 ) =AX*T/P  *SG 
G (I , 12) =AT*EA/P*SG 
G (I , 13) =AU*EA/P*SG 
G(I ,14) =AV*EA/P*SG 
G(1 ,15) =AW*EA/P  *SG 
G (I , 16 ) =AX*EA/P*SG 
52  G (I , 17) = (P-R*T/V) /P*SG 
DO  300  I=N2,N 
READ (5,53)  H,T,V,SG 
AT=CON/V 
AU=AT/V 
AV=AU/V 
AW=AV/V 
AX=AW/V 
EB=5. 475*T/TC 
EC=EXP (-EB) 

ET= (2. 0+EB) *EC 
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EU=ET+EC 
EV=EU+EC 
EW=EV+EC 
EX=EW+EC 
G (I , 1) =0 . 0 
G ( I , 2 ) =2 . 0  *AT*SG 
G (I , 3) =1 . 5*AU*SG 
G  ( 1 , 4 ) =4 . 0 *AV/3 .  0  *SG 
G(I , 5) =1. 25*AW*SG 
G (I , 6) =6 . 0*AX/5 . 0*SG 
G  ( 1 , 7) =T*AT*SG 
G  (1 , 8) =T*AU*SG 
G  ( 1 , 9 ) =T*AV*SG 
G(I,10)=T*AW*SG 
0(1,11) =T*AX*SG 
G  ( 1 , 12 ) =ET*AT*SG 
G  ( 1 , 1 3) =EU*AU/2 ,0*SG 
G ( 1 , 14 ) =EV*AV/3 . 0  *SG 
G  ( 1 , 15 ) =EW*AW/ 4, 0  *SG 
G(I,16)=EX*AX/5. 0*SG 
300  G(I , 17) =H*SG-C0NST 
NN=N+1 
DO  70  J=1,N 

70  DIV ( J) =ABS (G ( 1 , J)  ) 

DO  71  J=1,N 

DO  71  1=1, N 

71  G ( I ,  J)  =G (I ,  J)  /DIV ( J) 

K=N 

60  MAX=1 

Q=ABS(G(1,K) ) 

DO  65  1=2, K 

IF(ABS(G(I,K) ) .LT.Q)  GO  TO  65 
MAX=I 

Q=ABS ( G (MAX, K) ) 

65  CONTINUE 

DO  62  J=1 ,NN 
AST=G (MAX , J) 

G (MAX, J) =G (K , J) 

62  G (K , J) =AST 

K=K- 1 

I F ( K . GE  c  2 )  GO  TO  6  0 
9  K=N 

11  IF(G(K,K) .NE.O. 0)  GO  TO  14 

16  M=K- 1 

22  IF(G(M,K) .EQ.0.0)  GO  TO  20 

18  DO  19  J=1 , K 

S=G (M, J) 

G (M, J) =G (K , J) 

G (K , J) =S 
S=G (M,NN) 

G(M,NN)=G(K,NN) 

G ( K , NN ) =S 
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23 

14 

26 

30 

27 

29 

28 

32 

34 

35 

36 

38 

40 

400 

63 

55 

56 

57 
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GO  TO  14 
M=M-1 

IF(MoLT.O)  GO  TO  99 
IF (M. GT . 0 )  GO  TO  22 
WRITE (6,4) 

GO  TO  401 
D=G (K  ,K) 

DO  26  J=1 , K 
G  (K ,  J) =G (K, J) /D 
G(K,NN)=G(K,NN)/D 
I=K-1 

IF  ( G  ( I , K)  * EQ .0.0)  GO  TO  28 
B=G(I,K) 

DO  29  J=1 , K 

G(I,  J)=G(I,J)=B*G(K,J) 

G ( I ,  NN) =G ( I , NN) -B*G (K ,  NN) 

1=1-1 

IF(I.LT.O)  GO  TO  99 
IF(I.GT.O)  GO  TO  30 
K=K-1 

IF(K.LT.l)  GO  TO  99 
IF (K . GT. 1)  GO  TO  11 
I F  ( G  ( 1 , 1 ) . NE .0.0)  GO  TO  36 
WRITE (6, 4) 

GO  TO  401 

G(1,NN)=G(1,NN)/G(1,1) 

DO  40  1=2, N 

S I G=  0  e  0 

NI=I- 1 

DO  38  J=1,NI 

S I G=S I G+G ( I , J ) *G ( J ,NN) 

G(I,NN)=G(I,NN)-SIG 
DO  400  1=1, N 
Y (I) =G ( I ,NN) 

DO  63  1=1, N 
X(I) =Y (I) /DIV(I) 

WRITE  (6,55) 

FORMAT  ( 1HO ,  10X,  16HFINAL  SOLUTIONS.  ) 

WRITE (6,56) (X(I) ,1=1, N) 

FORMAT ( 1H  ,E15 . 8) 

CHECK  OF  SOLUTIONS. 

READ (5,57)  MM, Ml 
FORMAT (1H  ,214) 

M11=M1+1 
AM=M1 
Al=MM-Ml 
SR=0 . 0 

WRITE (6,213) 

FORMAT (1H  , 15HV (CU . FT . /LB . M. ) , 6X, 9HT (DEG. R. ) , 

1  9X,  6HP  .  EXP  .  ,  8X,  7HP.CALC.  ,7X,8HABS  .  DEV.  ,  7X,  8HPER.  DEV. ) 
DO  58  J=1 , Ml 
READ (5,59)  P,T,V 
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59  FORMAT (1H  ,3F15.8) 

AT=1 . 0/ ( V*V) 

AU=AT/V 

AV=AU/V 

AW=AV/V 

AX=AW/V 

EA=EXP (-5 . 475*T/TC) 

PK=R*T/V+X ( 2 ) *AT+X ( 3) *AU+X(4) *AV+X(5) *AW+X(6) *AX+X(7) 

1  *AT*T+X ( 8 ) *AU*T+X ( 9 ) *AV*T+X(10) *AW*T+X(11) *AX*T+X(12) 

2  *AT*EA+X(13) *AU*EA+X ( 14 ) *AV*EA+X(15) *AW*EA+X(16) *AX*EA 
DEV=ABS (P-PK) 

PERC=DEV*100. 0/P 
WRITE (6,61)  V,T,P,PK,DEV,PERC 
61  FORMAT ( 1H  , F15 . 4 , F15 . 2 , 4F15 . 4 ) 

SR=SR+DEV 
58  SP=SP+PERC 
SRA=SR/AM 
SPE=SP/AM 

WRITE (6,212)  SRA, SPE 

212  FORMAT (1H  , 9HAVE . DEV. = , F10 . 6 , 13HAVE . PER. DEV. = , F10 . 6 ) 

S  R=  0 . 0 
SP=0.0 

WRITE (6,214) 

214  FORMAT ( 1H  , 15HV(CU. FT . /LB .M. ) , 6X, 9HT ( DEG. R. ) , 9X, 6HH . 

1  EXP. , 8X, 7HH. CALC. , 7X , 8HABS . DEV. , 7X , 8HPER. DEV. ) 

DO  301  J=M11 , MM 

RE AD (5, 59)  H,T,V 

AT=CON/V 

AU=AT/V 

AV=AU/V 

AW=AV/V 

AX=AW/V 

EB=5 .4  75  *T/TC 

EC=EXP (-EB) 

ET=(2. 0+EB) *EC 

EU=ET+EC 

EV=EU+EC 

EW=EV+EC 

EX=EW+EC 

HK=X ( 2 ) *2 . 0  *AT+X ( 3) *1.5*AU+X(4) *4 . 0 *AV/3 . 0+X ( 5 ) *1.25 

1  *AW+X ( 6 ) *6. 0*AX/5. 0+X (7) *T*AT+X(8) *T*AU+X(9) *T*AV+X(10) 

2  *T*AW+X ( 11 ) *T*AX+X (12) *ET*AT+X(13) *EU*AU/2 . 0+X ( 14 ) *EV 

3  *AV/3 . 0+X ( 15 ) *EW*AW/4 . 0+X ( 16 ) *EX*AX/5.0 
DE V=ABS (H-HK) 

PERC=DEV*100 . 0/ABS (H) 

WRITE  (6, 61)  V,  T ,  H  ,HK  ,  DEV ,  PERC 
SR=SR+DEV 
301  SP=SP+PERC 
SRA=SR/A1 
SPE=SP/Al 

WRITE (6,212)  SRAf  SPE 
GO  TO  401 
99  WRITE (6 ,5) 

401  CONTINUE 
STOP 
END 
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P-V-T  AND  ENTHALPY  FITTING  PROGRAM  -  PART  II 

SAMPLE  OUTPUT 


OBJECT  PROGRAM  IS  BEING  ENTERED  INTO  STORAGE. 

FINAL  SOLUTIONS. 

0 . 76716072E  00 
-0 . 104 85025E  05 
0.50240222E  04 
-0 . 48803061E  04 
0 . 32307928E  04 
-0. 1097780 3E  04 
0 . 12119481E  02 
-0 . 45553829E  01 
0 . 16332147E  02 
-0 . 15146344E  02 
0.7463576 IE  01 
-0 . 21888554E  06 
0 . 82097618E  06 
-0 . 77795798E  06 
-0.89290670E  05 
0 . 21571796E  06 


V  ( f  t/lbm) 

T(°R) 

P.EXP. 

P.CALC. 

ABS.DEV. 

PER.  DEV. 

265.4247 

365.69 

14.6960 

14.6905 

0.0055 

0.0371 

11.5931 

365.69 

293.9200 

292.3259 

1.5941 

0.5424 

7.0565 

365.69 

440.8800 

437.4977 

3.3823 

0.7672+ 

4.6960 

365.69 

587. 8400 

584.4301 

3.4099 

0.5801 

3.1701 

365.69 

734.8000 

735.1740 

0.3740 

0.0509 

2.0095 

365.69 

881.7600 

888.5245 

6.7645 

0.7672- 

1.6047 

365.69 

955.2400 

952.6669 

2.5731 

0.2694 

1.3473 

365.69 

1028.7200 

1020.8281 

7.8919 

0.7672+ 

1.1444 

365.69 

1175.6800 

1169.5712 

6.1088 

0.5196 

1.0058 

365.69 

1469.6000 

1478.7901 

9.1901 

0.6253 

0.9477 

365.69 

1763.5200 

1749.9913 

13.5286 

0.7671+ 

0.9065 

365.69 

2057.4400 

2043.4027 

14.0373 

0.6823 

0 .8765 

365.69 

2351. 3600 

2336.8584 

14.5016 

0.6167 

0.B533 

365.69 

2645.2799 

2626.0901 

19.1899 

0.7254 

0.8339 

365.69 

2939.2000 

2922.9564 

16.2436 

0.5527 

0.7976 

365.69 

3673.9999 

3653.4062 

20.5938 

0.5605 

0.7711 

365.69 

4408.8000 

4386.2751 

22.5248 

0.5109 

0.7330 

365.69 

5878.4000 

5868.4759 

9.9241 

0.1688 

0.7069 

365.69 

7347.9998 

7301.8981 

46.1017 

0.6274 

0.6858 

365.69 

8817.5999 

8811.6422 

5.9578 

0.0676 

0.6686 

365.69 

10287.1999 

10343.2196 

56.0197 

0.5446 

0.6551 

365.69 

11756.7999 

11786.4904 

29.6906 

0.2525 

0.6429 

365.69 

13226.3997 

13305.2336 

78.8339 

0.5960 
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V ( ft  /lbm) 

T(°R) 

P.EXP. 

P . CALC . 

ABS.DEV. 

PER. DEV 

0.6324 

365.69 

14695.9996 

14808.7418 

112.7422 

0.7672- 

291.9279 

401.69 

14.6960 

14.6964 

0.0004 

0.0025 

13. 2568 

401.69 

293.9200 

293.1151 

0.8049 

0.2739 

8.3471 

401.69 

440.8800 

438.4718 

2.4082 

0.5462 

5.8627 

401.69 

587.8400 

583.8883 

3.9517 

0.6722 

4.3478 

401.69 

734.8000 

730.1023 

4.6977 

0.6393 

3.3139 

401.69 

881. 7600 

878.9353 

2.8247 

0.3204 

2.5887 

401.69 

1028.7200 

1025.8795 

2.8405 

0.2761 

2.0614 

401.69 

1175.6800 

1174.4077 

1.2723 

0.1082 

1.1012 

401.69 

2057.4400 

2073.2240 

15.7840 

0.7672- 

1.0300 

401.69 

2351. 3600 

2366.5034 

15.1434 

0.6440 

0.9799 

401.69 

2645.2799 

2665.5737 

20.2938 

0.7672- 

0.9437 

401.69 

2939.2000 

2953.4209 

14.2209 

0.4838 

0.8800 

401.69 

3673.9999 

3687.5097 

13.5098 

0.3677 

0.8388 

401.69 

4408.8000 

4403.4374 

5.3626 

0.1216 

0.7835 

401.69 

5878.4000 

5862.1262 

16.2738 

0.2768 

0.7474 

401.69 

7347.9998 

7309.2213 

38.7785 

0.5277 

0.7206 

401.69 

8817.5999 

8766.1270 

51.4729 

0.5838 

0.6990 

401.69 

10287.1999 

10269.4102 

17.7897 

0.1729 

0.6826 

401.69 

11756.7999 

11666.6064 

90.1935 

0.7672+ 

0.6679 

401.69 

13226.3997 

13150.9639 

75.4358 

0.5703 

0.6550 

401.69 

14695.9996 

14666.8120 

29.1877 

0.1986 

325.0748 

446.69 

14.6960 

14.6967 

0.0007 

0.0044 

9.7708 

446.69 

440.880 

439.7150 

1.1650 

0.2642 

7.0505 

446.69 

587.8400 

585.2897 

2.5503 

0.4338 

4.3242 

446.69 

881.7600 

876.1652 

5.5948 

0.6345 

2.9781 

446.69 

1175.6800 

1166.9556 

8.7244 

0.7421 

2.2091 

446.69 

1469.6000 

1460.5191 

9.0809 

0.6179 

1.7543 

446.69 

1763.5200 

1757.8892 

5.6308 

0.3193 

1.4843 

446.69 

2057.4400 

2057.1726 

0.2674 

0.0130 

1.3189 

446.69 

2351.3600 

2354.0887 

2.7287 

0.1160 

1.2098 

446.69 

2645.2799 

2650.4379 

5.1579 

0.1950 

1.1318 

446.69 

2939.2000 

2951.1535 

11.9535 

0.4067 

1.0124 

446.69 

3673.9999 

3681.8965 

7.8965 

0.2149 

0.9408 

446.69 

4408.8000 

4415.2249 

6.4249 

0.1457 

0.8562 

446.69 

5878.4000 

5862.3538 

16.0462 

0.2730 

0.8039 

446.69 

7347.9998 

7337.7340 

10.2658 

0.1397 

0.7679 

446.69 

8817.5999 

8794.5684 

23.0315 

0.2612 

0.7403 

446.69 

10287.1999 

10273.0477 

14.1522 

0.1376 

0.7185 

446.69 

11756.7999 

11740.3426 

16.4573 

0.1400 

0.7008 

446.69 

13226.3997 

13183.8579 

42.5419 

0.3216 

0.6845 

446.69 

14695.9996 

14754.2316 

58.2319 

0.3962 

358.1860 

491.69 

14.6960 

14.6954 

0.0006 

0.0041 

11.1073 

491.69 

440.8800 

440.0976 

0.7824 

0.1775 

8.1174 

491.69 

587.8400 

586.3886 

1.4514 

0.2469 

6.3277 

491.69 

734.8000 

732.3523 

2.4477 

0.3331 

5.1406 

491.69 

881.7600 

877.8263 

3.9337 

0.4461 

3.6710 

491.69 

1175.6800 

1168.3064 

7.3736 

0.6272 

2.8128 

491.69 

1469.6000 

1459.2208 

10.3792 

0.7063 
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V ( ft  /lbm) 

T(°R) 

P.EXP. 

P . CALC . 

ABS.DEV. 

PER. DEV. 

2.2697 

491.69 

1763.5200 

1751.2328 

12.2872 

0.6967 

1.9079 

491.69 

2057.4400 

2047.5735 

9.8665 

0.4796 

1.6622 

491.69 

2351.3600 

2345.1412 

6.2188 

0.2645 

1.4889 

491.69 

2645.2799 

2646.3940 

1.1140 

0.0421 

1.3667 

491.69 

2939.2000 

2940.0527 

0.8527 

0.0290 

1.1726 

491.69 

3673.9999 

3681.4495 

7.4495 

0.2028 

1.0609 

491.69 

4408.8000 

4422.5828 

13.7828 

0.3126 

0.9374 

491.69 

5878.4000 

5877.9048 

0.4951 

0.0084 

0.8658 

491.69 

7347.9998 

7358.8032 

10.8034 

0.1470 

0.8184 

491.69 

8817.5999 

8831.6737 

14.0737 

0.1596 

0.7564 

491.69 

11756.7999 

11797.3793 

40.5795 

0.3452 

0.7164 

491.69 

14695.9996 

14751.4786 

55.4790 

0.3775 

391.1897 

536.69 

14.6960 

14.6967 

0.0007 

0.0046 

18.8915 

536.69 

293.9200 

294.0320 

0.1120 

0.0381 

9.1315 

536.69 

587.8400 

586.5437 

1.2963 

0.2205 

5.8924 

536.69 

881. 7600 

877.4296 

4.3304 

0.4911 

4.2841 

536.69 

1175.6800 

1168.1666 

7.5134 

0.6391 

3.3378 

536.69 

1469.6000 

1458.4344 

11.1656 

0.7598 

2.7239 

536.69 

1763.5200 

1749.9909 

13.5291 

0.7672+ 

2.3033 

536.69 

2057.4400 

2042.9281 

14.5119 

0.7053 

2.0033 

536.69 

2351.3600 

2338.3028 

13.0572 

0.5553 

1.7849 

536.69 

2645.2799 

2633.3880 

11.8919 

0.4496 

1.6203 

536.69 

2939.2000 

2930.4448 

8.7552 

0.2979 

1. 3566 

536.69 

3673.9999 

3662.2813 

11.7187 

0.3190 

1.2013 

536.69 

4408.8000 

4395.9807 

12.8193 

0.2908 

1.0296 

536.69 

5878.4000 

5857.9833 

20.4167 

0.3473 

0.9358 

536.69 

7347.9998 

7313.9586 

34.0412 

0.4633 

0.8751 

536.69 

8817.5999 

8772.9465 

44.6534 

0.5064 

0.7970 

536.69 

11756.7999 

11782.0466 

25.2467 

0.2147 

0.7483 

536.69 

14695.9996 

14808.7408 

112.7411 

0.7672- 

424.2649 

581.69 

14.6960 

14.6942 

0.0018 

0.0122 

13.6252 

581.69 

440.8800 

440.6240 

0.2560 

0.0581 

10.1043 

581.69 

587.8400 

586.7092 

1.1308 

0.1924 

6.5786 

581.69 

881.7600 

879.6659 

2.0941 

0.2375 

4.8381 

581.69 

1175.6800 

1170.7000 

4.9800 

0.4236 

3.8096 

581.69 

1469.6000 

1460.4964 

9.1036 

0.6195 

3.1335 

581.69 

1763.5200 

1751.9757 

11.5443 

0.6546 

2.6628 

581.69 

2057.4400 

2044.5743 

12.8657 

0.6253 

2.3208 

581.69 

2351.3600 

2338.8712 

12.4888 

0.5311 

2.0641 

581.69 

2645.2799 

2635.3995 

9.8805 

0.3735 

1.8692 

581.69 

2939.2000 

2929.9942 

9.2057 

0.3132 

1.5439 

581. 69 

3673.9999 

3663.7606 

10.2394 

0.2787 

1.3475 

581.69 

4408.8000 

4397.8073 

10.9927 

0.2493 

1.1290 

581.69 

5878.4000 

5849.6770 

28.7230 

0.4886 

1.0105 

581.69 

7347.9998 

7291.6287 

56.3711 

0.7672+ 

0.9344 

581.69 

8817.5999 

8749.9542 

67.6458 

0.7672+ 

0.8409 

581.69 

11756.7999 

11694.6655 

62.1344 

0.5285 

0.7824 

581.69 

14695.9996 

14751.2942 

55.2946 

0.3763 
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V ( ft  /lbm) 

T(°R) 

P.EXP. 

P . CALC . 

ABS.DEV. 

PER. DEV. 

490.1648 

671.69 

14.6960 

14.6959 

0.0001 

0.0007 

16.0105 

671.69 

440.8800 

441.6202 

0.7403 

0.1679 

9.5039 

671.69 

734.8000 

735.4067 

0.6067 

0.0826 

7.8823 

671.69 

881.7600 

882.0753 

0.3153 

0.0358 

5.8638 

671.69 

1175.6800 

1175.0034 

0.6766 

0.0576 

4.6624 

671.69 

1469.6000 

1467.5900 

2.0100 

0.1368 

3.8695 

671.69 

1763.5200 

1760.3747 

3.1453 

0.1784 

3.3101 

671.69 

2057.4400 

2053.9768 

3.4631 

0.1683 

2.5856 

671.69 

2645.2799 

2639.9325 

5.3475 

0.2022 

2.3416 

671.69 

2939.2000 

2931.3838 

7.8162 

0.2659 

1.9145 

,  671.69 

3673.9999 

3668.1355 

5.8645 

0.1596 

1.6457 

671.69 

4408.8000 

4407.2346 

1.5654 

0.0355 

1. 3367 

671.69 

5878.4000 

5869.6043 

8.7957 

0.1496 

1.1659 

671.69 

7347.9998 

7328.3234 

19.6764 

0.2678 

1.0583 

671. 69 

8817.5999 

8786.3510 

31.2489 

0.3544 

0.8509 

671.69 

14695.9996 

14789.5172 

93.5176 

0.6363 

556.0288 

761.69 

14.6960 

14.6967 

0.0007 

0.0046 

18.3493 

761.69 

440.8800 

442.0574 

1.1774 

0.2671 

10.9730 

761.69 

734.8000 

736.1498 

1.3498 

0.1837 

9.1358 

761.69 

881.7600 

882.6932 

0.9332 

0.1058 

6.8264 

761.69 

1175.6800 

1178.3572 

2.6772 

0.2277 

5.4571 

761.69 

1469.6000 

1472.2877 

2.6877 

0.1829 

4.5488 

761.69 

1763.5200 

1766.7704 

3.2504 

0.1843 

3.9069 

761.69 

2057.4400 

2060.7122 

3.2723 

0.1590 

3.0638 

761.69 

2645.2799 

2648.9890 

3.7090 

0.1402 

2.7764 

761.69 

2939.2000 

2940.9802 

1.7802 

0.0606 

1.9381 

761.69 

4408.8000 

4415.9431 

7.1431 

0.1620 

1.5461 

761.69 

5878.4000 

5901.3717 

22.9717 

0.3908 

1.3282 

761.69 

7347.9998 

7362.3222 

14.3224 

0.1949 

1.1900 

761. 69 

8817.5999 

8810.8502 

6.7498 

0.0765 

1.0220 

761.69 

11756.7999 

11759.3624 

2.5625 

0.0218 

0.9239 

761.69 

14695.9996 

14760.3882 

64.3885 

0.4381 

621. 8571 

851.69 

14.6960 

14.6973 

0.0013 

0.0087 

20.6665 

851.69 

440.8800 

442.0497 

1.1697 

0.2653 

12.3985 

851.69 

734.8000 

737.0665 

2.2665 

0.3084 

10.3303 

851.69 

881.7600 

884 .9898 

3.2298 

0.3663 

7. 7549 

851.69 

1175.6800 

1180.5149 

4.8349 

0.4112 

6.2164 

851.69 

1469.6000 

1475.9148 

6.3148 

0.4297 

5.1974 

851.69 

1763.5200 

1770.7612 

7.2412 

0.4106 

4.4730 

851.69 

2057.4400 

2065.7655 

8.3255 

0.4047 

3.5154 

851.69 

2645.2799 

2656.9557 

11.6758 

0.4414 

3.1804 

851.69 

2939.2000 

2956.9020 

17.7020 

0.6023 

2.5940 

851. 69 

3673.9999 

3699.1318 

25.1319 

0.6840 

2.2127 

851.69 

4408.8000 

4442.2686 

33.4686 

0.7591 

1.7537 

851.69 

5878.4000 

5923.4963 

45.0964 

0.7672- 

1.4901 

851.69 

7347.9998 

7401.7185 

53.7187 

0.7311 

1.3210 

851.69 

8817.5999 

8872.6271 

55.0272 

0.6241 

1.1169 

851.69 

11756.7999 

11821.3652 

64.5653 

0.5492 

0.9978 

851.69 

14695.9996 

14808.7414 

112.7417 

0.7672- 
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AVE . DEV . =  16 . 92149  7AVE .PER. DEV.=  0.369233 


V ( f t  /lbm) 

T(°R) 

H.EXP. 

H. CALC. 

ABS.DEV. 

PER.  DEV. 

5.6493 

399.69 

-33.5000 

-35.1763 

1.6763 

5.0040 

3.7850 

399.69 

-48.8000 

-49.8439 

1.0439 

2.1392 

2.6558 

399.69 

-68.3000 

-66.6788 

1.6212 

2.3737 

1.3493 

399.69 

-108.1000 

-110.1000 

2.0000 

1.8501+ 

1.0930 

399.69 

-125.8000 

-125.6945 

0.1055 

0.0839 

6.4592 

429.69 

-27.9000 

-29.1963 

1.2963 

4.6461 

4.4959 

429.69 

-39.2000 

-40.4390 

1.2390 

3.1607 

3.3553 

429.69 

-51.6000 

-52.0589 

0.4589 

0.8894 

1.8502 

429.69 

-81.7000 

-83.7000 

2.0000 

2.4480+ 

1. 3557 

429.69 

-104.9000 

-104.2450 

0.6550 

0.6244 

7.1911 

459.69 

-23.9000 

-25.0490 

1.1490 

4.8077 

5.1070 

459.69 

-33.0000 

-34.2858 

1.2858 

3.8963 

3. 8987 

459.69 

-42.5000 

-43.5769 

1.0769 

2.5339 

2.3353 

459.69 

-65.6000 

-66.9045 

1.3045 

1.9885 

1.6539 

459.69 

-86.7000 

-86.8768 

0.1768 

0.2039 

7.8870 

489.69 

-20.9000 

-21.9468 

1.0468 

5.0085 

5.7034 

489.69 

-28.5000 

-29.6650 

1.1650 

4.0879 

4.3987 

489.69 

-36.2000 

-37.5284 

1.3284 

3.6695 

2.7203 

489.69 

-55.2000 

-56.7617 

1.5617 

2.8292 

1.9481 

489.69 

-72.8000 

-73.9324 

1.1324 

1.5555 

AVE.DEV.  =  1.166175  AVE. PER. DEV.  =  2.690023 
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Example  Problem: 

Use  of  the  two  types  of  constraint  linear  program  is 
best  illustrated  by  a  simple  example.  Statements  to  be  changed 
for  a  particular  problem  are  indicated  in  the  program  listings 
by  **. 

Approximating  Function: 

2 

P  =  a  +  a,T  +  a~T 
o  1  2 

Derived  Property: 

dP 

—  =  a,  +  2a9T  =  DP 

dT  1  Z 


Data: 


i  =  1,  30 


i  =  1,  10 


Fitting  Constraints : 

Minimize  the  maximum  percentage  deviation  in  pressure 
subject  to  the  maximum  absolute  error  in  DP  being  less  than 
or  equal  to  some  specified  value. 


Part  I  Program 

M  =  40  (total  number  of  data  points) 

Ml  =  30  (number  of  P  data  points) 

M2  =  10  (number  of  DP  data  points) 

N  =  3  (number  of  function  coefficients) 

=  0.02  (arbitrary  maximum  error  for  DP) 


CONST 


:  Ci  bJ  d  c  j.  1  9 1  qmfiXtf 

srtqano  >  trc  3;  .  : 

I  V'  *i  G 

•r  - .  =  .  .  •.  v  •  r.  I  '  •  Xt  n 

'  -•  •  <  '  •  •  ••  ) 

,  adnioq  t'ie.  >  <?  '  •  ■ ;  0£ 

.  vb  <K  lo  ,t  t)  \ 

a?.  V.  '--o  vo^i'i  •••  *0 

■  . 7  V  -•*-  ' 
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The  program  statements  with  **  beside  them  can  be 
replaced  by  the  following. 

COMMON  A( 4 , 80) ,  B(4),  C(80),  F(80),  K(4),  RX(4),  M, 
MM,  N,  NA,  BB 
READ  (5,1)  M,  Ml,  M2,  N 

1  FORMAT  ( 2X,  414) 

READ  (5,2)  CONST 

2  FORMAT  (1H  ,F15.8) 

NA  =  N  +  1 

Mil  =  Ml  +  1 
MM  =  M  +  M 
DO  102  I  =  1,  Ml 
READ  (5,101)  P , T 
WRITE  (6,400)  I ,P , T 

101  FORMAT  (1H  ,2F15.8) 

400  FORMAT  ( 1H  ,I3,2F15.8) 

F ( I )  =  1.0 

A  (2, 1)  =  1.0/P 
A ( 3 , I )  =  T/P 

102  A ( 4 , 1)  =  T*  *2/P 

DO  401  I  =  Mil ,M 

READ  (5,101)  DP ,  T 
WRITE  (6,400)  I , DP , T 
F  ( I )  =  DP 

A ( 2 , 1 )  =  0.0 


(S,  .  )  U3M 

]'.  )  :  -iMHO 

i  +  i  ■ 

(IOI,fi)  <3  AS' 
GO,-,  )  -IT  .3  • 
(,:>.<  x  Hi':  tat* 

C  .r :  .  ct,  hh  ' 

q\o.l  -  (T  <  £ ' 

T,qa  U01,«)  QA:  v[r 
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A(3,I)  =  1.0 
401  A(4,I)  =  2 . 0*T 

Note:  The  COMMON  statement  for  the  subroutine  must  be  made 

identical  to  that  of  the  main  program. 

Required  Output  Information: 

THIS  IS  THE  OPTIMUM  SOLUTION 

OBJECTIVE  FUNCTION  =  X.XXXX 

BASIC  VARIABLES  AND  VALUES  ARE  AS  FOLLOWS 


50 

X.XXXX 

25 

X.XXXX 

75 

X . XXXX 

32 

X.XXXX 

For  basic  variable  indexes  greater  than  M(40) ,  attach 
a  minus  sign  and  add  M.  Also,  segregate  the  controlling 
data  points  into  their  respective  types. 


BASIC  VARIABLES 


P  points 


DP  points 


The  maximum  errors  in  P  and  DP  occur  at  the  data  points 
having  the  basic  variable  indexes  above.  The  maximum  per¬ 
centage  error  in  P  is  the  objective  value  and  the  maximum 
absolute  error  in  DP  is  equal  to  the  specified  constant  (0.02). 


0 ,1  iI%C)iA 


X6C  r  )/••• 

■  .•  •  .  ■  ,  ■  '  ■  .  •  r 


CB 

xxxx.x 

X  'X.X 

til  Iodine  o  r,iU  aJ&'2-:30-  y  :  *>•  r  c-  nE  - 
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Note:  If  the  specified  value  for  CONST  is  too  small  no 

solution  may  be  possible.  In  the  event  of  this  occurr¬ 
ing,  the  program  will  print  UNBOUNDED  SOLUTION  and 
cease  calculating.  This  print  out  may  also  occur 
if  a  mistake  is  made  in  setting  up  the  input  section 
(**)  of  the  program. 

Part  II  Program 

This  program  is  simply  Gaussian  Elimination  solution 
of  a  set  of  linear  equations.  A  section  is  also  included 
for  checking  the  data  fitting  procedure. 

The  set  of  linear  equations  to  be  solved  for  the 
example  problem  are ; 


0.01 

-1.°/P10 

“Tio/pio 

-(Tio>  /P10 

A 

-1.0 

0.01 

+1.0/P25 

+T25//P25 

+(T25)  /P25 

ao 

+1.0 

o 

• 

o 

o 

• 

o 

-1.0 

-2.0T35 

al 

-DP35-0.02 

o 

• 

o 

o 

• 

o 

+1.0 

+2 . 0T32 

a2 

+DP32-0. 02 

... 

The  use  of  the  indexes  and  their  signs  from  the  Part 
I  program  becomes  obvious  from  the  above  example.  They  are 
used  to  solve  for  the  equation  coefficients  by  the  above 
method. 

The  program  given  can  be  simply  modified  to  solve  the 
above  set  of  linear  equations  and  the  check  of  solutions 
section  can  be  modified  and  included  if  desired. 
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X  should  be  equal  to  the  maximum  percentage  deviation 
in  P  given  by  the  Part  I  program  if  the  equations  have  been 
solved  accurately. 
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EXAMPLE  ENTHALPY  DEPARTURE  CALCULATIONS 
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Program  I 

This  program  utilizes  the  given  variables  V  and  T  to 
calculate  P  and  (H-Hq)t  values.  For  the  proposed  equation 
of  state  both  pressure  and  enthalpy  departure  are  explicit 
functions  of  V  and  T,  therefore  the  calculations  are  simple. 

Program  II 

This  program  utilizes  the  most  useful  method  of  cal¬ 
culating  enthalpy  departures.  The  given  variables  are  P  and 
T  and  the  calculated  quantities  are  volume  and  enthalpy  de¬ 
parture.  The  proposed  equation  of  state  is  explicit  in 
pressure  and,  therefore,  must  be  solved  by  iteration  for 
volume.  The  calculated  volume  is  then  used  to  solve  explicit¬ 
ly  for  enthalpy  departure. 

A  search  interval  method  of  solving  for  volume  was 
used  rather  than  Newton-Raphson  since  (9P/9V)T  values  ap¬ 
proach  zero  in  the  critical  region.  The  method  proved  quite 
efficient  but  the  Newton-Raphson  method  is  undoubtedly  better 
provided  no  calculations  are  undertaken  in  the  critical 
region. 
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ENTHALPY  DEPARTURE  CALCULATION  PROGRAM 


TEMPERATURE  AND  VOLUME  GIVEN 


$IBFTC  CLARE  NOLI ST, NODECK 

READ (5,60)  N , L, WT ,TC 
60  FORMAT (1H  ,  2l3,2F15.8) 

AL=L 

R=10, 73147 

CON=144 . 0 (WT*778 . 173) 

DO  80  K=1 ,N 

READ (5, 2)  XI , X2 , X3 , X4 , X5 , X6 , X7 , X8 , X9 , XI 0 , Xll , X12 , 

1  X13 , X14 , X15 
2  FORMAT (1H  ,3E15.8) 

WRITE  (6,52) 

52  FORMAT ( 1HJ, 13HEQUATION  CONSTANTS) 

WRITE (6,2)  XI , X2 , X3 , X4 , X5 , X6 , X7 , X8 , X9 , XI 0 , Xll , X12 , 

1  X13 , X14 , X15 
WRITE (6,30) 

30  FORMAT (1HJ,4X,9HT (DEG. R. ) , IX, 12HV(CU. FT. /M. ) ,6X, 

1  7HP (PS I A) , 6X, 7HP . CALC. , IX, 12HABS . ER.  IN  P,1X, 

2  12HPER.ER.  IN  P , 2X , 11HH (BTU. /LB . ) , 6X , 7HH . CALC . , IX , 

3  12HABS.ER.  IN  H , IX , 12HPER. ER.  IN  H) 

SA=0 . 0 

SB=0 . 0 
SC=0. 0 
SD=0 . 0 
DO  20  1=1, L 
READ (5,1)  P , T ,EXDH, V 
1  FORMAT (1H  , 4F15 . 8) 

EA=EXP (-5  o  475*T/TC) 

PC=R*T/V  +X1/V  **2+X2/V  **3+X3/V  **4+X4/V  **5+X5/V 

1  *  *6+X6 *T/V  **2+X7*T/V  **3+X8*T/V  **4+X9*T/V  **5+XlO* 

2  T/V  **6+Xll*EA/V  **2+X12*EA/V  **3+X13*EA/V  **4+X14* 

3  EA/V  **5+X15 *EA/V  **6 
AT=CON/V 

AU=AT/V 

AV=AU/V 

AW=AV/V 

AX=AW/V 

EB=5 .4  75  *T/TC 

EC=EXP (-EB) 

ET= (2, 0+EB) *EC 

EU=ET+EC 

EV=EU+EC 

EW=EV+EC 

EX=EW+EC 
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DH=Xl*2 . 0*AT+X2*1. 5*AU+X3*4. 0*AV/3. 0+X4 *1 . 25 *AW+X5* 

1  6  0  0*AX/5 . 0+X6*T*AT+X7*T*AU+X8*T*AV+X9*T*AW+Xl0*T* 

2  AX+Xll*ET*AT+X12*EU*AU/2 . 0+X13*EV*AV/3 . 0+X14 *EW* 

3  AW/ 4 . 0+X15  *EX*AX/5 . 0 
ERH=ABS (EXDH-DH) 

PERH=ERH*100 . O/ABS (EXDH) 

ERP=ABS (P-PC) 

PERP=ERP*100 . 0/P 
SA=SA+ERP 

SB=SB+PERP 

SC=SC+ERH 

SD=SD+PERH 

20  WRITE (6,3)  T,  V,  P , PC, ERP , PERP ,EXDH , DH ,ERH, PERH 
3  FORMAT (1H  ,10F13.4) 

SA=SA/AL 

SB=SB/AL 

SC=SC/AL 

SD=SD/AL 

80  WRITE (6,70)  SA,SB,SC,SD 

70  FORMAT ( 1H J, 16HAVE . ABS . ER. IN  P= , F15 . 6/1X , 16HAVE . PER. ER. 

1  IN  P=,Fl5o6/lX,16HAVE.ABScER.IN  H= , F15 . 6/lX, 16HAVE . 
PER. ER. IN  H=/F15. 6) 

STOP 

END 


*0  ’  •  *  '  ■ 

n\o  co. 


Z 


D-5 


• 

CM 

H 

rH 

r-' 

CM 

00 

CN 

ft 

ft 

CO 

CO 

ft 

o 

ft 

00 

ft 

ft 

r- 

r- 

CN 

ft 

00 

rH 

rH 

cn 

CM 

rH 

ft 

00 

ft 

o 

CM 

ft 

CO 

ft 

CO 

CO 

o 

<T\ 

ft 

ft 

a 

cn 

00 

cn 

r-  rH 

00 

00 

00 

00 

r^- 

00 

ft 

CO 

CO 

r- 

o 

i"- 

o 

ft 

O'* 

• 

CM 

ft 

cn 

o 

ft 

cn 

ft 

CN 

ft 

00 

CD 

ft 

Ch 

o 

o 

ft 

r- 

00 

cn 

ft 

ft 

a 

ft 

H 

CN 

ft 

r" 

ft 

CO 

ft 

ft 

ft 

o 

O 

O 

o 

o 

rH 

CO 

in 

ft 

CN 

CN 

ft 

i — 1 

H 

• 

<n 

CO 

cn 

CM 

00 

o\ 

o 

00 

CM 

CN 

C" 

CD 

CM 

CD 

<n 

co 

rH 

rH 

00 

ft 

ft 

cn 

i — 1 

ft 

ft 

CTv 

CO 

ft 

ft 

ft 

CD 

O 

CN 

ft 

00 

ft 

rH 

r- 

f" 

w  a 

<T> 

ft 

cn 

ft 

CM 

ft 

00 

cn 

CO 

rH 

CO 

rH 

ft 

ft 

ft 

ft 

ft 

CM 

00 

• 

r^» 

rH 

o 

CO 

ft 

cn 

rH 

CO 

CO 

O' i 

.H 

rH 

ft 

O 

cn 

ft 

CO 

r- 

CD 

00 

ft  a 

• 

9 

H 

in 

CO 

ft 

ft 

rH 

CN 

r — 1 

ft 

o 

O 

o 

o 

O 

o 

rH 

rH 

rH 

rH 

rH 

OV 

CO 

cn 

CM 

00 

CTC 

O 

00 

CM 

00 

ft 

-3* 

CN 

ft 

i — 1 

ft 

cn 

cn 

CN 

• 

ft 

CT> 

rH 

ft 

ft 

O' \ 

CO 

ft 

ft 

CO 

ft 

O'* 

CN 

CO 

rH 

ft 

00 

CM 

CN 

u 

CT> 

ft 

CTi 

ft 

CM 

ft 

00 

cn 

CO 

rH 

ft 

CM 

00 

ft 

CO 

CO 

ft 

ft 

r- 

CM 

ON 

ft 

r- 

CN 

00 

ft 

CN 

ft 

00 

00 

o 

CO 

r-~ 

ft 

CO 

ft 

ft 

cn 

u 

G\ 

ft 

CM 

ft  i — i 

cn 

rH 

ft 

ft 

ft 

ft 

CM 

CM 

ft 

ft 

OC 

CO 

ft 

ft 

o 

• 

ft 

ft 

i — 1 

ft 

CM 

ft 

ft 

CO 

o 

CM 

ft 

ft 

CO 

00 

rH 

CN 

ft 

in 

*-r* 

l-M 

1 

I 

i 

rH 

1 - 1 

1 

1 

1 

1 

1 — 1 

1 

1 

1 

1 

1 

i 

1 

l 

1 

1 

i  i  i 


S 

a 

Bn 

D 

O 


W 

i-3 

ft 


CO 


W 

o 

g 

o 

H 

CO 

O 

a 

M 

Q 

8 

a 

Eh 

a 

W 

O 

a 

M 

a 

ft 

ft 

H 

s 

g 

o 


ft 

Eh 

U 

ft 

O 


m 


a 

in 

00 

ft 

H 

CO 

cr> 

2 

6 

7 

cn 

cn 

o 

ft 

co 

r- 

cn 

ft 

CM 

CN 

CO 

\ 

D 

ft 

00 

CO 

00 

ft 

r-~ 

cr>  rH  rH  ft 

ft 

ft 

CN 

ft 

CO 

o 

co 

CO 

ft 

CM 

Eh 

ft 

ft 

CO 

o 

CN 

CN 

ft  ft  00 

o 

CM 

ft 

ft 

co 

00 

CM 

CN 

ft 

ft 

f" 

a 

i 

1 

1 

H 

H 

l 

1  i  1 

H 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

a 

• 

o 

ft 

cn 

o 

CO 

H*  ft  CM 

cn 

00 

o 

f" 

o 

ft 

o 

ft 

ft 

r" 

CO 

a 

ft 

r" 

00 

cn 

Ch 

00  CO  CO 

CD 

CO 

ft 

CN 

ft 

ft 

o 

ft 

cn 

ft 

o 

a  a 

H 

00 

ft 

rH  CN 

ft 

ft  O  rH 

OC 

1 — 1 

ft 

H 

CO 

CN 

ft 

o 

co 

1 - \ 

ft 

• 

ft 

ft 

H 

ft 

CO 

00  rH  (Tl 

CN 

CO 

H 

co 

co 

ft 

r-« 

t" 

ft 

ft 

r- 

&  a 

a  H 

a 

o 

o 

o 

ft 

CN 

o 

0 

0 

0 

o 

o 

rH 

o 

o 

o 

o 

o 

o 

o 

o 

• 

cn 

cn 

00 

CM 

o 

CO 

CM 

r-~ 

co 

00 

ft 

ft 

CN 

CM 

o 

r~ 

CO 

CM 

CM 

ft 

H 

CN 

ft 

co 

a 

ft 

cn 

co 

ft 

VO 

VO 

ft 

ft 

r-' 

o 

o 

CO 

ft 

VO 

o 

CO 

CM 

H 

CN 

o 

o 

o 

o 

o 

a 

a 

ft 

cn 

CO 

00 

cn 

VO 

CO 

ft* 

ft 

o 

ft 

CN 

O 

o 

ft 

cn 

CM 

H 

• 

00 

o 

ft 

ft 

in 

cn 

CO 

o 

r- 

cn 

cn 

CM 

H 

cn 

ft 

ft 

CO 

00 

CM 

CD 

a 

a 

a 

a 

a 

a 

a 

ft 

CM 

ft 

VO 

CO 

cq 

H 

H 

ft 

H 

cn 

H* 

ft 

VO 

H 

ft 

ft 

ft 

cn 

CO 

cn 

O 

ft 

ft 

ft 

CO 

ft 

ft 

ft 

cn 

H 

00 

< 

rf 

ft 

I — 1 

H 

H 

cn 

ft 

CN 

in 

00 

cn 

ft 

VO 

00 

VO 

o 

H 

H 

00 

1 — 1 

cn 

cn 

00 

CN 

o 

VO 

CN 

r^- 

co 

00 

ft 

ft 

00 

00 

o 

CO 

00 

00 

00 

o 

CO 

ft 

cn 

• 

ft 

cn 

CO 

ft 

VO 

VO 

ft 

ft 

r- 

o 

o 

vo 

ft 

CO 

o 

co 

CM 

CO 

t"- 

CD 

ft 

CO 

VO 

CN 

u 

ft 

cn 

00 

00 

cn 

"S’ 

VO 

VO 

ft* 

ft 

o 

ft 

CM 

CM 

cn 

o 

ft 

cn 

00 

CM 

VO 

rH 

ft 

ft 

co 

o 

ft 

ft 

ft 

cn 

00 

o 

r- 

cn 

cn 

CN 

H 

o 

ft 

ft 

CO 

CD 

r- 

rH 

O 

o 

o 

o 

o 

u 

H 

ft 

H 

cn 

ft 

vo 

H 

ft 

ft 

ft 

cn 

CO 

o 

cn 

ft 

ft 

ft 

CO 

ft 

1 

• 

o 

o 

o 

ft 

o 

o 

o 

H 

o 

o 

o 

o 

cn 

00 

o 

o 

o 

<n 

00 

ft 

CN 

VO 

ft 

a 

CO 

00 

o 

ft 

O 

VO 

00 

o 

ft 

o 

vo 

00 

o 

ft 

cn 

CO 

00 

o 

ft 

cn 

O 

o 

O 

o 

o 

H 

H 

CN 

H 

H 

CN 

H 

H 

l — I 

H 

H 

H 

a 

a 

a 

a 

a 

C"- 

CM 

vo 

CN 

r"- 

00 

H1 

ft 

ro 

ft 

ft 

00 

r-' 

ft* 

cn 

VO 

cn 

< 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CO 

cn 

o 

CM 

cn 

ft 

ft 

ft 

ft 

H 

VO 

ft 

ft 

vo 

ft 

ft 

CN 

VO 

cn 

a 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

H 

co 

r> 

H 

ft 

00 

ft 

00 

co 

a 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

00 

ft 

ft 

ft 

VO 

H* 

CM 

CO 

00 

o 

ft 

O 

VO 

00 

o 

ft 

o 

vo 

00 

o 

ft 

o 

co 

00 

o 

ft 

o 

• 

• 

• 

• 

ft 

00 

o 

cn 

ft 

a 

1 — 1 

rH 

CN 

H 

H 

CN 

H 

H 

CM 

H 

H 

CN 

H 

o 

CN 

ft 

H 

ft 

ft 

ft 

H 

cn 

eh 

o 

o 

o 

o 

o 

— 

ii 

II 

II 

II 

% 

1 

ft 

ft 

1 

CN 

1 

H 

ft 

S 

ft 

o 

CO 

co 

O 

CM 

cn 

ft 

CM 

r^- 

» — 1 

o 

r>- 

CN 

cn 

o 

ft 

f" 

CO 

H 

a 

a 

a 

a 

E-i 

O 

o 

o 

o 

o 

cn 

ft 

ft 

cn 

ft 

cn 

in 

ft 

O 

ft 

f — 1 

OO 

ft 

CO 

CO 

00 

o 

00 

cn 

ro 

H* 

00 

ft 

H* 

cn 

ft 

cn 

ft 

ft 

ft 

cn 

o 

cn 

CO 

ft 

00 

o 

cn 

CN 

ft 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

Eh 

CO 

ft 

CO 

ft 

o 

ft* 

ft 

00 

ft 

H 

H 

CD 

CO 

vo 

00 

CO 

r" 

cn 

w 

H 

H 

H 

o 

00 

KT 

CN 

CN 

ft 

a 

a 

a 

a 

a 

u 

(N 

ft 

CO 

ft 

ft 

" — 

ft 

ft 

CN 

H 

H 

co 

ft* 

ft 

H 

H 

r" 

in 

CO 

CM 

rH 

["> 

ft 

ft 

CM 

H 

O 

cn 

o 

ft 

ft 

> 

a 

a 

a 

a 

00 

ft 

CN 

cn 

ft 

• 

• 

• 

• 

o 

CO 

H 

CM 

VO 

H 

a 

a 

a 

a 

H 

o 

ft 

O 

cn 

ft 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

a 

a 

a 

a 

Eh 

ft 

cn 

ft 

rH 

ft 

a 

VO 

VO 

VO 

VO 

VO 

co 

co 

VO 

vo 

co 

vo 

vo 

vo 

vo 

CO 

co 

CO 

vo 

vo 

vo 

< 

a 

< 

a 

< 

VO 

H 

CM 

H 

ft 

o 

• 

ft 

• 

• 

V - - 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

O'* 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

a 

a  a 

a 

o 

o 

o 

O 

o 

o 

<n 

<n 

cn 

cn 

cn 

CM 

CM 

CN 

CM 

CN 

ft 

ft 

ft 

ft 

ft 

CD 

00 

00 

00 

00 

> 

>  > 

a 

1 

1 

i 

1 

EH 

ft 

ft 

ft 

ft 

ft 

ft* 

ft* 

ft* 

ft* 

ft* 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

< 

< 

< 

< 

X-* 

QJ 

rrv 

00 

’ o 

CO 

A 

on 

• 

fii 

IT. 

CO 

% 

•  - 

> 

n 

■ 

» 

A 

y~i 

r 

cv 

V 

• 

■ 

m 

r  tt 


:  i  »  ^  l  l 


o»  '-1  t-)  "X  .  o 

r  '  .  r'  "■> 


,  -  i  O  CC  % T 

4."  1  iB  H»  O  C  V'  •  ■  c: 

(i 


r  :  1  O  r 


,  CO  *t'  ts  >-»  O'  ’ 

J  1 


W 


)  '  ■  v  O  CO  C  "I 

....  v-.  c  v  a/h  -•  o  .j 


OJCt  t.'  1  <  C-  i  -  t  «  •  Cv  >  •’  W  C  C  ' 

r.««*  *  •  ...  *  '  * 


< 


: 


-  5  J 

:>  )  c  ■  > 

►« 

W  fc  <  *  w 

C>  o.1  ^  ^  fu 

on  v  -.  r?  m  *r> 

,o  ru  V 

OJ 

>  C  1  > 

«t  .;  i 

VO  Q.  0  >  tr  V 

1  ‘  ■  ■  ■> 

c  O  O  O  O 

-  .  tv  H  V 

o  o  o  o  o 

t  «i  -a  u4  ui 

vj  a i  oo 

r>  •  o> 

-  ■  -in  r>  -  j 

t/o  cv  vo  c  a  n 

<V  C  •  f- 

p  "l  'n  **• 

>  *  ♦  W 

-  .  .*>(..  o»4 

w 

'I  h  i  Ol  Y' 

o  o  o  <  •  »*3 

W  M  j-  l'  ’ 

on  n  i  *v  *» 

.  O'- 

po  ,  • :  r 

•  IP.-;'  •  v  : 

: 

M  VWH V  O 

O  ,t>  O  "!  O  M 

- 

i  w-  c.1 

•  •  •  •  •  G 

o  ©  ©  t  >0 

i  .  ;  ... 

D-6 


Pi* 

W  S3 

Pi*  53 

W  M 

ph 


c^(Nh-HcnHr'^,o^r^y3c^inc^Ln(T\inHio 

no^nonwoa^cD^MnnconcohounLn 

ontNiocn^u3cot)'(NO(Tirooooocovo(Nin 

O  H  <n  03  OWHOO^UJODOOlOCrtlNOOVOCOLn 


intMfNHO'jnotMO^fOfNHOin^nNH 


Pi 

W  S3 

• 

CO  53 


3 


H 


roCTocMoinrQOCToooouoCToinoocoo^rr'-M4 

^nHOinvDCTio3omounvo^y3U)incoHfM 

h^Mooo\mnoLn^cohoh^iD(NUDM 

VOOVDOHfMMrrouJHOJOmHOHnmH 


HHHtMOrlHOOJOHHHHOrlHHHH 


u 

3 

u 

a 


n<j)cx)Oinnoo^oooiAo^incocoo'jtN^ 

vonooorrvDOicooincMn^^u)ioincoH(N 

is^ir'ocr(crifoinorj'^,cor^or'^u)fs)u)n 

HOO^OrUOr »-*34Ot''-rMOCMUOCTOOOCTOVOUOr''CTO 


inc7i\ooir)a\ocjnrrin^ruou)HCiM^^n 
n^u>H(NCM^incoo(Nfo^vDcooj(Nnint^ 
III*— >•— llltlr— lllllllllll 
I  I  I 


CQ 

A 

\ 

D 

H 

CQ 
' — ^ 

a 


• 

pi 

w  Pi 


p« 


CM 

CM 

VO 

VO 

.ER 

P 

o 

O 

O 

o 

o 

cn  s 

W 

W 

w 

w 

W 

SH 

rH 

rH 

M4 

r* 

VO 

VO 

00 

rH 

<T\ 

o 

M4 

ro 

VO 

r- 

ro 

co 

VO 

rH 

o 

rH 

M* 

co 

00 

rH 

rH 

o 

uo 

• 

00 

CM 

in 

CM 

rH 

u 

rH 

rH 

CO 

CM 

a 

o 

o 

o 

o 

O 

U 

1 

m4 

m4 

I 

CM 

VO 

in 

• 

Pi 

o 

o 

O 

o 

o 

W 

W 

W 

w 

W 

CM 

ro 

r- 

M4 

co 

CM 

o 

in 

VO 

< 

CM 

00 

rH 

in 

VO 

H 

O 

r* 

CM 

00 

o 

U) 

i"* 

ro 

00 

co 

PM 

CM 

oo 

ro 

00 

CM 

o 

o 

VO 

rH 

CTO 

Pi 

in 

rH 

rH 

CM 

CO 

CO  • 

Eh  o 

o 

o 

o 

o 

3  in 

1 

M4 

rH 

1 

rH 

1 

VO 

• 

s 

Eh  O 

O 

o 

o 

o 

CO 

53  W 

w 

w 

W 

w 

ro 

Eh 

O  m 

00 

oo 

o 

P-i 

L>  CM 

CM 

CM 

VO 

co 

S-> 

o 

CTO 

CO 

> 

53  m 

r* 

ro 

in 

uo 

O  oo 

o 

in 

ro 

CO 

M  ^ 

ro 

uo 

VO 

r- 

Eh  o 

CM 

in 

Pi 

rH 

ro 

O 

D  • 

• 

• 

• 

• 

O  o 

o 

o 

o 

o 

W  1 

» 

1 

Eh 

UO 

00 

ro 

rH 

00 

CO 

CM 

VO 

t"' 

CO 

CTO 

o 

uo 

VO 

r-" 

CTO 

LT0 

CM 

CM 

00 

CO 

00 

oo 

00 

uo 

r- 

CO 

rH 

i — f 

M4 

ro 

ro 

CM 

uo 

VO 

o 

00 

VO 

UO 

CM 

ro 

rr 

VO 

o 

CM 

CM 

ro 

uo 

00 

O 

CM 

ro 

M4 

VO 

00 

CM 

CM 

ro 

UO 

t"- 

1 

l 

1 

rH 

1 

rH 

1 

1 

1 

i 

1 

rH 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

rH 

UO  M4 

CTO 

VO 

ro 

o 

in 

r" 

ro 

ro 

VO 

M4 

00 

r- 

rH 

M4 

CTO 

o 

o 

O 

vo  ro 

OO 

UO 

r- 

CTO 

VO 

00 

CO 

r» 

CM 

00 

o 

vo 

CO 

rH 

ro 

UO 

ro 

co  uo 

rH 

CM 

o 

ro 

CM 

o 

ro 

1 - 1 

rH 

oo 

rH 

CTO 

o 

CTO 

r" 

VO 

ro  rH 

VO 

m 

CM 

VO 

CM 

VO 

ro 

o 

ro 

o 

CM 

CM 

ro 

VO 

M4 

o 

o  o 

3. 

2. 

0. 

o 

o 

1. 

0. 

0. 

0. 

0. 

*0 

o 

0. 

0. 

0. 

0. 

0. 

CTO 

rH 

rH 

r- 

r* 

CM 

CTO 

rH 

t"- 

00 

CTO 

O 

r"> 

CTO 

vo 

uo 

VO 

O 

CM 

ro 

CM 

CO 

CM 

"S’ 

rH 

uo 

ro 

O 

M4 

rH 

CO 

rH 

00 

uo 

oo 

00 

ro 

00 

rH 

^J4 

rH 

VO 

CO 

vo 

vo 

00 

VO 

ro 

o 

CO 

ro 

CM 

ro 

r- 

rH 

uo 

CM 

in 

o 

CTO 

r- 

rH 

in 

CO 

CM 

rH 

vo 

CTO 

ro 

ro 

o 

M* 

uo 

ro 

ro 

rH 

rr 

ro 

rH 

vo 

CTO 

CM 

rH 

O 

ro 

o 

rH 

1 — 1 

CM 

ro 

O 

CTO 

U0 

in 

rH 

rH 

rH 

1 - 1 

rH 

CTO 

CTO 

r-» 

ro 

00 

rH 

rH 

r" 

CM 

oo 

o 

ro 

vo 

1 — \ 

CO 

in 

CTO 

r- 

UO 

ro 

CM 

vo 

r- 

in 

rH 

UO 

VO 

o 

vo 

CO 

UO 

rH 

00 

rH 

O 

rH 

CM 

VO 

00 

rH 

uo 

CO 

ro 

00 

vo 

uo 

vo 

1 — 1 

ro 

VO 

CTO 

VO 

vo 

r"* 

VO 

CM 

00 

M1 

CM 

uo 

CTO 

o 

CM 

1 - 1 

in 

1 — 1 

CM 

oo 

ro 

o 

vo 

VO 

CTO 

uo 

M4 

VO 

vo 

OO 

vo 

CO 

ro 

CTO 

CM 

00 

o 

vo 

oo 

CO 

00 

VO 

CTO 

O 

CTO 

CTO 

CTO 

in 

in 

CTO 

CTO 

CTO 

rH 

i — i 

oo 

o 

CTO 

OO 

CTO 

CTO 

CTO 

CO 

CO 

CTO 

UO 

r- 

CTO 

uo 

o 

UO 

!"• 

CTO 

in 

o 

U0 

oo 

CTO 

OO 

UO 

I-" 

CTO 

CTO 

rH 

CM 

rH 

CM 

rH 

1 — 1 

rH 

rH 

,0 

o  o 

O' 

,0 

O 

,  0 

,0 

o 

,0 

,  0 

O 

,  0 

,0 

o 

,  0 

,  0 

,0 

o 

o 

o 

o  o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o  o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

vo 

00  o 

in 

o 

VO 

oo 

o 

uo 

o 

VO 

00 

o 

uo 

o 

vo 

00 

o 

uo 

o 

rH 

l — 1 

CM 

rH 

rH 

CM 

rH 

1 — 1 

CM 

rH 

rH 

CM 

MAO 
00  CO  VO 
OOI  H 
n  in  vo 
vo  H  vo 
N  >  H 
•  •  • 

O  O  rH  CM 
i — I 


II  II  II  II 


ro 

o 

00 

ro 

O 

CM 

OO 

ro 

CM 

> — 1 

o 

CM 

oo 

o 

M4 

l-" 

ro 

rH 

PU 

Pi 

H 

H 

oo 

uo 

uo 

oo 

ro 

CTO 

in 

uo 

O 

uo 

rH 

CO 

uo 

ro 

r* 

ro 

00 

o 

00 

CO 

in 

CTO 

UO 

oo 

in 

in 

uo 

o 

o 

CTO 

ro 

UO 

CO 

o 

CTO 

CM 

M4 

S 

2 

vo 

t—  VO 

ro 

o 

M4 

ro 

CO 

ro  rH 

rH 

00 

ro 

vo 

00 

ro 

f" 

CTO 

l-H 

H 

H 

H 

in 

ro 

CM 

rH 

rH 

vo 

'O4 

ro 

« — I 

rH 

UO 

ro 

CM 

rH 

UO 

'34 

CM 

rH 

Pi 

Pi 

Pi 

Pi 

W 

W 

W 

w 

U) 

Pi 

CO 

CO 

OO 

CTO 

oo 

oo 

CTO 

CTO 

oo 

CTO 

CTO 

CTO 

oo 

CTO 

oo 

CTO 

CTO 

CTO 

CTO 

CTO 

CTO 

CTO 

p} 

H 

PQ 

CQ 

vo 

VO 

vo 

vo 

vo 

vo 

vo 

vo 

vo 

VO 

vo 

VO 

vo 

VO 

VO 

vo 

VO 

VO 

VO 

VO 

< 

Pi 

< 

3 

CTO 

CTO 

CTO 

oo 

oo 

CTO 

CTO 

CTO 

CTO 

OO 

CTO 

CTO 

oo 

OO 

CTO 

CTO 

CTO 

OO 

CTO 

CTO 

H  W 

w 

w 

CO 

OO 

CTO 

CTO 

CTO 

CM 

CM 

CM 

CM 

CM 

in 

U0 

UO 

U0 

in 

00 

00 

00 

00 

00 

>  > 

£ 

> 

ro 

ro 

ro 

ro 

ro 

"34 

M4 

•^f4 

M4 

M4 

M* 

M4 

M4 

M4 

M4 

M4 

M4 

M4 

< 

< 

< 

< 

689979 


ft  »'  to 


n  l-i  ^1  '.r.  a 

-  r  m  o  < 

>  ">  c 

c.'  O  U'  n  M  p 

■  •  c.  t  to  to  ru 


(.  i  T  h  t  ^  ■ 


rn  ->  "» 
tr  T- 

tO  -i  l-*  - 


*  A< 


V  ■  -T 

'  • '  .  ^  "1 

O-  c.»  .  V*t 

.  J  n  r  j 

,  ‘S*  f> 


..  y  5  ;mon  "1  t  >  O  '  >  ^  ° 


5  O  .  v  C  >  O  < 


c:  QQ  >  i  ••  •  :-i>  co  - 1 

1  l ':  1  '  V  *>  (  •  CV 


-0  o  ‘  v  ft-  •  ;  O  O  t — »  tO  AJ  -  •  ■  :•">  *V  t-~  ’ 


'  O  "'1  i~T  >--» 

6/ 

OJ  -"i  (Vi 

ru 

u 

to 

» 

3  ftJ  fO 

ro  *'■  ru 

1— * 

o 

r> 

c.» 

o 

ru 

.13  ■ .  c  ■ 

r~* 

H* 

vO  fit 


O  O  o  3  0 

ca  w  i  « 

-■»  >- 

o  (V  o  r  o 

>  c-  •*>  • 

ft1  O  »-*  CO 

to  mun 

IO  05H  H  *  ' 
„  '  O  O  t -> 

c  :■  o 


tv  ru  v  c  t^ 
tv  ru  t~  to 

r»  >  ro  --t  v 

v  ••  v%  co  m  v«  to 


ft> 

V.' 

I 

t~,  tv  o  o 

v-  t~,  t-  ru  > 

tv  to  to  to  O 

•t  ru  to  o  ^3 

t-  nK  > 
•  *  *  r  C 


D-7 


ENTHALPY  DEPARTURE  CALCULATION  PROGRAM 


PRESSURE  AND  TEMPERATURE  GIVEN 


$IBFTC  CLARE  NOLI ST , NODECK 
READ (5,60)  N,L,WT,TC 
60  FORMAT ( 1H  ,213, 2F15 . 8) 

AL=L 

R=10. 73147 

CON=144 . 0/ (WT*778. 173) 

DO  80  K=1 , N 

READ (5 ,2)  X1,X2,X3,X4 ,X5,X6,X7,X8,X9,X10,X11,X12, 

1  X13 , X14 , X15 
2  FORMAT (1H  ,3E15.8) 

WRITE (6,52) 

52  FORMAT ( 1HJ, 18HEQUATION  CONSTANTS) 

WRITE (6,2)  X1,X2 ,X3,X4,X5,X6,X7,X8,X9,X10,X11,X12, 

1  X13 , X14 , X15 
WRITE (6,50) 

50  FORMAT ( 1H J, 6X, 7HP (PSIA) , 4X , 9HT ( DEG . R. ) ,1X,12HV(CU. 

1  FT./M. ) ,6X,7HV.CALC. , IX , 12HABS . ER.  IN  V,1X,12HPER. 

2  ER.  IN  V, 2X, 11HH (BTU. /LB . ) , 6X, 7HH .CALC. , IX, 12HABS . 

3  ER.  IN  H, IX, 12HPER.ER.  IN  H) 

SA=0  o  0 

SB=0 . 0 
SC=0. 0 
S  D=  0 . 0 
DO  20  1=1, L 
READ (5,1)  P , T , EXDH , VE 
1  FORMAT (1H  , 4F15 . 8 ) 

10  VO=VE 
DELV=0 . l*VO 

11  Vl=VO-DELV 
V2=VCH-DELV 
EA=EXP (-5.475  *T/TC) 

PO=R*T/VO+Xl/VO**2+X2/VO**3+X3/VO**4+X4/VO**5+X5/VO 

1  **6+X6*T/VO**2+X7*T/VO**3+X8*T/VO**4+X9*T/VO**5+X10 

2  *T/V0**6+X11*EA/V0**2+X12 *EA/V0**3+X13*EA/V0**4+X14 

3  *EA/VO**5+X15  *EA/VO**6 

Pl=R*T/VlfXl/Vl**2  +  X2/Vl**3+X3/Vl*M+X4/Vl**5+X5/Vl 

1  **6+X6*T/Vl**2+X7*T/Vl**3+X8*T/Vl**4+X9*T/Vl**5+ 

2  X10*T/V1**6+X11*EA/V1**2  +  X12  *EA/Vl**3+X13*EA/Vl** 

3  4+X14 *EA/V1 **5+Xl5 *EA/Vl **6 

P2  =  R*T/V2+Xl/V2  **2  +  X2/V2  **3+X3/V2  **4  +  X4/V2  **5+X5/V2 

1  **6+X6  *T/V2**2  +  X7*T/V2  **3+X8*T/V2  **4  +  X9*T/V2  **5+X10 

2  *T/V2  *  *  6+X1 1  *E A/V2  *  *  2+X1 2  *EA/V2  *  *  3+Xl  3  *EA/V2  *  *  4  -t  X14 

3  *EA/V2**5+X15*EA/V2**6 
EO-PO-P 

El-Pl-P 

E2-P2-P 


j;  ,t:  UO1.'. 

ct  :  (-••  .  :•  cia~ 

,£,i.v  ,  :  U  tv-hc 


..  •  •  • 

:  ,  ,  ^  ; 

-xx, i-  x, ex/. 

<  ,  .51  T,  HX)  TACO'S 
,  X  r  .  :■  "IT  . 

ATS YX  >  :  ■  :T.fi  T  "  IT  1  : 

c  l}  , .  xx  e*; 
rcG^sr/:.^ 

. 

■  :r  ...  ...  T  *  \.'- 

•/»dSC.xx,.:'JAv.i  '  • 

.  ,  H  : : :. 


0. 

C  ,OOl 

£  cx 

:v,  :ctx  (x  •?/'/•  .. 

OV’  i  ♦  0  -*VJ  a 

vVf  a-ov-  v 

.  v.  .  ^  .£~)  '  :•■  :■£  ••  : 

<  **  \  x\  i  -  •  \  ;  ■  -/  . :  •  :  ■ 

.  •  .  :  •  Vv'\T  3X  -  >  17\,T*'  K+X  r*.IV\T-  -'•:••••/' 

TV  ’•  i  *  !  £  I  .-i-.  XT'  ■  V  i •,  ■  ■■ 

/  *  1  V\AS  •  2  IX+  5  *  *  XV\A X: ;t  1  ’  /.+*  f 

r*SX+*** f;V\T* 8 X  >€ *  *  SV  \T* O, 
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D-8 


W=ABS (EO) 

X=ABS (El) 

Z=ABS  (E2 ) 

RR=EO/W 

S=E1/X 

A=E2/Z 

IF (W-l . 0E-02 )  22,23,23 
2  2  V=  VO 

GO  TO  15 

23  IF(X-1. 0E-02)  24,25,25 

24  V=V1 

GO  TO  15 

25  IF ( Z-l . 0E-02 )  26,27,27 

26  V=V2 

GO  TO  15 

27  IF (RR+S)  28,12,28 

12  DELV=0 . 5  *DELV 

VO=VO-DELV 
GO  TO  11 

28  IF ( RR+A)  29,16,29 

16  DELV=0 . 5  *DELV 

VO=VO+DELV 
GO  TO  11 

29  DELV=0 . 8  *VO 
GO  TO  11 

15  CONTINUE 

AT=CON/V 
AU=AT/V 
AV=AU/V 
AW=AV/V 
AX=AW/V 
EB=5 .4  75  *T/TC 
EC=EXP (-EB) 

ET= (2 . 0+EB) *EC 

EU=ET+EC 

EV=EU+EC 

EW=EV+EC 

EX=EW+EC 

DH=X1*2. 0*AT+X2 *1.5 *AU+X3 *4.0 *AV/3. 0+X4 *1.2  5  *AW+X5 

1  *6 . 0*AX/5. 0+X6*T*AT+X7*T*AU+X8*T*AV+X9*T*AW+X10* 

2  T*AX+Xll*ET*AT+Xl2*EU*AU/2. 0+X13 *EV*AV/3 . 0+X14*EW* 

3  *AW/4 . 0+X15  *EX*AX/5 . 0 
ERV=ABS  (VE-V) 

PE  RV=E  RV*  1 0  0 . 0  / VE 
ER=ABS (DH-EXDH) 

PER=ER*100. 0/ABS (EXDH) 

SA=SA+ERV 

SB=SB+PERV 

SC=SC+ER 

SD=SD+PER 


, 1  ,  ,  -  -V 

-7.  •  - 

t?\o: i  a  ■ 

<so  -?C0  . 

■ 

<v:  :  -v.  • 

sv-v 


\j  [2  7-  F/*  >V 
f: :  .  ,  s?  : 

r.'j*-  , 

3TJ  .  •  •' 


*<f  •  A+V.  '  •  •  ‘V;  +  .,rT*  ./  -0.2\XA  j  ‘  ■ 

(  -j- /)  aa a 

\o  :c  i*v*.r--va;~q 

(iiaxt;  ?  AQ  .ooi*HS«;t 
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2  0  WRITE (6,3) P,T,VE,V,ERV, PERV,EXDH,DH,ER, PER 

3  FORMAT ( 1H  , 10F13 . 4 ) 

SA=SA/AL 

SB=SB/AL 

SC=SC/AL 

SD=SD/AL 

80  WRITE (6,70)  SA,SB,SC,SD 

70  FORMAT (1HJ,16HAVE.ABS.ERJN  V= , F15 . 6/1X , 16HAVE . PER. 

1  ER.IN  V= , F15 . 6 /IX, 16HAVE . ABS .ER. IN  H=,F15.6/1X, 

2  16HAVE. PER. ER.IN  H=,F15.6) 

STOP 

END 
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